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Various problems arising in the theory of the excitation of a perfectly reflecting wedge 
or corner by a plane, evlindrical, or spherical wave, are dealt with. The incident wave is 
represented by a line source (acoustic or electromagnetic) parallel to the edge. The spherical 
wave is emitted by an acoustic point source or by a Hertz dipole with its axis parallel to the 
edge. The case of an incident plane wave field is obtained as the limiting case (for large dis- 
tances of the source from the edge) of the evlindrical or spherical wave excitation. 


1. Introduction 


Sommerfeld’s classical problem of the excitation of a perfectly reflecting half plane by a 
plane (acoustic or electromagnetic) wave field is based upon the theory of multivalued solutions 
of the wave equation [{34].2. Subsequently this method was extended to the problem of the dif- 
fraction of evlindrical or spherical waves by a half plane [3] or by a wedge [44] (in the latter 
case also for a plane incident wave). Here the incident cylindrical wave is emitted by a line 
source parallel to the edge (of the half plane or wedge) and the spherical wave is emitted by a 
point source (in the acoustic case) or by a Hertz dipole (in the electromagnetic case) parallel 


to the edge. The solutions for the half plane problem are given in the form of integral expres- 
sions involving a real variable, while the corresponding wedge problem solutions are given as 
contour integrals. 

The problems described above have been treated ever since with the aid of various techniques. 
The wedge problem as treated by MacDonald [22] requires the conventional methods of nine- 
teenth century mathematical physics only. Here the solution is given by an infinite series and 
this is transformed into Wiegrefe’s complex integral. For the special case of a half plane (angle 
of the wedge equal to 2 z), reductions to integrals involving a real variable are given. Another 
approach is a generalization of the method of the images [1, 16, 43]. 

At first the solution for a wedge of an angle +/m with positive integer m is established as 
the sum of the effects of the original source and its 2 m-1 images. This sum is replaced by a 
contour integral which is extended so as to obtain a formula valid for any positive m and hence- 
forth for any angle @ of the wedge. Again, the solution is given by a contour integral. In 
more recent times a number of attempts have been made to apply integral transform methods 
(20, 17] and integral equation methods [23, 14, 13] to the half plane diffraction problem. 
Finally, the problem of the diffraction by a perfectly conducting half plane of electromagnetic 
waves emitted by an arbitrarily oriented electric or magnetic dipole has also been solved [32, 
39, 15). 

Sommerfeld’s solution (plane waves incident on a half plane) is given in a form readily 
available for numerical computation inasmuch as it contains only the tabulated error fune- 
tion besides elementary functions. The solution of the other cases can either be given in the 
form of infinite series (usually suitable for the numerical evaluation of the near field) or in the 
form of integral expressions. A reduction of the latter to an asymptotic series has to be carried 
out if approximate expressions are desired for distances that are large compared to the wave 
length. Expressions of this form have been derived for the case of the diffraction of a plane 
wave by a wedge [44, 29]. (An asymptotic expansion which fails to give the proper behavior 
in the neighborhood of the boundaries of geometric optics is ignored here [36].) The main 
feature of these asymptotic series is that their terms are generalized error functions. The work 


referred to above is concerned only with time harmonic problems. 
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A second group of investigations centers around problems in which the dependence upon 
time is not periodic. This is the case of the diffraction of a pulse or a transient wave by a wedge 
or a half plane. This group of investigations was initiated by Sommerfeld [35]. The following 
contributions are listed here: Normal incidence of a plane pulse on a half plane [21,11]; two- 
dimensional pulse incident on a half plane [12]; plane pulse (unit step function) incident on a 
wedge or a corner [19,24, 25]; arbitrary pulse incident on a wedge or a corner [18]; evlindrica] 
Dirac pulse incident on a half plane [38]; and spherical pulse (unit step function) incident on a 
half plane [40]. This paper will deal with various problems arising in the theory of the exci- 
tation of a perfectly reflecting wedge or corner by a plane, cylindrical, or spherical wave field. 

The incident cylindrical wave is represented by an acoustic or electromagnetic line source 
parallel to the edge. The spherical wave is emitted by an acoustic point source or by a Hertz 
dipole with its axis parallel to the edge. The case of an incident plane wave field is obtained as 
the limiting case (for large distances of the source from the edge) of the cylindrical or spherical 
wave excitation. 

In section 2 of this paper, various representations for the time harmonic solution of the 
problems outlined before are given with special emphasis to a representation upon which the 
investigations of sections 3 to 6 are based. A straightforward method is applied based upon the 
representation of the incident field in the form of a Kantorovich-Lebedeyv transform. 

Section 3 contains a detailed discussion of the results obtained in section 2. Relations to 
geometric optics are discussed. 

Section 4 contains an asymptotic expansion for the far field excited by the incidence of a 
plane wave on a wedge. This expansion is not of the Wiegrefe-Pauli type inasmuch as it con- 
sists of two error function terms and an asymptotic series with inverse powers of the variable, 
Remarkable is the fact, that this series represents asymptotically (for large distances from the 
edge) the deviation of the wedge field from the corresponding half-plane field. 

Section 5 contains expressions for the energy radiated from a Hertz dipole in the presence 
of a perfectly reflecting wedge or cornér (antenna and corner reflector). 

Section 6 gives the pulse or transient solutions corresponding to the time harmonic prob- 
lems investigated in section 2. It is shown that previously known results can be obtained as 


special cases of the expressions given here. 
2. Time Harmonic Wedge Diffraction Problems 


The case of the excitation of a perfectly reflecting wedge by a (acoustic or electromagnetic) 
line source parallel to the edge of the wedge is considered first. (The solution for the excitation 
by a point source can then be obtained by a further integration.) Cylindrical coordinates 
(p,o,2) are introduced such that the 2 axis coincides with the edge of the wedge and the walls of 
the wedge are bounded by ¢=0 and ¢=a. Let the locations of the line source and the point 
of observation be characterized by ((p’,¢’) and P(p,¢), respectively (fig. 1). Let the field 
inside the wedge be periodic of the form 





b= (p,d)exp(iat), (1) 
pra 
Ppp) 
r 
Q(p, p) 
—> p-0 
FiGuRE 1. 
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where [ satisfies the two-dimensional wave equation 


, O.1W,10U , 
AU +k ; 5 sar thU=0. 2 
| a. dp tp? dee t (2) 


A solution of (2) representing divergent waves is of the form 


U exp (4 ivo) HH, (kp). (3) 


H,® (kp) is the second Hankel function, and y is an arbitrary parameter. Following a pro- 


cedure used elsewhere [26], the substitution 
k=22r/\¥=w/ce iy (4) 
is introduced. This means that instead of (2) and (1) 
AlU—yl=0, &=lU'(p,d)exp(yet) (5) 


is considered with a particular solution corresponding to (3) 


['—exp(+ivo)K, (yp), (6) 
where 
K, (z) = — }irexp(—i}vr)H,® (—iz) (7) 


denotes the modified Hankel function. The transition from k to y amounts to the transition 
from a wave problem to an exponential decay problem. At first the exponential decay prob- 
lem is solved, and the wave problem solution is established by returning from y to k. This 
will not only simplify questions with respect to convergence but the investigations in sections 
4 and 6 are based upon the exponential decay solution and not on the time harmonie solution. 
The incident field in P(p,¢) due to the line source ((p’,¢’) is given by 


, 27 - ° 9 , 9 
U‘= H,™ (kr) Kyi y|p?+p’2—2pp’ cos (¢— ¢’)|'/*}- (8) 
T 
Let the total field [” be represented as a sum of the incident field LU‘ and a reflected 
field (°, where UT” has no singularities within the wedge: 
U=U'+U". (9) 
In order to represent ("' in a form suitable for the problem considered here it is advisable 


to express the modified Hankel function on the r. h. s. of (8) not by means of its well-known 
addition theorem [9, p. 44] but by the equivalent expression ((3A) appendix): 


9 


, 27 - © 9} 9 
Uy Ky vlp?+ i 2 pp’ COS (¢—¢’)]'?} 
T 


dix ‘| K (yp) Ky yp") cosh [v(r—|o—¢’|) |dv, (10) 


which is of the form (6). Let the reflected field be represented in the same form: 


U'=4i9r | Ki, (¥p) Ky (¥p’) [fi (v) exp (vd) + fo(v) exp(— v9) |dv. (11) 
0 


The so far unknown functions f, and f, are to be determined by the boundary conditions. 
There are two cases of boundary conditions: (a) 7=0 and (b) 0U/0¢=0 on the walls 
o=0 and ¢=a of the wedge. With electromagnetic waves on a perfect conducting wedge, 
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case (a) refers to the electric field, case (b) refers to the magnetic field. With acoustic wayes 
case (a) refers to free wall, and case (b) refers to rigid wall. 
This leads to 


fi +h= —cosh[v(r— 9’ 
fiexp(va + frexp| — va) coshly(4r a+’ )| |?) 
and 
hh, ~sinh[v(r— 9’ )] 
fexp(va)—frexp va) =sinh[y(r— at ¢’)]. 13) 
for cases (a) and (b), respectively. It may be noted that the boundary conditions have been 


satisfied by putting the integrand equal to zero. This can be justified by verifying that the 
solution thus obtained satisfies all necessary conditions [31]. But it may also be remarked 
that lis represented as a Kontorovich-Lebedev transform (9, p. 75). After f and h, have been 
determined by (12) and (13) the solution is then, by (9), (10), and (11), 


ee oe , sinh (av) ; ; 
l 4in | kK »(Y¥p K P Yp’ : coshl{p a o @’ )| + cosh{y a Oo—o® )] dy. 14) 
Jo sinh (av) 
The upper and lower sign is valid for cases (a) and (b), respectively. The relation ((9A 


appendix) shows that the convergence of the integral (14) at its upper limit for a real and 
positive y is secured when 
Os + ¢’ S2a, 


and this region contains the whole region of the wedge. If one returns now to the wave problem 


replacing y by 7k according to (4) one arrives at the result 


. my , sinh (ay ; , 
l ‘| exp (xv) HT? kp) H? kp cosh lv la Qo QO . cosh Via QD 2 dy. 15) 


sinh (av) 


The behavior of the integrand at its upper limit leads by [9, p. 8S, form. 22] to the condition 


! 


, . 
rs QO QO <a T 


This restricts the validity of (15) to the case of a wedge of an angle @ larger than x with 
the point of observation P in the geometric shadow of the source ¢4 The expression (14 
shall now be transformed into an infinite series. For this purpose the elementary relation 
[2, p. 368} 


coshily a—o@ ~~ €,COS (76 / a) 
oe , 


: » > € i. @4=2 torae 
sinh av) a n= ni--- (nwa): 


valid for 0S@S 2a is used. Term by term integration vields 


— ‘ ,, cosh [v(a—s)] . l q y sinh (av) ,- — 
Ky (Vp) Ky. (Vp ; sinh mvidyp pe COS (nT a@) = A »(Yp A AY¥p)de. 
‘ sinh (ve) a < Jo w+ (nr/a)? 
The integral is known [10, vol. 2, p. 176]. (Note that this formula is not only valid for positive 
integers but for Re n>0). Hence 


"= » sinh (rv) ,- ‘ : l . . - 
J id + (naw/a)* A (Ye) Ki (Yp ly » wD ys a(Yp)K, rial Vp ) 
o — 


for p<_p’ and for p> p’ the same formula with p and p’ interchanged. Therefore, the field 77 can 
be written as 


. 2. nr i } - 
[ Sy e{ cos[ " (o—¢) [Feos| 8 co 6)” |} Lon 10) Kwa (16) 
@ n=0 a a ; 
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for p<p’ and for pp’ the same formula with p and p’ interchanged. The expression (16) is 
valid for arbitrarily complex y and if y is replaced by ik the solution for the corresponding wave 


problem hecomes 


l = > ie of cos | (6 6") | } cos| m “(¢ ' ) |}on alkp) HH  (kp’) (17) 
Q net c c 


for p< p’, and for p> p’ the same expression is used with p and p’ interchanged. The preceding 
expressions (16) and (17) represent rapidly convergent series if the arguments of the Bessel 
functions are small. The expressions (16) and (17) shall now be transferred back into another 
integral expression which is the basis of the investigations in sections 3 and 6. For this purpose 
the product /,,.:a(yp)Ayx/a(yp’) in (16) is replaced by an integral representation ((1A) appendix). 


This leads to the following expressions. 


. dje cos (nO lo) Tn eialV¥p) Ky eial¥p’) 


mS) €, cos (n16/a) K,|v(p?+ p’?—2pp’ cos r)'?] cos (nrr/a)dr 
A at 


us p > €, COS (n7O/a) sin (n/a) | CXP (-—-n7ws a) Ky [¥(p? t p” +-2pp’ cosh x) */2] dz. (18) 
0 J0 
The second sum vanishes when a=ar/m, m=—1, 2,3, .... The two sums on the r. h. s. 


of (18) are denoted by s; and sy. 
Interchanging the order of summation and integration one finds 


*r \ 
S » ” ,* , ” T T 
8, (2r lim K,l¥(p?+ p’?—2pp’ cos s)'?| >) e, { cos| (at 0) |+eos| n wo |} dz; 
Vv Ji n=O Qa a 


(19) 


ie F Ky |y(p? + p’*+2pp’ cosh r)'/?| 


+ . us : T 
- exp (—nrsr/a) { sinf (3+ 0 | + sin| x 0) | bar. (20) 
n=l Qa a 


The sum in the integrand of (20) can be summed [2, p. 203] and the result is 





S (27) Ky |y(p?+ p’?+2pp’ cosh s)'?] 


. Tv . T 
( sin | r—0) | sin (7-4 | 
a 
um ; dy. (21 
Tv Tv 
cosh (rr/a) - cos (7 | cosh (7?sja)—Ccos | (r- | 
a a 


s, as given by (19) ean be expressed as a sum of a finite number of terms. The result is (appen- 
lix) 


s;=ar SSK He + p’2?—2pp’ cos(2ar+8,]' \. (22) 


r=r 


The limits of the summation in (22) are 


x+0 x—@ — 
‘ | 2a } ra E } (23) 


rt6 — : , 
Where the well-known symbol | = ] means the largest positive or negative number which 
~@ 
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: Tr om , . 
is smaller or equal to—>~- The sum (22) is zero when 7, >r, and the first or last term of the 
x 


sum has to be halved when 7; or r, is an integer. With (18) to (23) the expression (16) for the 
field can also be written as 

l= F(o—¢’) + F(¢+-¢’), (24) 
with 


F(@) Qin . 2, the y\p* 1 p”* pp’ cos (Jai 1 6) |! 2 lor , Ko ly (p? t p” t 2pp’ cosh r)'/* 
r r e hi 


, T ; T 
sin (r—@) sin (7+) 
a Qa 
+ dr, (25) 
T T 
cosh (9r/a) cos (rw—@) cosh (mr/a) cos (aw -+-@) 
a a 


and r; and r, defined by (23). 


2.1. Spherical and Plane Field Excitation 


The formula [9, p. 95, formula 53] 


24 2%) —lexp[—v(r2+ 22)!2]—2e ‘| K, [r (2+ 2?) '2J eos (zed, (26) 
. it] 


oe Ss 


expresses the generation of a spherical wave function by means of a cylindrical wave function. 
To obtain the field (” corresponding to an incident spherical field multiply the r. h. s. of (8) 
and (25) by cos(zr), replace y by (7?++*)'? and integrate the resulting expression with respect to 
vfrom0to ©. This operation performed an (8) gives, using (26), 


Qin! Ko |r (¥?+- 0)! 2] eos (zr) de=i (r?+- 2?) Sexp[—v(r? 4+ 22) '2] (27) 
i.e.,an incident spherical field originating from a source point ((p’,¢’, 0). The same operation 
performed on (25) gives similar expressions. (For similar procedures, see [42, 5]). It follows 


then that in the case of an incident plane, cylindrical or spherical field the field 1” within the 
wedge of an angle a can be written in the form 


l'—Fi¢—?d’)+ Fi¢+¢’). (28) 


F(@) is expressed as the sum of two functions F, and F;: 


F'(6) = F\ (0) + F418), (29) 
where 
fF (6) ; > t(0 +- Dear), (30) 
r r 


(31) 


; 7 rT 
. sin (w—@) sin (w+ 6) 
: a a 
fF, (0) (2a) | I(x) + dy. 
Jo 1 T 
cosh (rr/a) —cos (a3 —@) cosh (rr/a) — cos (r +0) 
a a 


The representations (28) are also valid for complex values of y with Re y=0 and conse- 
quently for y=7k with Im ks0. The functions f and // in (30) and (31) depend on the char- 
acter of the excitation and are listed below. (The incident plane field case is obtained as the 
limiting case of the incident spherical field case for p’—>@ .) 
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Functions f and H for the exponential decay problem 


hio—¢’ H (x) 
(incident field) 


exp[ypcos(¢@ — o’)} exp(— ypcosh x) 
plane field 


Qin 'Koly Zim! Koly(p? + p’? + 2pp’cosh x)!?] 
(evlindrical field) 


I 2exp[— y(r2+ 22)!2 (p? + p’2+ 22+ 2pp’cosh xr)-!2 


spherical field) -exp[— y(p? + p+ 22+ 2pp’cosh x)! ]. 


Functions f and H for the wave problem 
fio’ H (x) 
incident field) 


exp[tkpcos(@ — ’)] exp(—tkpcosh x) 


plane wave 


Hj” (kr) Hs? (k(p? + p’2 + 2pp’cosh x)! 2] 


(evlindrical wave) 


(r2+ 22)-! 2 exp[—tk(r?+ 2)!” (p?+ p’2?+ 22+ 2pp’cosh x)~!? 
(spherical wave) -exp|— tk(p?+ p’? + 2pp’cosh «+ 22)"?]. 
2 2_ on ry y1/2 26 
r=|p*+p 2pp’ cos (¢—¢’ )|'”. (32) 


The function /(@—¢’) represents the analytic character of the incident plane, cylindrical, or 
spherical field. The upper or lower sign in (28) has to be taken according as an incident elee- 
tromagnetic or acoustic field is considered. The total field ” is now represented in three dif- 
ferent forms, (14), (16), (28) and (15), (17), (28) for the exponential decay and the wave prob- 
lem, respectively. Since in the following, (28) will be dealt with exclusively, the expressions 
for U of the form (14), (16), (15), and (17) corresponding to an incident plane or spherical 
field are not listed here. They are obtained in the same manner. The field as represented by 
(28) consists of the sum of two functions /; and F,. Each term r=0 of F, represents a field of 
the same character as the incident field (;-=0) but originating from one of the images of the 
source 4 with respect to the walls of the wedge. The number of the images and their location 
depend on 7; and r, and henceforth on the configuration. F, is given by an integral expression 
which vanishes when the angle of the wedge is of the form x/m, m=1,2,3,. . .. Itis easily 
seen that in this case 1’=F\(@—¢’) + F\(@F¢@’) is composed of the effect of the original source 
and its 2 m—1 images. In case of an incident plane field and an angle 22 (half plane), 7, can 
be reduced to the error integral (39). 


3. Discussion of F, and F, for Particular Cases 


It is obvious that the number of terms of the sum (30) can be made arbitrarily large if 
a is arbitrarily small. The case 0< a< 7 shall from now on be referred to as a corner. It will 
be shown now that for a wedge (r<a<2r) the sum (30) consists of at most one term. Two 
cases must be distinguished here: case a, a>x, 0<¢’<a—r, figure 2; and case b, a>z, 


a—nr<¢'< ha, figure 3. The restriction ¢< }a@ means no loss of generality. 
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These cases represent the two possible configurations. 

(asea. The direction @’ of the source { is such, that a geometric shadow region exists 
The boundary of the geometric reflection is ¢,=—x2x—@’ and the boundary of the geometric 
shadow is @=—2+ 9’. 

Cave b. The direction @’ of the source Y is such that no geometric shadow region exists. 
The boundaries of the geometric reflection with respect to the walls ¢=0 and ¢=a are ¢,= r—¢' 
and @)=2a—r—a’, respectively. 

Let the region 0S @Sa be subdivided into three regions I, Il, LLL bounded by the walls and 
the boundaries of geometric optics ¢=¢, and ¢=¢@. Thev are defined by: 

Region I. O<¢oS¢, 
Region II. >, SoS 
Region Ill. @S¢Sa 
and 
o—17—?¢', g=1r+¢’ case a, >) 
o=—17—¢', g&=2a—r—@ case b. 34) 


In these cases the limits 7,, 7, of (30) are easily found to be as listed below 


(‘ase b 
Region I II Ill I II lil 
r,.? 0.0 0.0 0.-—1 0.0 0.0 0.0 
02 0.0 0, | 0, | 0.0 0, l # l 


The upper or lower line for (7,72) is valid according as @=¢@—¢’ or 6=¢+¢'. With these 
results one obtains by (30) for both cases a and b, 


Fi (¢—¢') = F\(o+ ¢')=f(o—¢') Ff(@+¢’) valid in I, (: 


~~ 


h\(@¢— ¢') + Fi(¢+ 6’) =flo—¢’) valid in IT, 36) 
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while in region Ii] 
F\(¢—¢’) $ Fi\(o+-¢’)=0 case a (37) 
k\(o— 6’) = Fi (+ 6’) =f(o— 0’) F£f(+¢’—2a) case b. (38) 


Equations (35) to (38) together with (31) give the representation for the field LU” by (28) 
and (29) in agreement with the results obtained elsewhere [26]. The case of a source within 
a corner (0<a< mr) (fig. 4) and referred to as case c is treated the same way. At first r; and 
r, are determined by (23) and then F, by (30). Together with F, by (31), the field U" is known 
by (28). 


p:a p=Y, 






p=" 


epgyre™ P=, 
r 





Figure 4. 


Analogous to the cases a and b there exist two directions ¢, and @ in which a plane wave 
field incident on the corner emerges after multiple reflections on the walls. They are given in 
(50) and (51). 

The second function F, (@) given by the integral (31) can be given explicitly when a=2x 
half plane) and when the incident field is a plane field. Since this result will be needed later 
it is represented again although known [34]. Now, for a plane incident field, f(@—¢’) 


exp lyp cos (¢— ¢’ )] and /7(r) exp| yp cosh s). For a=2r, 


' l 
F(8) r' cos ( | cosh ( 7: r) (cosh #+-cos 6)~' exp (—yp cosh «)dz, 
by (31). The substitution cosh += 1+? gives 
’ l 1/2 I 1/2 I 
F,(0) a cos (54) exp (—vyp) t-'?(1+- cos 6+1)~! exp (—yopt) dt. 


This integral is known [10, vol. 1, p. 136]. Thus 


F’,(0) bcos (4 6) sec (4 6) Erfe [(2yp)'?, cos § @|| exp (yp cos 8), (39) 
where 
Erfe (z)=29°'? | exp (—®)dt. (40) 
Be 

Since the half plane case is represented by case a, as described in figure 1, the matching 
expression for /, in the three regions is known by (35) to (37), and the field U7 (28) excited by 
an incident plane field on a perfectly reflecting half plane can be written down explicitly. It 
can easily be verified that U7 can be represented by an expression which is valid over the whole 

range of the half plane observing that by (40) 


Erf (— z)=2—Erfe (2), (41) 
479794 3s » 35] 








and so [” becomes 
’=exp [yp cos (@—¢’ )|' 1— 3 Erfe [(2yp)'” cos (4¢—3¢’)| 
Eexp [yp cos (¢+¢’)|' 1— 4} Erfe [(2yp)'? cos(5 + 3¢’)}}, (42) 

for the exponential decay problem. The solution for the wave problem is obtained as 
U’=exp [ikp cos (¢@—¢’)|: 1— § Erfe [(2:kp)'” cos (4 @—3 9’)! 

=exp likp cos (@+¢’)| 1— 4 Erfe [(2ikp)'* cos (4 4+ 4¢")|}- (48 
(For tables of the error function for complex argument see [6].) 

In general case, F, consists of a sum of two integrals of the form 


F(0)=I(a,r—0)+ (a,x +8), (44) 
with 


Ta) (2a)~' Hr) sin (rp/e)| (cosh (r/a) —cos (rbia)|~ dr. 


pba 


Obviously the denominator of the integrand (45) is zero at its lower limit »=0 when 
¥=v,—=7+(¢+¢’)=2na, n=0,41,42,... . (46) 
But the numerator of the integrand vanishes likewise for the same value. It is easily shown 
27] that 
lim J(aw) = = FH (0) (47) 


and the behavior of (45) is the same for y+ 2na. It may be noted that (47) does not depend 
on the angle of the wedge. For any configuration is always 


0<¢o+¢'< 2a, 0< o—¢' | <a. (48) 
This condition restricts (46) to only two possible directions ¢, and @ which define the 
boundaries of geometric optics as outlined before. They are defined by 


, 


, 
case a: @=—7—@ ,(n=0); hH=t-+ @ ,(n=—0), 
ease b: ; T o’. (n=0): @&=—Za—z—oO.~7 (n 1) 
in agreement with (36). A choice of values for the parameter n other than given above is 
incompatible with the configurations under consideration. In order to treat case ¢ it is 


practical to subdivide this case according as the angle @ of the corner can be described as 


a=a/(2i-+e¢«), (=1,2,3,... casec, 
(49) 
a=-/(2l+-1+.e), /=0,1,2,... case ¢, 
where ¢ is a number between 0 and 1. Then the relation (46) leads to 
, 
Q, €a—@® 
case ¢; (50) 
¢ = Min [a(2 e)—¢’ .¢’ + ec] 
and 
’ , 
2, a(l ej @ 
case ¢, (51) 


$ 


Minfa(l t-€)—¢' .a(1 é)- ¢'| 
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and ¢ is connected with @ by (49). In case e=0, this leads to g,=g, and particularly, 
Q, = ¢’ when a= 2x/(2/) 
o=®=a—¢ when a=/(2/+1). 


An elementary expression for the upper limit of Z(a,y) can be given. Integration of (45) by 
parts gives 


—rl (ay) =H (r) arctan [cotan (4 m/e) tanh (3 rr/a)] 


0 


~- | H’(r) arctan [cotan (3 m/e) tanh (4a2r/a)\dz. 


Since, with the exception of the plane wave case, /7(7)->0 when z—->@ it becomes 


l(a) =r | H'(a) arctan [ cotan (Srv a)tanh( 5 xz/a )] dy, (52) 
and hence 
Taw) Sao H’ (x) | arctan | cotan ( : ry a) tanh (5 wr /c )] dy, 


and, after replacing tanh ($2r/a) by its upper limit 1, 
(aw) S H(O) arctan[cotan(4ay/a)| x7. (53) 


The equal sign holds when y has one of the values defined by (46) and the accuracy of the 
inequality (53) increases with decreasing a. Since, in general }ay/a will not lie between 0 


and z the ‘‘reduced” property 
dey/a=try/ainl,, !,=0,1,2,..., (54) 
is introduced and /, is such that }ay/a@ is between 0 and «Thus 


(aw) Se" (0) arctan[cotan(4m/a)] ‘ 
and finally 


I (a) <41/H(0)||1—y a). (55) 
4. Asymptotic Expansions 


The derivation of an asymptotic expansion of the field U’ for large distances from the edge 
of the wedge is equivalent with the reduction of the integral (31) to an asymptotic series. Com- 
paratively simple is the case of a plane field incidence (fig. 5). Then for the exponential decay 
problem 

H(r)=exp(—vyp cosh x). (56) 


7-\p yp:r p:p' 


ra 


NZ 


ii — Frovnn 5. 
J 
ra psa 


= y’ 


a 


\ , 
! 
| 
| 
| 
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In this case F(@) is represented as a Laplace transform with yp as parameter. Methods 
for the asymptotic development of integrals of this type are known [8, p. 45]. Since by (44) 
F,(@) is composed of two integrals of the form (45) it is sufficient to investigate 


I (ap) (2a) ~' sin (w/a) leosh (4 /a)— cos (rp/a)|~' exp (—yp cosh sr)dz, 57) 
. ai 
and yw has one of the values r+ (¢+¢’). The substitution 
cosh r,=1-+¢ (58) 


transforms (57) into 


. » 


12 
I (a) exp ( ve) | t <K(t) exp (—ypt)dt, (59) 


with 


} 


K(t)= (2a)7' sin (#W/a)(2+1)7'? | cosh [w/a log (1 + t+ (f+ 21)}| cosh (mp/a))~!. (60) 


Thus (59) is represented as a Laplace transform. For the following, the discussion shall 
be restricted to case a, as described in section 3 and fig. 2. Cases b and c¢, as deseribed in 
section 3, can be treated in the same manner. Now, A(¢) as given by (60) represents a regular 
function in the vicinity of t=0 provided that ¢+¢’ #-x (corresponding to y#0 by (46)). The 
singularities of A(?) lie on the negative axis and the smallest of their moduli is equal to 


Min [1 +cos (@4 ¢’)|. 


The Tavlor expansion of A(t) near the origin is 
‘ | 


A(t) >) «A, (ay)t", t< Min]l -cos(o+’)], (61) 
with 
A, (aw) (n!) ld" h(t) (/t” |, nm (HD) 


The first two coefficients are 
Ay (ayy (2a 12 2ceolan! Sry a 
A,(a,y) (Sa)~'2-'®cotan(4 ay a) /1+2(9/ a)? [sin($rp/a)]~*}. (0:3) 


The asymptotic expansion of J (ay) for large yp is obtained by inserting the series (61) 
into (59) and integrating term by term [8, p. 48, theorem 6]. The result is 


te 


, & (64) 


- 
IS 

~~ 

x 

! 


(ap) ~ —exp(—yp) 2A, (ea) B (n+ 4) (ye) "7, 
0 


The coefficients A, tend to infinity when ¢+¢’ tends to # (and hence y to 0) and the expan- 
sion (64) becomes useless. This is the case when the configuration is such that one of the 
boundaries of geometric optics is approached. In this case the point t=0 becomes a singular 
point of A(¢) and no expansion of the type (61) exists. In order to obtain an asymptotic ex- 
pansion which is valid over the whole region one proceeds as follows. The integrand of (57) 


at z=0 is equal to —(2a)™! 


cotan ($7~a@). In the vicinity of y=0 (the boundaries of geo- 
metric optics) 


(2a@)~'cotan(} my a) (rp) '+ Op) (65) 


and the singular part of (65) is independent of the angle a. Let the field (7 be now written 
in the form 
U=(U) eure t+[(U—(U) cnae] 66) 


and by (28), (29), and (35) to (37) 
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ds 


4) 


‘) 








U= (UJ an2e t+ Fi o—¢') —[F2O—#' Janae F {F(O+ 4) —[F2O+ Janae} (67) 
Now, (()a-0, is given by (42), and by (44) 
Fi(@+¢') —|Fi(O49' )lan2r=1(a,n—F 9’) —1(24,4—6F9') + (ae +O+9')—(U(24,4+6+9'). 
The difference 
liaxr—o +’) —1(24,47—OF 9’) 


: . us 
(2a) sin (x—o +9") | @ 
: (4r)~' cos (46+ 40’) | 
a COS (99> 2 exp (—Yp cosh rs) dr (68) 


T 7 cosh (4.7) —sin ($¢+ $¢’) 
. cosh (rr/a) —cos (r—O+®@ ) 
a 


is by (65) such that the integrand of (68) vanishes identically when @¢+¢’=2, and the 
asymptotic expansion of (68) is given as the difference of two series of the type (64) whose 
coefficients A,(a.7—oF¢’)—A,(24,7—OF9’) tend to zero when ¢+¢’ tends to x. Using (64) 


one arrives at the asymptotic expansion 


F.(o¢—¢’ [F',(6—’) la-ort | F2(6+06’) —[F2(64-’) Janos S(Yp) 
~—exp(—Yp) >> B,(yp)~"—",  —3r/2< arg ¥<3n/2. (69) 
2=8 
The coefficients B, are defined by 
B,=T(n+3)'A,(awr—o4+ 6’) —A,(24,2—64+ 06’) +A,(a,7+ o—9')—A,(24,4+ 6-9’) 
L.A, (ar—o—o’)—A,(24,n—6—¢')+A,(a,xr+6 +0’) —A,(22,7+04+ 9’)]} (70) 
With the A, given by (62), (63). It follows then from (67), (42), and (69) for the field excited 
by an incident plane field on a perfectly reflecting wedge of an angle a for case a (fig. 5) 
[’=exp lyp cos (¢—¢’)]{ 1— }Erfel(2yp)'? cos ($¢—}¢’)] 
Eexplyp cos (¢+ 6)! 1— }Erfe|(2yp)'” cos (46+ $¢")]) +S(yp), 
where S(yp) is asymptotically given by (69). It follows then for the wave problem 
1’ exp |ikp cos (@—9')|, 1 — 4 Erfe|(2ikp)'”? cos (4¢—}3¢’)] 


-exp likp cos (6+ ¢’)|{1— $Erfe[(2ikp)'? cos (46+ $¢")]} +S(ikp), (71) 


where for large kp 


S(ikp) ~ exp | ikp —im/4) » Bi" (kp) ~"—* (72) 
n= 
and these expressions are valid over the whole region of the wedge. The plus or minus sign 
has to be taken according as the acoustic or the electromagnetic case is considered. Obviously 
S represents the deviation of the wedge field from the half-plane field and the expansions (69) 
and (72) represent asymptotic expansions for this deviation term. If (yp)'?cos(4¢+ $¢’) and 
(kp)'*?cos(4¢ + 43¢’) are large, i. e., for the far field at large distances from the boundaries of 
geometric optics one can replace the error function by its asymptotic expansion [9, p. 147}, 


Erfe (z) ~~! exp (—2*) 3) (—1I)"T(n4+-3)2e-""',  —39/4< arg 2< 32/4, (73) 
n= 


and the asymptotic expansion of (” is of the same type as the one for S. On the boundary of 


the geometric reflection (¢+¢’=7), (71) becomes 
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U=exp [ikp cos (¢@—¢’)]{1—4 Erfe [(2ikp)'? cos (4¢—4¢’)|} = 4 exp (—ikp)+S(ikp). (74) 


On the boundary of the geometric shadow (¢—¢’=7), 


U=texp (—ikp)+ exp [(ikp cos (¢+-¢’)]-{1— 4 Erfe [(2ikp)'” cos (4¢+-4¢’)|} + SUkp). (75) 
Using (73) these results can be written as 


l’=exp [ikp cos (@¢@—¢’)| F 3 exp (—ikp)+Ol(kp)-'? |, (74a) 


U'=}exp (—ikp)+ Ol (kp) 3). (75a) 


These equations express the fact that the far field on the boundaries of geometric optics is 
independent of the angle of the wedge. 


5. Energy Radiated From a Hertz Dipole and From an Acoustic Point Source 
in the Presence of a Perfectly Reflecting Wedge or Corner 


Let the source Q(p’, ¢’, 0) be a Hertz dipole with its axis parallel to the edge. The fune- 
tions f(@) and //(r) are then 


£(0) = (2 + 2)-'2 exp[— ik? + 2)" I, (76) 


so ° 


r=(p*+p 2pp’ cos 6)'?, 


,* 


2 exp[—ik(p?+ p’*+-2pp’ cosh r+ 2*)'/?]. (77) 


H (x)= (p?+ p’*+-2pp’ cosh r+ 2*) 
The function U° given by (28) (taken with the minus sign) represents the vector potential 
A of the electromagnetic field. Its components are 
pe ee ae es 2 (78) 
With this relation the components -,, £,, /, of the electric field vector / are known and 
the radiated energy can be determined by Poyntings method. The application of this method 


is very difficult in view of the complexity of the expression for 7”. It is therefore advisable 
to use the so called E. M. F. method for the determination of the radiation resistance R of 
the dipole [37, p. 457] and [33, p. 134]. FP is given by 

R -(L/J)Re(E,)p. (79) 


Here L and J are the length and the current of the dipole respectively. Re(/,)p is the 
real part of the z component of the electric field at the location of the dipole Q(p’, ¢’, 0). It 
follows from (78) for the field, 

E.=—il ( ~ +l ): (80) 
Cis a constant (real > 0) and is of no importance here since the “reduced” radiation resistance 
R/Ro will be determined. This is the quotient of the radiation resistance ? for the complete 
system (dipole plus reflector) and the radiation resistance /?, of the dipole in free space (with- 
out the presence of the reflector). At first the contribution of F, (which represents a set of 
dipoles consisting of the original source and a number of images) to £, is considered. From 
(30) and (76) 


U; Fi (@ -’) Fi (@ +-") 


pr ) » j ° , ° € {2 > 9° 9 ® , ” 9 ) 
= 55 (r2+-2*) —! exp [(—tk(r2+4-2*)'|— 0 (r24- 2%) —! exp [—ik(r24+-22)'“I, (81) 
p=” g=4: 


r3=p*+ p’*—2pp’ cos (6—¢’ + 2ap) 
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4) 














r2=p*+ p’*—2pp’ cos (6+ ¢’ + 2aq) (82) 


ro-@ r™—o+q 
- a si me 2a 

r—o+9’ r™—o— 9’ 
1 a , 2 2a , 


It is then easily found that 


. 2] r : n2 m2 
Re | - i '( a hel 1 ) ] -— ( ke > atkd,) a p aiken) | (83) 
- D n=ny m= ™m, 


with 
d,=2p’ sin (an), d,,=2p’ sin (¢?’+ am) 


n,=—I[$x/al, nN =[}x/al, 


r+29’ r—2¢’ 
my -| = ? | Mos | = ? } (84) 
a 2a 


2 $ 99 


ao) =o"! sin ¢£+¢°% cos f¢—F"3 sin ¢,  -g(0) = 2/3. 


The term »=0 in the sum (83) represents the contribution of the original dipole source. 
It follows then for the radiation resistance of the Hertz dipole in free space by (83) and (79) 


Ry =2L/(BI)CKk. (85) 
The contribution of F, to U" is 
(',= F,(@—¢’) — F:(o+ ¢’) (86) 


and it follows by (31) and (77) 


; 2/ ; . Ma 2kp’ cosh x/2) 
Re| ic : +h ») ] (2a) "Ch 4 2 sin (r*/a) | —.- : le 
Db 0 


4 ( 
cosh (mr/a)—cos (4?/a) 


: T — . T 94/ 
: sin (r—2¢') sin (r+2¢') 
di y Qa - 
< q{ 2k ‘osh 5 ) + dx >* (87) 
0 2 T F 
. cosh (ra) —cos | (r—2o | 
a 


lt follows, then, for the quotient of the radiation resistance P? of the dipole in the presence 
of the reflector and the radiation resistance /?, of the dipole in free space by eq (79), (83), (87), 


and (S85). 


1 
a 
= 
Fe 

~ 
=~ 
.~) 
~ 


T ») , 
cosh (rr/a)—cos (r+ 29’) 
a 


z 
: m e : qd ( 2kp’ cosh = ) 

R/R,=3/2) >> a(kd,) — >* g(kd,,) |—= <4 2sin (4/a) —~—— dr 

_ ’ i — 4 , la » cosh (rr/a)—cos (r?/a) 


” m m 


° T . T 
. sin | (3x— 26) | sin | (r+ 29) | 
: r a a 
— | a( 2kp cosh ) + : : ; = \|dxr * (88) 
vi ; 2 T ‘ , ul i_« , 
cosh (mr/a) —cos (r—2¢’) cosh (arr/a)—cos | —- (r+-2¢’) 
a a 
The occurring parameters and g(¢) are given in (84). Since g(¢) vanishes for large ¢ it is 
clear that if the distance of the dipole from the walls is large (p’-> ©) only the term n=0 in the 
sum gives a contribution, and then 


lim R/Ry=3/2g(0)=1. 
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r, i. e., in the case of.a wedge 


The formula (S88) is valid for an arbitrary a. But if a 
By (84) the 


reflector (figs. 2 and 3), the two sums on the r. h. s. of (S8) can be given explicitly. 


limits ”,, 2, m,, My, of the sums in case a>-z are 


ny =Ng=mM,=— M,—0 for 0< ¢’ ir, 
ny =N—m,=—0, mM, | for $r<¢’< ha. 
Thus 
SS) gtkd,)— >> g(kd,,) =2/3—q(2kp’ sin ¢’), if 0<¢’< $n, 
ii " m m 
2/33, if In<¢’< }a. 


The analysis is exactly the same when the source Q(p’,@,0) is an acoustic point source, 


In this case the plus sign in (28) has to be taken. The function g(¢) in (84) which character- 


izes the coupling term between two Hertz dipoles (parallel to each other and the line con- 

necting their centers perpendicular to their axis) has to be replaced by the corresponding 

coupling term 4,(¢) for two acoustic point sources [30]. 

n(o) =o sin &. (89) 

In this case the corresponding formula for the relative radiation resistance as defined 
before is 

4 , d 
. nf 2kp’ cosh =) 


n> m) 
R/Ryo=>5 alkdn)+ SS gtkd,) — (2a) -!< 2 sin (#/a) » dr 
F j 0 cosh (ar/a)—cos (2a) 


an m m 
° T A ° T 
sin (r—2¢’ ) sin (xr+-2¢’) 
a a 
ly 


T 7 


,. r 
+ orf 2kp’ cosh = ) 
0 : - ; T 9,’ us ‘ , 
cosh (mr/a)—cos \7—<zp ) cosh (ar/a)—cos (x-+-2¢') 
a a 


“ 


(90) 


The definition of the n,, ne, m,, m., d,, d, is the same as in (84), while g,(¢) is given by 


(89). If a>-x, then as before 


n; = No—=m,—mMm,=—0 for 0 <@'’< $x, 
nN, =No—m,=—0, Ms | for } 3: ¢’ : Sa. 
and 
ne ms 
‘ \ ‘ ; ° , ; 
S gi(kd,) + >) a(kd,,)=1+ gq, (2kp’ sin $’), 0< @ } mr, 
n=n m= m4 
, , 
l ir c+) >a 


6. Transient Solutions 


The results obtained in section 2 shall now be used to investigate the corresponding 
The majority of the contributions quoted at the beginning of this paper 


pulse problems. 
Integral transform 


are based upon a direct attack on the time dependent wave equation. 
methods were used in [38] (Kontorovich-Lebedev transform) and [40] (Fourier transform). 
In the cases considered here the application of Laplace’s transform leads rapidly to the de- 
Let the excitation of the source be represented by a pulse starting at the time 


sired results. 
Let the pulse function be such 


t=0. Let the pulse function be g(t), and g(t)=0 for t<0. 
that it admits its representation by Laplace’s integral formula: 
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q(t)= (21) '| | | g(r) exp ( vr)dr exp (thd. (91) 


In section 2, various expressions l’ for the exponential decay problem were given. This 


case corresponded according to (5) to the time dependency exp (yt), (put e=1). It can then 
be concluded that for a pulse excitation g(f) as represented by (91) the expression for the field 


is represented by 


&(t)—(2ri)~! | | g(r) exp (—Yr) dr |U exp (yt)dy- (92) 


The solution of the pulse problem is therefore given as the inverse Laplace transforma- 
tion of 


v | g(r) exp (—Yr)dr, (93) 


with respect to y. 
6.1. Special Cases for g(t) 


For a unit step function g(t) =1, the field (92) becomes 


& (¢) = (27) | l'y~' exp (vt)d¥. (94) 


For a “Dirae”’ pulse excitation imposed at the time t=, >0: 
a(t) =6(t—t,), one finds first 


| g(r) exp (—Yr)dr | 5(r—?t,) exp (—Y¥r)dr=exp (—Y1t)). (95) 
The field due to a “Dirae” pulse excitation becomes then 
%*c4 i= r 
®,,(t—?t,) = (29) | U” exp [y(t—t)] dv. (96) 


If the order of integration in (92) can be interchanged, it results in 


(1) | ain { Qn) ‘| U exp [v(t— liv bar. (97) 


Since by (96) the inner integral is the effeet of a ‘Dirac’ pulse, excitation (97) can be 


written as 


P(t) | Y(t) Pp(t—r)dr, (98) 


i. ¢., the field due to an excitation g(t) of the source Q can be obtained by an integration proc- 
ess involving the field due to a “Dirac”? pulse excitation. The general equations (96), (98) 
shall now be applied to the problem of the diffraction of arbitrary pulses (plane, cylindrical, 
and spherical) incident on a perfectly reflecting wedge or corner. The case of “Dirac’’ pulse 
excitation imposed at the time t=0 shall be considered first. The Laplace inversion of a 
function h(y) with respect to y shall be denoted by A(#), i. e., 


h(t) = (2ri) “| h exp (vt)d¥. (99) 


It follows from (96) and (28) for the field &, due to a plane, cylindrical or spherical “Dirac” 


pulse ; 


359 








®,(t)=F(o—¢') =F (o+o’) 


>> f(o—¢'+2an) = SS f(o—e'+2am)+ H (x) |G(o—’) =G(o+ o’) dz, 


” n m my 


™—-Q--@ T—-Q—-@® TQ QD ™—-QO—@Q 
thy — b. Vi . Mi — . Mle a 
2a : 2a 2a 2a 
T . Tr , 
sin (r—@) sin (3-0) 
a a 
+ 


(7 (A) - (2a) 
r © 
cosh (rr/a)—cos (r—@) cosh (ar/a)—cos (r+) 
a a 


The respective Laplace inversions f(@) and H(+) are for a plane pulse: 


where 


£(0) =6(t+ p cos @) 
H1(.2) =5(t—p cosh vr). 
For a evlindrical pulse one has 


{(0) =2i9-'(f — p*— p’*+2pp’ cos 9) '?, for t>(p*+ p’*—2 9’ cos 8)'”, 
0, otherwise. 


— 
oO 


H(2) =2ix-' (fF — p?— p’?—2 pp’ cosh s)~', for 2pp’ cosh r< tf — p?— p” 
0, otherwise. 


t—(p*+ p+ 2°+2pp’cosh x)'?]. 





Hx) =(p*? + p+ ~+2pp’cosh 2) 


(100) 


(101) 


(102) 


(1035) 


(104) 


(105) 


(106) 


(107) 


The Laplace inversions for /(@) and //(2) are obtained from (95) for the plane and for the 
spherical pulse, while [10, vol. 1, p. 277, formula 8] gives the inversion for the cylindrical 


pulse. It does not seem possible to express the integral | H(x)G(@)dz in a closed form in the 
0 


case of a cylindrical pulse. To obtain the field due to a pulse excitation g(t), the equations 
(98) and (100) have to be used. For a plane and a spherical pulse the integration process in 


(98) can be performed immediately. The results are 


ne m 


&(t)= >) glt+pcos(o—¢' + 2an)| — >) glt+p cos (6+¢’+2am)|] 
a=" mt m 
*cosh tp 
- g(t—p cosh r)(G(o--0') =G(o+ ¢’)|dz, 


for the plane pulse. 


©(t)=>5 R;'gt—R,) F DS Ra'gt—h,) 


— 
n nn m m 
"y 
+{ (p?+ p’2+- 22+ 2pp’ cosh r) ~'? g[t— (p? + p’?2+- 22+ 2pp’ cosh x)'/?|] { G(@—¢’) F G(e4 
0 


for the spherical pulse. 
The upper limit z, of the integral is given by 
2pp’cosh z,=—C—p*— p’”*— 2 for (109), and cosh x,=t/p for (108). 
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(108 


(109) 


(110) 














RR, =|p° + p’*+-27—2pp’ cos (¢—¢’ +2an)]'” 
(111) 
100 R,=[p?+ p+ 27—2pp’ cos (6+ $9’ +2am))"”. 
} 

The other properties are given in (101). Again in the case of a wedge (a>2), the prop- 
erties mm), Me, mm), Mm. are given in section 3. The representation of the cylindrical pulse field 
leads to a double integral. For the special case of a half plane, however, the integral occurring 
in (100) can be solved. In the following some special cases are considered in which the solution 

) for the transient field become elementary. 
01 6.2. Special Cases 
) 
a. Cylindrical Pulse and Half Plane 

The case of a line source Q parallel to the edge of a perfectly reflecting half plane is con- 

sidered. The regions I, I, III as deseribed in section 3 are indicated in figure 6. For a half 
(2) plane (a Ir), one obtains from (101), 
(3) G(0) = — 2"! cos (46) cosh (42) (cosh 2+ cos8)~!. (112) 
It follows then for a evlindrical “Dirac” pulse bv (105). 
4) aor ee) 
Pa 
il 
5) X 
I y 
0 i. 
6) bis rae 
= +p Q (p,-y") 
i) FiGure 6. 
he "? ae ' — "I (t?— p*?— p’"*—pp’ cosh «)~'? 
al (4) GQ(0)dr=—2ix-? cos (46) cosh (4.) dr. 
0 a . cosh r+cos 6 
he ros . . » . . . , 9 » so ~ 2 
lhe upper limit 2, of the integral is given by 2pp’ cosh 2,=—p?—p”, and from (105) it 
ns is obvious that the integral is only different from zero when t>p+p’. The substitution cosh 
in r=1+2*, cosh (4r)dr=2' "dp gives for the integral 
Lin\G@de ix~' cos (40); sec (46)|(t2—p?—p’?+2pp’ cos 0)~'”, if t>p+p’, 
0, otherwise, (113) 
: It follows then from (100) by means of (35), (36), (37), (104), and (113) for the field due 
to a cylindrical “Dirae’’ pulse incident on a half plane 
®,(t)= ix! { (P— Rj)? + (P— R3)-'2@— 3 [ (@— R})-'* F (P— R3)-'7)} ~—sfor region I, 
\ 2ia—'{ (P®— R})-'4— 3 (P@ — R?)-"47 +. (A@— R32)" J} for region IT, 
| 
ix '[((P@— Ri) 2 = (P@— R3)-"], for region III, (114) 
) where 
Ri= p?+-p’?—2pp’cos(¢—¢’), Ri=p?+ p’*—2pp’cos(o+ ¢’). (115) 


In (114), (@—Ri,)-2=0 for t<p+p’ when this expression is within the brackets; other- 
wise zero for t< Rj,.. The expression (114) is elementary in all the three regions. Former 
investigations [38] established @),(¢) as an infinite series involving Legendre functions. (It 
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can be shown that this series can be summed and leads to the expressions in (114).) From 
(114) and (98) it follows for the field @ due to an excitation g(f), 


"se Ry *? Ry 
(t)=2ir {| [(¢—7)?— Ri] -' g(r) dr [(¢—7)?— R53]! g(r)dr 


°*t—p—p’ 
4] \((t—r)? I?) '2=((t—r)? I?) "alg(sidr \ region I. 
J 0 


"6 Ry 
(t)=2ix { | \(t—7r)?#— Ri]? q(r)dr 


°¢—p»—p’ 
—4{ [((¢—7)?— Rj) -'?#+ ((t—7)?#— F3) 'aly(aide be region II. 


gion IIT. (116) 


°g=—p—9/ 
b)t)=i9 ‘| [((t—7)?— RF) - 2 = ((t— 7)? — R3) —"2 J] g(r) dr, r 


b. Plane Unit Step Function Pulse and Wedge or Corner 


Now, g(t)=1 for t>0, and g(f)=0 for t< 0. Then the terms of the sums on the r. h. s, 
of (108) are either zero or one according as the argument of g is negative or positive. The 


well-known formula 


. : i ] j 1/2 
leosh (62)—a]~‘dr=26-'(1—a?)~'? aretan ( : ) tanh (3b. | 
: l—a s 


cives for the integral in (108) by (101 
*cosh tip ‘ ; “/- (?? 2) 1/2) F/a Tia T 
| G(@dr——-x7' arctan . a I A cotan| . (r—@ ] 
0 er (&—p) +e <a 
t+ (#2 2) 1/2] */a F/a z 
—r arctan | = A : I A cotan (a4 +0) 
[¢+ (t?— p*)'4|*/-+ p*/* 2a 


for t>p and equal to zero for t< p. Thus the field (108) for a unit step function excitation 
| | 
(19, 24] is elementary. 
c. Spherical Unit Step Function Pulse and Half Plane 


The center of the spherical pulse is Q (fig. 6). Again, g(t)—1. The excited field is 
given by (109) and G(@) is given by (112). The integral occuring in 109) is of the form 


| (p* | p’? J 2pp cosh r) '?@G(0) da 


—r ' cos (30) | cosh (4.7) (cosh 2+-cos @)~'(p?+ p’?+ 27+ 2pp’ cosh «) 7! 2dr 
—2'?7-' cos (30) | (1+cos 0+-2*) —'[(p+ p’)?+ 22+ 2pp’r*] -!2 de, 


using the substitution cosh r=1+ 7, cosh ($4)dr—2'? de. The upper limit 7; is given by 
r= |P—(p+p’)?—2]'. Again this integral is elementary. Each sum in (109) reduces to at 
most one elementary term (see the tables in section 3), and the resulting field [40] is elementary. 

It should be remembered that the previous investigations were based on the assumption 
that the velocity of the field is unity. The time parameter ¢ has therefore to be replaced by 
ct when the velocity of the field is equal toc. The expression U’ for the time harmonic field as 
given in section 2 can easily be regained as the limiting case of a transient field ® due to an 
excitation g(t)=exp (lot). 
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8. Appendix 


To prove 


rl,(a) K,(b) | K,| (a2+ 6?—2ab cos t)'?] cos (vt) dt—sin (xv) 


exp (—vt) Ky|(a?+ 6?+2ab cosh t)'*) dt, ab, (Al) 


K, (a) K, (6) K,| (a?+ 6?+ 2ab cosh ¢)'*] cosh (vt) dt, (A2) 


» 


Kk, | (a?+ 6? —2ab cos ¢)'?]=227! K,,(a)K;,,(b) cosh [v(r—'@ | dv, OS ReoS2e,  (A3) 
. ti 


one uses the formulas 


rly(z) | exp (2 cos f) cos (vt) dt—sin (rv) exp (—z cosh t—vt) dt {9, p. 21], (A4) 
K, (z) aI s- ' exp(—2—! t/a dx 19, p. S2] (A5) 
v - 2 q - e | 4 ~ dy ls } + Owls ake 
' 2+ h? 


) I, ( : ab/r )de, (9, p. 53}. (A6) 


- ] a- 
1, (a) Ky, (6) 5! q ‘exp(- r— 4y 


8.1. Proof of (Al) 


Starting with (A6) and replacing the modified Bessel function in the integrand by (4), 
one obtains, interchanging the order of integration, 


» yw . 2 h?—2ab cos 
rl, (a) K,(6) >|, cos on! [ r ‘ exp/( :— —— ‘Yate Jat 


’ P 24 f2i 2 be s| / 
— sin (re) | exp (—vt) | | "ie exp ( nn i —_ Jar ar, 
~ 0 0 Jf 


The inner integrals on the r. h. s. are known by (A5). This proves (Al 


8.2. Proof of (A2) 
The relation 
K,(2) = 4 sec(av)[I_,(- 1,(2)], (A7) 


together with (Al) gives (A2). 


8.3. Proof of (A3) 


Replace in (A2) » by i and obtain by Fourier’s inversion theorem 


Kyl (a +. §24 ab cosh ft)! 2] Qr ‘| iad a) K,,(b) cos (vt) dy. (AS) 
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4), 
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(7). 
gen 


4) 








Replace here ¢ by i (#—@) and obtain (A3). ‘The restriction 0 ¢Re ¢<2x is necessary to 
ensure convergence of the integral in (A3) at its upper limit, since for fixed real z and large 


positive » >< (9, p. 88, formula 19}, 


l , 
K.,(z) = (22)'?("?—2*) -" exp ( —br){ sin[ b cosh~!(»/r) — (v?—2")! +7 ]+00 y}. (A9) 


8.4. Evaluation of s, Defined in (19) 


Using the formula 


as sin |(2.V—+1)y/2| 
€, COS (ny) .; i, 
n= . sinh \y é) 


one obtains instead of (19) 


s,— (22) ~' lim K\¥(p?+ p’?—2pp’ cos r)' | 
V+ Jo 
: T : , T : 
sin} — (2N+1)(4+@) | sin} — (2N+1) (4—8@) 
2a Za 
+ dy. 
° T ° Tv 
sin (r+ 6) sin (r—6) 
2a 2a 
Substitute in the first integral #(2+@) 2at and in the second x(r4—@)=2at. It follows 


then that 
sin (2N+1)¢ 


2 ; dt, 
sin ¢ 


x; =ar*lim | Ky {y(o?+ 0’2—2pp’ cos (6+ 2at/x)|}! 
Ne of 


T Tr . ° a ° ° 
where f, - (+6) and t,=5- (r—8). Now s, is a Dirichlet integral and its value follows 


~a =a 
immediately [2, p. 492] as represented in (22). 


Wasnineron, D. C., January 14, 1958. 
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Vo!. 61, No. 5, November 1958 Research Paper 2907 


A Liquid-Helium Cold Cell for Use With an X-ray 
Diffractometer 
Igor A. Black,’ Leonard H. Bolz, Frank P. Brooks,' Floyd A. Mauer, and H. Steffen Peiser 


\ liquid-helium cold cell for use with the General Electric X-ray goniometer has been 


designed and constructed. 


It is used to obtain X-ray diffraction patterns of polverystalline 
solids deposited from the gaseous phase on a helium-cooled surface. 


This work is part of 


the program for studying the stabilization of free radicals. 


1. Introduction 


{ 3-vr program of basic research has been under- 
taken at the National Bureau of Standards to study 
the preparation of free radicals and the preserva- 
tion of these unstable molecular fragments by freez- 
ing. The usual technique involves dissociating a 
gas in a furnace, an rf field, or a spark discharge, and 
depositing the products on a surface cooled by liquid 
helium. 

The freezing technique is widely used to store 
molecular fragments for study by spectroscopy, 
calorimetry, paramagnetic resonance, and other 
techniques. A knowledge of the structures of the 
solids deposited and the transformations they under- 
go during warmup is important in interpreting the 
results of these experiments. In calorimetry, for 
example, it may be difficult to determine whether an 
observed heat effect should be attributed to re- 
combination of free radieals or crystallization of an 
amorphous deposit. To aid in settling such ques- 
tions, a study of the structures and transformations 
of solids encountered in free-radical research was 
initiated. Special equipment was used for freezing 
gaseous samples at liquid-helium temperature and 
obtaining X-ray powder diffraction patterns. This 
paper describes the construction and use of the 
liquid-helium cold cell or’ cryostat developed for 
this work. 


2. Selection of Techniques and Equipment 


Reviews of the literature on low-temperature X- 
ray diffraction by Ruhemann [1]? and by Steward [2] 
as well as the additional references at the end of 
this paper [8 to 8] show that a considerable number 
of X-ray cryostats have been described. Most of 
these, however, are for use with film cameras and 
are not well adapted to observing transient effects 
that occur on warmup of samples containing trapped 
free radicals. Of the diffractometer cryostats de- 
scribed, only the one by Barrett [7] can be used at 
liquid-helium temperature, and it is designed to fit 
X-ray equipment. that is no longer available. 

In order to take advantage of the rapid-recording 
feature of modern electronic-counter ditfractometers, 
design of a cryostat that could be substituted for 


the room-temperature specimen mount on one of | 


a 
Arthur 1). Little, Ine 
' 


Figures in brackets indicate the literature reference at the end of this paper. 


these instruments was undertaken. The X-ray 
equipment selected for modification was the General 
Electric diffractometer, the advantage of this in- 
strument being that it has a vertical goniometer 
axis. The cryostat can be made to turn on this axis 
as required without causing liquid refrigerants to 
shift position. 


3. Design and Construction 


Onee the X-ray technique and equipment had 
been selected, design of the cryostat followed along 
fairly conventional lines. It may be considered to 
have evolved from the optical cold cell described by 
Duerig and Mador [9]. X-ray windows and other 
special fittings had to be provided, and many changes 
were made to take advantage of modern manufac- 
turing techniques, but the basic design is quite 
similar. 

In constructing the cryostat, stainless steel was 
used wherever possible. Joints that do not have to 
be demountable were “heli-are’’? welded and leak 
tested before assembly. This type of construction 
complicates modification or repair of the apparatus 
but almost completely eliminates troubles from 
vacuum leaks. Of the demountable joints, only two 
gasket seals and one soft-solder joint have to be 
vacuum tight. 

The X-ray specimen cryostat is shown (in section) 
in figure 1. It is made up of two concentric Dewar 
vessels. Liquid helium is contaitied in the inner 
vessel (8) which has a hemicylindrical copper block 
(17) at the bottom. The specimen deposited on the 
plane, vertical surface of this block replaces the 
usual room temperature specimen. Beryllium win- 
dows (19) in the outer wall permit the X-ray beam 
to enter and leave the cryostat so that diffraction 
patterns of the material deposited on the surface can 
be obtained. 

The helium vessel is surrounded by a liquid-nitro- 
gen tank in the shape of a hollow cylinder (9) with 
copper radiation shields, (7) and (15), attached to it. 
The X-ray window in the radiation shield is covered 
with 0.00035-in.-thick nickel foil (18) which shields 
the specimen from room-temperature surfaces and 
filters the X-ray beam to eliminate copper K@ radia- 
tion. Except for small openings for the specimen 
tube and optical windows, the helium vessel is 
entirely surrounded by surfaces which are at or near 
liquid-nitrogen temperature. 

Vacuum spaces around the helium-cooled specimen 
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Figure 1. Section through the \-ray-diffraction cryostat. 


1, Valve for evacuating helium Dewar; 2, helium fill-tube stopper; 5, helium 
vent: 4, electrical connector for helium-level indicator; 5, nitrogen vent or fill tube 
one not shown); 6, helium fill tube; 7, upper radiation shield; 8, helium vessel, 
9, nitrogen vessel; 10, valve for evacuating specimen chamber and nitrogen Dewar; 
11, resistors for helium-level indicator (2); 12, thermocouple leads; 13, demount 
ible flange; 14, inner-vacuum wall; 15, lower radiation shield; 16, brass ring and 
soft-solder joint; 17, copper block; 18, nickel foil X-ray windows 0.00035 in. thick; 
19, beryllium X-ray windows; 20, thermocouple junction; 21, translational-ad 
justment vernier; 22, translational-adjustment slide; 23, rotational-adjust ment 
worm gear 


block and around the outside of the nitrogen con- 
tainer are connected and are evacuated through a 
single valve (10) at the back of the apparatus. 
Another valve (1) is used to evacuate the space 
around the helium container. This space is sep- 
arated from the specimen chamber by a wall (14) 
and seldom requires pumping. The wall probably 





accounts for as much as 75 percent of the heat 
transferred to the helium vessel once the specimen 
has been deposited, and reduces the operating time 
available on one filling to about 2 hr. However. 
separating the vacuum space around the helium 
vessel from that around the specimen prevents the 
sudden loss of helium and specimen that would resylj 
if the pressure in the specimen chamber should rige 
excessively during deposition or warmup of the 
sample. 

X-ray windows (19) in the outer wall of the ery. 
stat are covered with beryllium strips 0.025 in, thick. 
The two windows cover 85° of are each and permit 
measurements to be made in the angular range from 

2° to + 168° 26. The beryllium strips are held jn 
place with wire bands tightened by screws and are 
sealed to the stainless steel walls with Glyptal. The 
windows have not leaked during 6 months of use and 
have even survived a small azide explosion inside 
the crvostat. 

A thermocouple (20) is used to measure the tem- 
perature of the copper block. For the alloys used 
(Ag plus 0.37 atomic percent Au versus Au plus 
2.1 atomic percent Co) the sensitivity is approxi- 
mately 10 wv/°K at 4.2 °K. The junction is inserted 
in a well at the bottom edge of the specimen surface 
and the leads are cemented to the copper block te 
keep them at the temperature of the junction. The 
presence of gas in the specimen chamber affects heat. 
transfer conditions, and if this precaution is not 
taken, temperature readings are pressure dependent, 

The resistance of two carbon resistors (11) mounted 
in the helium vessel is recorded to indicate when 
liquid helium reaches the \-liter and the 2-liter 
levels. A current of about 50 ma flows through the 
two 15-ohm resistors and the rate at which heat is 
carried away (by the gas in one case or the liquid in 
the other) determines the temperature and, hence, 
the resistance. 

The outputs of the thermocouple and the liquid- 
level gage are recorded continuously on a chart 
that is synchronized with the diffractometer chart. 

Pressure Measurements are made with a Philips 
cold-cathode ionization gage. Although the modi- 
fication is not shown in the figures in this paper, the 
gage tube has recently been mounted on the cryostat 
itself so that pressure measurements can be made 
after the pumping manifold has been detached. 

A slide (22) is provided for moving the cryostat on 
its base in a direction normal to the specimen surface 
in order to aline that surface with the axis of the 
goniometer. A worm gear (23) on the base is used 
to rotate the cryostat around a vertical axis to make 
the normal to the specimen surface bisect the angle 
between the incident and diffracted beams. 

The dimensions have been chosen so that the slide 
can be removed from the base shown in figure | 
and inserted in another base that fits the North 
American Philips high-angle goniometer. However, 
that instrument has to be modified to operate with 
its axis vertical. 

The lower chamber can be detached at the flange 
(13) to make the specimen block, thermocouple and 
lower radiation shield accessible for cleaning and 
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repairs. The flange is also intended to permit the 
lower chamber to be replaced with u film cassette. 
if modification of the copper block is required, it can 
be removed by unsoldering the soft-solder joint (16). 
These provisions are intended to make it possible 
1o modify the cryostat (ec. g., for single-crystal work) 
if necessary. 

~The cryostat is shown mounted on the X-ray- 
diffraction unit in figure 2. X-rays from the source 
Jit (1) enter the eryostat through the beryllium 
windows (2). The diffracted beam enters the de- 
‘rector slit system (3). Its intensity is plotted by 
the recorder (4) against the diffraction angle, indi- 
sated on the protractor (5). Gaseous samples are 
introduced through the silica tube (6). The speci- 
men surface can be viewed through a glass window 
7) in the outer wall. A system of pulleys and a 
counterweight balance the weight of the cryostat 
68 Ib) while leaving it free to rotate with the 
yoniometer table. The supporting cables can be 
seen at the top of the ervostat. 


* 
< 





Figure 2. The cryostat mounted on the \N-ray goniometer. 
ug 


1, X-ray source slit; 2, beryllium window; 3, detector slit; 4, X-ray recorder; 5, 
protractor; 6, specimen discharge tube; 7, window for viewing specimen surface 


4. Operation 


The system used to evacuate the cryostat and 
deposit gaseous samples is shown schematically in 
figure 3. In preparing for a run at liquid-helium 
temperature, the specimen chamber is pumped down 
to 10° mm of Hg and liquid nitrogen is added to 
the cold trap. Within 15 min the pressure drops 
below 3 107° mm of Hg and the cryostat is ready 
to be filled. , 

Before refrigerants are introduced, the gas to be 
deposited is usually pumped through the specimen 
chamber to simulate actual deposition of a sample 
on the helium-cooled surface. This aids in flushing 
the specimen manifold to remove contaminants, and 
permits the flow rate to be adjusted. After this 
adjustment the stopcock (fig. 3) between the capillary 
and the sample-gas input tube is closed so that the 
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Ficgure 3. Gas train and pumping system used with 
the cryostat. 


train up to that point is evacuated and gas does not 
collect in the tubing. The cryostat is then isolated 
from the pumping manifold by closing the valve 
and liquid nitrogen is introduced. Because of the 
cost of liquid helium, care must be taken to cool 
the apparatus to the temperature of liquid nitrogen 
before the helium is transferred. We prefer to fill 
the inner vessel to overflowing with liquid nitrogen 
and then syphon the liquid into the outer vessel. 
This method brings liquid nitrogen into contact with 
ell the surfaces that must later be cooled by helium. 
After thorough precooling of the apparatus, approxi- 
mately half a liter of helium is required to cool the 
inner Dewar to liquid helium temperature. 

The procedure for producing free radicals and de- 
positing them on a helium-cooled surface has been 
described by Broida and Pellam [10]. The system 
shown in figure 3 is suitable for most gases. but dif- 
ferent equipment is sometimes used to deposit a 
liquid (through the vapor phase) or the gaseous 
products of a chemical reaction. 

When the gas being deposited is to be dissociated 
to produce free radicals, it is passed’ through a silica 
tube in the rf field of a waveguide or diathermy cup. 
Power is provided by a 125-w microwave-diathermy 
unit operating in the 2,400- to 2,500-me band. With 
the power on, the gas flow is increased relative to 
the rate of condensation on the helium-cooled surface. 
When the flow is sufficient to give a pressure of 0.01 
to 0.1 mm of Hg in the silica tube, a Tesla coil can 
be used to initiate the discharge. The gas (or free 
radicals) is allowed to deposit until the surface of the 
copper block loses its metallic appearance. Fifteen 
or twenty minutes is required to deposit the 0.001- 
to 0.005-in. layer ordinarily used. 


5. Performance 


Once the sample has been deposited, obtaining 
X-ray patterns is routine. Although 25 percent of 
the X-ray beam is lost in passing through the 
beryllium windows, this loss is partly compensated 
when the path inside the cryostat is evacuated. For 
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copper radiation, the net loss is 15 percent. Inten- 
sities are adequate when the samples are crystalline. 
With the widest slits (3° beam-divergence slit, 
medium-resolution detector Soller slit, and 0.2° de- 
tector slit) counting rates considerably in excess of 
10,000 eps are often attained. However, we do have 
evidence that crystallites deposited on the cold sur- 
face are often highly oriented so that the intensities 
of certain lines are greatly enhanced at the expense 
of others. 

Helium is transferred from a 25-liter storage flask 
by means of a vacuum-insulated transfer syphon 
with a k-in. inner tube that extends 12 in. below the 
vacuum wall and reaches the bottom of the cryostat. 
With a pressure of 0.75 psi, less than 5 min is re- 
quired to fill the helium Dewar to the 2-liter level. 
The Dewar will hold as much as 2's liters of liquid 
helium, enough to last for 3 hr under ideal condi- 
tions. Depositing samples always increases the 
boil-off rate and reduces the time at liquid-helium 
temperature to about 2 hr. This is sufficient for the 
measurements required. Much of the data is taken 
during normal warmup of the Dewar and additional 
helium capacity would only delay this important 
phase of the experiments. 

If care is taken in alining the cryostat to allow for 
the thickness of the film deposited, there is no loss 
in the accuracy of interplanar spacings. Expansion 
curves are easily obtained by measuring the lattice 
constant of the specimen during wermup. 

Resolving power of the diffractometer is not ap- 
preciably affected by the ervostat. The diffraction 
peak reproduced in figure 4 shows partial resolution 
of the Cu Ka;-Ka, doublet at 43.46° 26 where the 
separation is about 0.1° 26. Complete resolution of 
lines only 1.2° 2@ apart is shown in figure 5. 

The apparatus has been useful not only in deter- 




















Figure 5. 


«1 portion of the diffraction pattern of annealed solu 
C*O), showing one peak due to the gold plating on the COppe 


specimen hlock. 


The pattern indicates the amount « nstrumental line bre 


pected 


midening to be» 


mining the nature of the deposits originally frozey 
out at liquid-helium temperature, but in following 
the changes during warmup as well. The sudden 
crystallization of amorphous deposits with the evo- 
lution of considerable heat has been observed in 
several Phase changes inferred from other 
types of measurements have been confirmed, and 
additional phases have been discovered. For ex- 
ample, X-ray patterns of the three phases of oxygen 
are shown in figure 6 along with a pattern of ozon 
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prepared in a separate experiment by dissociating 
oxygen in a rf field. Phe dissociation products were 
deposited at 4.2° K and allowed to warm up until 
the residual oxygen sublimed at about 50° Kk. A 
temperature of 56 _K was maintained by pumping 
on nitrogen in the inner Dewar. These results will 
be reported in more detail in later papers dealing 
with specific materials. ; . 
The performance of the X-ray cryostat during 6 
months of use has been entirely satisfactory. How- 
ever, the natural warmup characteristics sometimes 
do not give enough time In a given temperature range 
+o make the measurements required. Swenson and 
Stahl [8] have described a Dewar for maintaining 
temperatures between 4° and 80° K by circulating 
nelium gas cooled by liquid helium. Plans are being 
-onsidered for adapting their technique or providing 
vy auxiliary Dewar from which very small quantities 
of helium ean be transferred at frequent intervals to 
arrest the warmup. Other modifications may be- 
ome necessary if emphasis changes to the investi- 
ration of single crystals or orientation textures. 


1] 
[2] 


a 


[5] 
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The Second Spectrum of Ruthenium (Ru 11) 
Allen G. Shenstone! and William F. Meggers 


Wavelengths and estimated intensities of 
The wavelengths range from 1054.684 to 6371.29 A. 


ionized ruthenium atoms are presented. 


The Zeeman effect has been investigated for 488 lines ranging from 2323 to 4010 A, 


of these basic data of Rui has resulted in 
between 68-even, and 140-odd energy levels. 


into designated spectral terms and assigned to electron configurations. 


arise from 4d7, 4d®°5s, and 4d55s?. All the 
found, and 4d? atF represents the ground state 
series a limit, 135200 em-!, has been calculated 
(1. P.) of 16.76 eleetron volts for Ru* ions. 


As far back as 1930 [1]? an analysis was under- 
taken of the first spark spectrum of ruthenium, 
element number 44. Although the fundamental 
multiplets were easily found, the complexity of the 
spectrum made further work unrewarding until 
more complete and accurate measurements could 
be made. In 1955 |2] Kessler and Meggers published 
new measurements, and in 1952 Meggers measured 
Zeeman effects on plates taken at the Massachusetts 
Institute of Technology with their large-field magnet 
and high-dispersion spectrographs. In addition, 
the spectrum was photographed on the Princeton 
vacuum instrument to provide wavelengths from 
500 A to 2200 A. This instrument uses a 30,000 
line per inch grating of 2m radius of curvature in 
a normal incidence mounting. All of the measure- 
ments of wavelength were made from spectrograms 
of the spark, a source which is not altogether satis- 
factory for first-spark spectra. For instance, 
Dobbie’s [3] excellent analysis of Fer was made 
entirely from the polar lines in an are. Therefore, 
some further observations were made with such a 
source, and a considerable number of lines were 
found in the longer wavelength region which were 
not observed in the spark. All of the are and spark 
photographs were made with buttons of pure com- 
pressed ruthenium powder. The metallic impurities 


were very faint and easily eliminated. In_ the 
vacuum region the gas and carbon lines, which 
are always present, were used as standards. In 


that region the differentiation of stages of ionization 
was made from an examination of the polarity of 
the lines on a number of plates. In some cases, it 
was even possible to distinguish lines due to the 
lowest levels from those due to levels only 20000 
em~' higher. As usual, lines due to transitions from 
high even levels to intermediate odd levels are very 
wide and displaced from the position they occupy in 
an are source. 

Although the Zeeman effeet on ruthenium lines 
was photographed and machine-measured at the 
Massachusetts Institute of Technology in 1940 only 
the results for Rui were exploited and published 


—_—_—_—_— 


Present address 
ton, N. J 


° 


? Figures in brackets indicate the literature references on page 378. 


Palmer Physical Laboratory, Princeton University, Prince- 
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spectral lines characteristic of singly 


Analysis 
the classification of 1,633 lines as transitions 
All but three of these levels have been grouped 
The low-even terms 
terms from the 4d? configuration have been 
of Rut. From two members of the 4d¢ ns*D 
which indicates a principal ionization potential 


| [4]. In 1949 new Zeeman spectrograms of ruthen- 
| ium were made at the Massachusetts Institute of 
Technology but the magnetic field intensity did not 
exceed 80,000 oersteds. Because some of the earlier 
spectrograms were made with higher field intensities 
which greatly enhanced Ruu, the best of these 
(Z-74) were loaned to the National Bureau of 
Standards where they were hand-measured and com- 
puted. A field intensity of 88,200 oersteds was 
derived from the magnetic splitting of impurity lines 
(Mg it 2795.53 and 2802.60 A, Cul 3247.54 and 
3273.96 A, Agi 3280.68 and 3382.89 A, and Ca u 
3933.66 and 3968.47 A). 

The Rutt lines were easily distinguished from 
Rut by comparing the total intensities of the Zee- 
man patterns with the intensities in the central 
comparison spectrum produced without the magnetic 
field; the Rum lines were greatly enhanced by the 
magnetic field. In addition, the fully recorded 
and resolved patterns of Rui lines were recognized 
by their types 4, 5, and 6 characteristic of even 
multiplicities according to Back and Landé [5). 
Type 4 exhibits two symmetrically least-displaced 
p-components of equal intensity, and the two 
strongest n-components are also least displaced, 
when the level with the larger J has the smaller g. 
Type 5 occurs when the larger J has the larger g, 
in which case the two strongest n-components 
exhibit the maximum displacement. Type 6. is 
produced by levels with equal J but unequal g, in 
which case the strongest p-components show the 
largest displacement, and the strongest n-com- 
ponents an average displacement. The observed 
g-factors given to 3 decimals in tables 1, 2, and 5 

/ are the average of all calculations from resolved 
patterns, indicated in table 6 by Ag. Only the 
displacements of maximum intensities of unresolved 
patterns on the MIT spectrograms were measured, 
and the types noted. Then after such a line was 
classified it was often possible to calculate a g-factor 
for one level if the other one is known. In many 
type 6 patterns the over-all width (w, of the un- 
resolved p-components was measured so that g- 
values could be calculated for both energy levels 
when their J-values were fixed by the classification; 
in these cases the Ag (in table 6) is inclosed in 
parentheses and the g-factors are given to two deci- 


373 








mals. For a few levels, the g-values were calculated 
to two decimals from unresolved patterns on plates 
taken 30 years ago at Princeton, according to a 
procedure described by Shenstone and Blair [61. 

As remarked above, the lowest multiplets were 
tound many vears ago; but even after the new 
wavelength measurements were made and Zeeman 
effects observed, it still proved very difficult to break 
into the higher groups of terms. This was surprising 
since in a spectrum of this kind, one expects many 
intersystem and interfamily lines. There are many 
such combinations but few of them are with the sex- 
tet system. In point of fact, the strongest inter- 
family combinations with that system are in the long 
wavelength region where the lines can be observed 
only at the poles of the are. This was unfortunately 
discovered after the analysis was practically complete. 

The break into the higher quartet terms was made 
through a search in the Schumann region for the 
differences occurring in the low ‘F term. From such 
differences the usual square array was built up, though 
the work was not made any easier by the fact that 
great reliance was placed on a comparison with Fe 11. 


TaBLeE 1. FE 


Configuration Desig- J Level qd obs q Landé 
nation 
We a ‘Fk }! 0.0 1. 333 
hi a ‘éF 3! 1523. 1 1.238 
td a ‘Fk 4 2493. 9 1 ovr 
bi a ‘Fk | 3104. 2 0. 419 
te? a ‘Pp 2 8256. 7 1. 53 1. 600 
i? a ‘P l S477. 7 1. OS 1. 733 
HOD)S5 a *D }! 151. 4 1. 53 1 556 
if? a 4p 0 0373. 9 ” 60 ? 667 
{PODS a 3 10150. 4 1 576 1. 587 
tof 5 a 3 LOS5S1. 7 1. 641 1 657 
ii? a ty } 1LOS60. 9 1 ttt 
HOD A a *D {! 11303. 8 1. 843 1. S67 
tf Ds a *D ()! 11604. 0 a. 263 3. 333 
We a? 3! 12293. 4 0. S89 
h/ a?P | 129056. 6 1. 333 
hi? a?D 2', 14581. 2 1. 200 
td a?Hu 5! 14663. 4 1. O91 
td a?P )! 147040, 5 0. 667 
bed? a?H {! 16125. 0 0. 909 
hel a?D 1's 17017. 6 0. SOO 
$9(5])) 5s aé‘*D 34 19378. 7 1. 402 1. 429 
1d°®(°D) 5s a‘D 24 20515. 0 1. 337 1. 371 
{®(°1)) 5s a‘*Dp 1's 21246. 4 1. ISS 1. 200 
te? a?F 2! 21557. 8 0. 863 0. 857 
PPED)S a‘D 0! 21645.6 O. O15 0. OOO 
Vk a ?F 3!. 22280. 0 1. 138 1. 143 
$®@H) 5s aéH 6 25952. 2 1. 235 1. 231 
1d®GH)5s aéH 5! 26100. 4 1. 16 1. 133 
1d7°(H)5s a‘éH }! 26118. 4 1. OS 0. 970 
1PCOH)Ses a ‘H 3h. ~H4168. 0 0. 759 0. 667 
42° @P)5s8 bh 4p 2', 26911. 4 1.578 1. 600 
id® AF )5s b 4B t's 27544.6 | 1. 200 1. 333 
4POCUR)5s ° h 4k 3'4 27948. 8 1. OOF 1. 238 
41P UR) 5s bh 4k 2'% 28138.8 0. O41 1. O29 
1PUk)5s h 4k 1's 28495. 3 | 0. 420 0. 419 





Only after a very arduous numerical analysis was j 
possible to assign term designations to the even levels 
and a reasonable identification of odd levels was jm. 
possible until the analysis was nearly complete, { 
warning must be given against too great reliance on 
the names chosen. The levels are so mixed in charae. 
ter that the most that can be hoped is that designa- 
tions have been selected which correspond to the 
particular Russell-Saunders type that is the chief eon. 
tributor to each level. The terms of 4d? and 4d°OD)3. 
are the most reliably named because they are reagop. 
ably pure. The remaining even levels are in genera] 
less pure the higher they are. 

The names of the odd levels are considerably Jess 
reliable than those of the even levels. The levels Were 
found chiefly by the numerical analysis mentioned 
above, assisted by rough predictions of position and 
intensity of combinations. Because the spectrum 
Ter, isoelectronic with Ru it, is only partially ana. 
lyzed [7], the chief reliance had to be on the homolo- 
yous spectrum Fe i [S| for predictions of position, 
Cnfortunately, the correspondence is much rougher 
than was expected. For instance, the terms a ‘H and 


en levels of ud 


Configuration Desig- J Level giobs q( Landé 
nation 

$7° (303) 5s a ) POOLS. 6 1. 232 1. 273 
1! EUP ds h 4P | 29001. 0 1. 646 1. 733 
$9 (3(4)5s a 4‘ }', 300909. 6 1. 167 1. 172 
$®3Q)5s a AG D1, 30293. 3 0. 651 0. 571 
1° (30,)5s a ‘sy 3' 30430. 5 1. O21 0. O84 
$PUP)5s bh 4p ()! 30489. 3 2. 541 ”. 667 
1 OH) Ss h 2H 5I, 32623. 0 1. O65 1. OO] 
H/T) Ss h 11 }'., 32686. 7 0. 955 0). 900 
HAD 5s h 4) |! 32888, 4 1. 200 1. 200 
td® AD )5s h 4D) 0! 329060. 8 0. 138 0. 00 
td yds bh *D 2 33018. 8 1. 338 1. 371 
{PAD 5s h 4p 3! 33332. 7 1. 307 1. 429 
1PEAUP)Ss h 2p l 33734. 8 1. ISD 1. 333 
HPORS bh 2k 3 34038. 3 1. O83 1. 143 
if? h 21) | 34703. | O. O65 0. S00 
hel 5h *)) 245 34820. 3 1. 102 1. 206 
1A) As h 2h 2! 352908. 4 0. 945 0). 857 
1/5 Ds a “SS 2 B35857 t) » OOO 
HWeCT)S a? 6! 35939. 6 1. 06 1. O77 
$d®(3C,)5s » 25 }! 36016. 0 1. O87 1. 111 
{7 CAP)5s h 2p 0 36095. 0 O. 875 0. 667 
td ('T)5s a 2] 5! 36229. 7 0.93 0. 923 
1d) 5s h Gy 3! 36515. 8 0. 885 0. 889 
$7® (1) 58 ec 2Q) } 37433. 0 1. O78 1. 111 
$d°®('G) 5s eG 3! 37981.5 0. 942 0. S89 
1d®(])\5s er ?pD 2 B89O81.8 1. 18 1. 200 
1 1D))5s 1) l 39711.9 > OL 797 0. 800 
11 ('S)\5s8 aes (! 10226. 2 1. 800 » OOO 
td®('D)5es d?)) }' 12750. | 0. SO3 0. S00 
$d®('T))5s d?p Dl, 13609, 9 | 1. 216 1. 200 
1/°(0OD) 6s e %]) {', R4510. 9 1. 556 
W/HOD)6s e 4p 315 SO440. 4 1. 429 
WPCD)6 '}) 214 87523. 4 1. 371 
19D), (Rui) Limit 135200 





Col 





SIS Was it 
ven levels 
S Was im. 
plete. A 
liance oy 
Nn charae- 
designa. 
d to the 
‘hief con. 
41° (5D) )5¢ 
© Teason. 
n general 


ably less 
Vels were 
entioned 
tion and 
Pectrum 
lly ana. 
homolo- 
position. 
rougher 


a*H and 


Landé 


200 \ 


SOO 
000 
S00 
200 


556 
129 
371 








Configuration 


jd®(5D) 5p 
1°(8D)S5p 
16D) 5p 
16D) 5p 
td GD)5p 


iP)5p 
Op 5p 


1d63P »p 
1d63P »p 
td ‘F) p 
ed »p 
{dF »p 
163k 5p 
1GH 5p 


4/( P) Sp 


479794 


Desig- 
nation 


y ip 


ip 


ip 


‘H 
‘H 
GG 
iH 


i(} 
Cy 


47) 
4]) 
‘I) 


‘|) 
i(} 
i(} 


'() 
‘]) 


J 


RS WO te Or oer Nw 


We SO OO 


m Co bo 


ee | 


we ah 


evel 


LG471.0 
JA7 11. 5 
47285. 1 
17708.6 
17983. 8 
5O7 58. 3 
508 Ld 3 
5OS6 2. 5 
51179. 7 
51317.0 
51380. 1 
51548. 8 
52820. 3 
52964. 4 
53317. 1 
53685. 2 
5 L065 \ 
94225. 5 
54663. 4 
4794. 3 
54981. 5 
5S 2 24 3 
55695, 2 
56664. 9 
57 263.8 
62446. 1 
B25 26. 6 
AH2906. 8 
63307. 7 
A3518. 5 
83588. 2 
63797.6 
A3851 3 
63941. 2 
64025. 3 
aa L385 ; 
64 163 S 
A4748.8 
64921.6 
65212. 9 
65244. 7 
65259.8 
65326 7 
a5 47 2. I 
B5508. 1 
B5546. 5 
65776. 1 
66012. 4 
H6198. 6 
G6249. 3 
BB338. O 
HAA, 3. ; 3 
466651. 0 
A7099. 3 


47103. 0 


6 
a) 
6 


6 


7164. 3 
7401. 4 
7601. 7? 
7646. 1 


67900. 1 


TABLE 2. 


g(obs) 
1. 530 
1. 569 
1. 634 
1. 840 
3. 280 
1. 44 
1. 42 
1. 457 
1. 312 
1. 058 
0. 62 
1. 589 
1. 835 
1. 346 
1. 391 
2. 365 
1. 348 
1. 274 
1. 145 
1. O77 
0. O11 
0. 468 
1.570 
1. 708 
2. 648 
1. 798 
0. O61 
1. 024 
1. OOS 
1. 24 
1. 14 
1, 205 
lL. 136 
0. 978 
0. 646 
959 
1. 022 
1. 045 
0. 428 
1. 032 
1. 324 
1. 205 
1. 152 
1. 046 
1. 85? 
1. 14 
0, 991 
1. 199 
1. 09 
1. O4 
1. 343 
0. 789 
1. 279 
O. 855 
1. 187 
1. 158 
1.13 
lL. 166 
1. O2 
1. 234 


g(Landé) 


1. 556 
1. 587 
1. 657 
1. 867 
3. 333 
1. 455 
1. 434 
1. 397 
1.314 
1. 067 
0. 667 
1.714 
1. S86 
1. 333 
1. 429 
2. 400 
1. 371 
1. 238 
1. 200 
1. 029 
0. 000 
0. 400 
1. 600 
1. 733 
2. 667 
2. 000 
0. 667 
0. 970 
1. 133 
1. 273 
1. 23) 
1. 600 
1. 172 
0. O84 
0. 571 
. 727 
1. 029 
1. 238 
0. 400 
0. SSO 
1. 200 
1. 200 
1. 333 
0. 965 
1. 733 
1. O77 
lL. 111 
1. 143 
1. LOS 
0. 923 
1. 429 
0. 000 
l. 371 
0. 857 
1. 200 
1.172 
1. O91 
1. 273 
0. O84 
1. 200 
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Odd levels of Ru u 


Configuration 


1d®(G)5p 
1d°(3F) 5p 
11°F) 5p 
1d®(H) Sp 
1d°(3F)5p 


GE)S5p 

(3G) 5p 
LEUP)5p 
L8P)5p 
16( 


®3G)5p 


td 
til* 
1y 
ty 
ty 
®8G)5p 
if FE) Sp 
‘G)5Sp 
GG) 5p 
GP) 5p 


id 
4d 
td 
td 
td 


AGG) 5p 
d®8GQ) 5p 
d® 3G) 5p 
d® 8G) 5p 
dG8GQ)5p 


1°CbGQG 5p 


1d®GFk 5p 
1°) Sp 
1d°GP)4p 


'D)5p 
1G) Sp 
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J-value in the region below 80000 em™'. 


given an incorrect J-value. 


levels, and the J 


Numbers of levels of each J-value of the structure 
,d® 5p below 80000 em~ 


TABLE 4. 


No. of levels ] 7! 6! 5! }! >} 21 1! Ol 
Predicted 2} 6 13) 20: 25) 26) 23) 14 
Identified 2) 6 13) 20) 25) 24) 20) 10 
Observed (not identi- 

fied) | 2 


> ‘F have quite regular intervals in Fe 1m but in Run 
thev are so irregular that parts of the terms at first 
eluded discovery. Figure 1 shows these unusual in- 
tervals in the terms a*H, 6*F, 6*P and a‘*G. Of 
considerable assistance, however, was the fact that 
the positions of the sextet and quartet terms of d® s in 
Fe u when multiplied by % give a rather close corre- 
spondence with the same terms in Rui. The differ- 
ences between the higher and lower multiplicities 
based on the same ion term is greater by an average 
factor of about 1.30 in Ru 1. 

The naming of the levels was based on relative in- 
tensities of combination, on the Zeeman effect, and 
on the comparison with Fe 1. The various criteria 
were frequently at variance and compromises were 
therefore necessary. It will be noticed that hardly 
any of the g-values are in close agreement with the 
Landé values, and that the values tend to have less 
variation from the mean than do the Landé values. 
The worst disagreements are for the levels 2 *D3,, and 
y*D§,,, and it is quite probable that these levels are 
either misnamed or unaccountably perturbed. 

In Fe u there are known a number of even levels 
based on the high terms *F, *P, 'F, 'G from d® in 
Femur. A thorough search has been made for the 
equivalent levels in Rui without success, although 
there remain unidentified a number of lines which 
could be combinations of such levels with known 
terms. “This failure is all the stranger when it is 
noticed that some odd terms of this same group have 
apparently been discovered. It is possible that the 
observations are not sufficiently exhaustive. Some 
plates taken at Princeton definitely show a consider- 
ably increased number of lines. 

Tables 1 and 2, respectively, list all of the known 
even and odd levels in order of magnitude, and Table 
$3 relates them to the atomic structure. The levels 
are all given Russell-Saunders names except three 
odd levels which are numbered. Table 4 lists the 
numbers of 4d° 5p levels to be expected with each 
It is satis- 
factory that phe numbers discovered are so nearly 
those predicted, but two discrepancies must be 
noted. The most striking is the presence of an extra 
level of J=3' 
all the J 
them appears in any way spurious or to have been 
The two unidentified 
1's do not present such a serious prob- 


levels of m 


lem since there is room for them numerically even 
» if they do not fit the names and positions of any of 


The J=2% column lacks two 
0 four levels, but the latter fact 


the missing levels. 


A quite thorough examination of 
3's levels has been made and no one of 


is not surprising because such levels are always 
difficult to find. 

It was satisfactory to find the 4d° 5s? °S level and 
three of its associated P° terms, y,z *P° and w*P®. 
The odd terms are confirmed by strong multiplets 
in combination with a °D, and the °S itself combines 
with 2°P°. The *P®° term of the odd structure was 
found from its combinations with a*D in the Schu- 
mann region. An unsuccessful search was made 
for 4d° 5s(7S)5p *P°. 

Most of the data in tables 1 and 2 were supplied 
in October 1957 for inclusion in Volume III of 
Atomic Energy Levels [9|, where all the levels are 
grouped in spectral terms and the terms themselves 
are arranged in order of magnitude of the lowest 
level of each. This style of presentation has been 
used consistently in Atomic Energy Levels; it has the 
advantage of displaying directly the intervals be- 
tween levels of terms and thus showing the actual 
splitting as well as whether any complex term is 
normal, inverted, or partially inverted. Since some 
slight revisions of level values, intervals, and g-factors 
for Rui have been made during the past year our 
final results are shown in table 5 in the same style as 
in Atomie Energy Levels |9|. A large majority of 
the Ru terms are inverted as expected, because 
they arise from configurations containing six or 
seven type-d electrons, that is, the d-shell is more 
than half filled. However, some terms produced by 
configurations containing six d-electrons are normal 
or only partially inverted, and even d’, expected to 
give inverted terms, yields two terms, a?F and 
b?D, that are erect instead of inverted. 

For pure LS coupling the intervals in complex 
terms are proportional to the larger J-values defining 
an interval. Thus, for 4d’ a ‘*F the intervals should 
be proportional to 4.5:3.5:2.5 whereas they are 
actually in ratio 4.5:2.9:1.8. Similarly, the in- 
tervals between levels of 4d° (°D)5s a *D should be 
proportional to 4.50:3.50:2.50:1.50, but their ratios 
are observed to be 4.50:3.16:2.04:1.36. The ir- 
regular distribution of levels in four 4d° 5s terms is 
illustrated in figure 1. The remaining Ru um terms 
with three or more levels all flagrantly violate the 


interval rule. These violations, as well as the 
departures of g-factors from the Landé values, 


indicate that LS coupling is far from perfect in the 
ruthenium ion but nevertheless strong enough to 
permit plausible designations of the spectral terms. 

The line list (table 6) is given complete although 
a considerable number of lines remain unidentified. 
The majority cf these are in the Schumann region 
and many are of large intensity. They probably 
have their origin in terms of the 4d° 5s 5p structure, 
a few of which have been found. At 1937 A there 
is a pair of very strong lines which have been worked 
over very completely, but unsuccessfully. In the 
longer wavelength region there is a scattering of 
unidentified lines of some intensity and even two 
with resolved Zeeman patterns. They may be due 
to terms built on high levels of 3d° or to terms of 
3d° 48°. A few lines have been remeasured on 
plates taken in Princeton on high dispersion and have 
been separated into pairs or threes. 
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Fiaure |. A plot of four Ru u terms from 4d® 5s. 


Table 6 contains the measured wavelengths and 
estimated relative intensities of 2,227 spectral lines 
ascribed to Ru um in the wavelength range 1054.684 
to 6371.29 A. The wavelengths shorter than 2000 
A are valid for vacuum, but the longer ones are 
proper for standard air. For purposes of analysis 
and interpretation all wavelengths have been con- 
verted to vacuum wavenumbers, ¢, in K (=em~') 
units in column 3. The 4th column contains the 
designations of two energy levels (tables 1, 2, 5) 
which differ in value by the amount shown in column 
3. 
For wavelengths greater than 2322 A six additional 
columns are required in table 6 to report the Zeeman 
data for 488 Ruri lines. Instead of giving all the 
observed components of each Zeeman pattern, only 
the strongest p- and n-components are reported in 
columns 9 and 10, respectively. This information 
together with the Zeeman-type numbers [5] in col- 
umn 5, and the g-factors in columns 7 and 8 corre- 
sponding to the J-values in column 4, permits the 
reconstruction of the complete Zeeman pattern. 
Actually the absolute values of J and g were first 
derived from fully observed and resolved Zeeman 





patterns indicated by Ag in column 6.° Later 


: ; » When 
many more lines were connected by the usual methods 


of interval analysis, most of the unresolve 
patterns could be interpreted as 
analysis and extending the system of g-factor 
Only a few Zeeman patterns were too unsymmetried 
(us) to be useful. ' 


ed Zeeman 
confirming th 


- 


There are only two resolved aal 
four unresolved Zeeman patterns that appear { 
disagree with the classification indicated in column 
4. These are marked with asterisks accompanying 
the type number; they may be ascribed to errops 6 
observation in a rich mixture of Zeeman patterns 
belonging to Ru, Ru 1, and to impurities, ; 

The total number of Ru 11 lines classified in tables 
is 1,633, and, excepting the spectral range 1054 ty 
1400 A, there remain relatively few lines of any 
importance to explain. , 

Series Limit. Because only two members of the 
4° ns °D series are known, an accurate limit eanno} 
be calculated directly because such series usualy 
obey the Ritz formula 





and the constant @ is by no means negligible. Ap 
estimate of the value of @ is, therefore, necessary and 
this can be obtained only by a general comparison of 
spectra of both the first and second long periods, 
Catalin and Rico {10} have made such a comparison 
and they deduce a value of a=4.09x107°, With 
this value and the correct value of the separation 
between a°D and e°D, the limit is at 135200 em™ 
to the nearest 100 em™', giving an I. P. of 16.76 y. 

In a private communication the late Professor 
Catalan furnished a term table of Rhu, a spectrum 
isoelectronic with Rui and Ter. The agreement ig 
excellent in the parts based on the °D ion, but some 
of the 4d‘ terms and many of the high odd and even 
terms show a considerable deviation from the posi- 
tions which would be predicted from our analysis 
of Rutt. 
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TABLE 6. The second spectrum of ruthenium (Ru u)—Continued 
l 2 3 } 5 6 7 8 iT) 10 
Intensity » (air) og Term combina- Type Ag g J Strong p| Strong I 
tion 
A AK 
8 2472. 99 10424. 7 a *Po 2 *F ix, 
10 2473. 57 10415. 1 
5h 2473. 83 10410. 9 
0 2474. 76 10395. 7 b4p w ADS, 5 L. 200 1. 30 r 1. 44 
1 2476. 15 10373. | b 2G w 2H 
10 2477. 17 10356. 4 b 4F r 4F3 
50 2477. 26 10355. 0 64D e4'Dp 
S 2477. 98 10343. 3 
20 2478. 36 10337. 2 
500 2478. 93 10327. 9 a *Ds 2 §F3u 6 0. 354 1. 635 1. 301 O. S51 1. 468 
P 2480 20 — b 4*Ba—z 2*Pix |) 
5 2480. 2¢ 10307. 0 b 4D r ‘Pi, 
30 2480. 81 10297. | a‘*H z 218. 5 IS] 1. 23 1. O5 09 2. 16? 
200 2481. 11 10292. 4 atG y 4H ! 1. 237 1. 20 . 20u 1. 02 
15 2481. 85 10280. 4 b *F t *Ding 
60 2482. 78 10265. 2 b4*D ‘Dp 6 0. O44 1. 345 1. 301 110 1. 323 
l 2483. 20 10258. 5 
3 2483. 54 10253. 0 
50 2483. 96 10246, 2 a‘H 2 ‘Lin, 6 0. 10 1. 23 1. O83 67 1. 13 
5 2484. 49 10237. 5 a ?} 2 {Hi 
10 2484. 67 10234. 6 b *Dy r *P5 6 1. 62 0. 14 1. 76 SI 0. 95 
8 2487. 20 10193. 7 b *D. y ®Phrg §* 0. 268 1. 306 1. 574 . 135 2. 244? 
tH] 2487. 51 1O1SS8. 7 b *P : 
3 2487. 90 10182. 4 a‘H xr 4] ? 
30 2488. 57 10171. 6 b 4F 2 *Hii¢ { 1. 093 1. 02 23u 0. 76 
> 2488. &S 1O1G6 bo 
50 2489. 34 10159. 2 a ?] oe? Hy | 0. 90 0. 91 ( . 87 
6 2490. 54 10139. 8 a‘éH z 2] 
50 2491. OS 10131. 2 a *Gyac—y 2D 
10 2491. 10 10130. 8 a*Hyc—z 71k, 
5 2491. 34 10127. 0 b 4F- y *Ghu 5 0. 948 1. OPS ' 1. 24 
10 2491. 56 10125. 4 b4*Ds w iF 7.4 1. 344 1. 32 0 1. 26 
5 2493. 59 10090. 9 b *F» r 4D 
300 2493. 68 10089. 2 a*Hsy—z ‘Tix ay 0 1. 124 
10 2494. 16 LOOS 1. 6 h‘*Dp o4F 
10 2494. 48 10076. 4 a*D z ®F5 5 2. 50 1. 83 0. 67 1. 25 3. O08 
2h 2495. 35 10062. 4 
200 2495. 69 10057. 0 b 4F, J *Ging j 1. 200 117 u 1. O38 
2h 2495. 94 10053. 0 
25 2496. 86 10038. 2 6 0. 26 0. 98 
l 2497. 36 10030. 2 b‘*D r *Hjy? 
. | 
200 2498. 41 10013. 4 a s—z *Fix ? 
200 2498. 58 10010. 7 a*D z 8 Fu, | 0. 180 1. 637 1. 457 O09] 1. OO7 
Ss 2499. 38 39997. 8 b AFau—y *G5 
3 2499. 56 39994. 9 
50 2500. 14 39985. 6 a 2I¢ y 715 5 0. OS) 1. OS 1. OO .4lu 1. 52 
100 2501. 95 39956. 8 b 4F, z *Ha, | 0. O06 1. 20 1. 14 25u 0. 87+ 
" b4D w iF 
9EnO 2 aa 2 2 ; 
1Oh ~~ 1 oe 30) 3995 1 » hb ‘T) , iD: 
5H 2503. 22 39936. 5 
10 2504. 92 39909. 5 b Gy —y 215 
OEnk aA - a a *Ciy z 8 
10 2505. 64 39897. 9 b 4D 2s 
4 
15 2505. 89 39894. 0 atHys—z 2F hy 
25 + 2506. 25 39888. 2 a *Gsue—r *Ging 
25 2506. 42 39885. 5 b4‘D w Fou 
300 2507. 00 39876. 3 a *Du—z ®*F iy 1 0. 536 1. 869 1. 333 267 0. 592 
1 2507. 50 39868. 4 
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TABLE 6. The second spectrum of ruthenium (Ru 11)—Continued 
. 
l 2 3 | 5 6 ra 8 9 10 
Intensity A (air a Term combina- Ty pe Ag q q Strong p Strong . 
tion 
A A 
75 2628. 75 38029, 5 b 2G), r 23 7* d? 1. O8 
5 2629. 37 38020. 6 a4, z 2Hi. 6 0. 159 1. 167 1. OOS 0. 716 ? 
150 2630. 04 38010. 9 a 2] r *H2,. 7,4 1. OS 1. 09 0 1. O51 
20 2631. 09 37995. 7 a‘H / IFS | 259 0. 758 1. O17 0. 126 0. 110 
25 2632. 73 37972. 0 aitG 4 2G 5 0. 38] 0. 652 1. O35 0. 19] l. 604 
5 2633. 50 37961. 0 
20 2633. 82 37956. 4 h 4#P 2 2p ! 1. 372 2. 541 1. 169 OS6 0. 483 
10 2635. 21 37936. 3 a‘éG r 4p 5 0. 305 0. 647 0. 953 154 1. 718 
10 2635, 34 37934. 4 bh 204 H 
70 2635. 84 37927. 2 bh 4} ‘| | 165 l. 207 1. O42 OS2 0. 300 
100 2636. 54 37917. 2 aéH g 4]: 5 200 0. 748 0. 948 100 1. 648 
20 2636. 84 37912. 9 c2q G3 6 0. 06 1. 09 1. O83 2d 1. 065 
10 2638. 35 37891. 2 h 4]) ¥ 2p 6 0. 457 lL. 205 0. 748 697 0. 977 
| 2639. 34 37877. 0 
15 2639. 58 37873. 5 h 4} | 5 254 1. 193 959 127 1. 828 
50 2641. 63 37844. 2 h 2] ‘] 7 0 1. O88 
5 2642. 10 37837. 4 Db *P j ‘I | 25 né2 1. O47 0. 262 
20) 2642. 80 37827. 4 b 4] (ij 
30 2642. 88 37826. 2 CG 4 7G) 
20 2643 14 37822 a) a ‘H > 4, 
3 2643. 40 37818. 7 64D 4?D 
3 2643. 60 37815. 9 r 2] r *Gij 
10) 2644. 62 STSOL. 4 aad r 4h 6 1. 167 1. 19 1. 18 
75 2645. 97 37782. 0 h 4] 4 IFS 6 1. 20 1 15 23 1. 18 
100 2648. 78 37742. 0 a‘+éH ‘(iu . l. 16us 
3 2650. 26 37720. 9 a 2 H 
10 2652. 94 37682. 8 h*D 4?D ! 1. 200 1 15 l. O7us 
S 2653. 10 37680. 6 até H; 7 1. OP 1 1 O16 0 1. OO6 
25 2653. 95 37668. 5 h 4] > 2G 5 170 l. 202 1. O52 0. OS6 1. 797 
l/ 2654. 524 37660. 33 
10 2654. 876 37655. 34 b ?P, ‘Di..? 
5 2655. OOS 37653. 47 a*s 4 4G) } 
300 2656, 235 37636. 07 a‘éH 'H 6 0. OS 22 1. 14 0. 51 1. 18 
25 2657. 198 37622. 43 
10 2657. 326 37620. 62 
Ld 2658. 244 37607. 63 h 4p ‘1D 6 2.53 0. 37 1. OS 1. 45 
l 2658. 672 37601. 58 
l 2658, 922 37598. O04 
l 2660. 152 37580. 66 e?D y ?P 
200 2661. 169 57566. 30 aéH 1G 7 1. 23 |. 24 0 l. 1438 
300 2661. 610 37560. OF a°*[) z 6 7 1. 54 1. 56 0 lL. 502 
20 2661. 82 ST DDT. 11 aéH 2 4G 
15 2662. 183 37551. 99 bh4D /?D 
25 2662 SSO) 5742 6 h ‘p / ‘T) 5 0. S56 l 652 () 796 0. 428 2. OSO 
10 2663. 754 37529. S84 b 2G r 2Q3 | 174 0. S&3 709 O87 hy 
2 26065. 376 37507. 00 a‘G 4 4G 5 330 1. O21 69! 162 ? 
100 2667. 390 37478. 68 aéH ‘H lus 1. 16 1. 14 ' 1. 02 
10 2667. 79 37473. O7 aéH (| 6 298 0. 752 0. 980 0. 798 0. 866 
? 2668. 632 37461. 25 ait r ifs 
75 2669. 43 37450. 05 2G r2G 6 0. O31 943 913 107 . 928 
LS 2670. 527 7 4AS4. 66 6 117 1. 120 
25 2672. 212 37411. 06 b *P; 44D 
25 2672. 354 37409. 07 a‘ > 44 | 
10 2672. S80 7401. 71 a‘, z 2H 
30 » 2673. 004 37399, 98 a‘H 2 4G 
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TABLE 6. The second spectrum of ruthenium (Ru) Continued 
10 2 3 4 5 6 7 8 sy) 10 
Strong » Intensity d (air) o Term combina- Type Ag g g Strong p| Strong n 
: tion 
1 kK 
1. O8 10 9674. 219 37382. 99 a*Hy.—z ‘Gig 5 0. 380 (0. 753) 1. 133 . 190 ? 
? 10 2675. 197 37369. 32 | b *Fas—y *Piy 2 
1. 05] 25 2675. 543 37364. 49 ¢ *G3u—w 2Giy ; ies 
0. 110 100 2676. 183 7355.55 fa *Hgy—z tH, 5 0. 143 1. 24 1. 10 0. 79w 2. 03 
l. 604 la Sas—y ®Phx. 
{ 2676. 775 37347. 29 b °Gys—ar AHGy 6 . 104 (1.067) 0. 96 47 ? 
0. 483 3 2678. 054 37329. 46 a*Hyc—y *Psy a. Z 
l. 718 800 2678. 759 37319. 63 | a®*Dyc—z Dix 7 5 1. 532 (1. 530) 0 1. 531 
5 2679. 450 37310. O01 0 1. 55 
0. 300 100 2680. 585 37294. 21 b?Gauy—ar *Ghy 6 . 221 0. 893 1. 114 0. 794 1. 003 
2 2681. 293 37284. 36 a So wiFs a 
1. 648 % : 
l. 065 | 2681. 417 37282. 64 ‘ on enghiaie 
0. 977 3 2682. 632 37265. 76 
5 2682. 754 $7264.06 | b *Fyas—z 2Ghy 6 (0. 05) 1. O9 1. O4 .19 1. O74 
1. 828 3 2682. 894 37262. 12 b *Dus—lix ical 
15 2684. 613 37238. 26 b*Dy—y ? Day 5 0. 249 1. 406 1. 157 125 2. 028 
1. O38 
15 2685. 152 37230. 78 a *G5ig—z 2T8x 6 . 197 (1. 237) 1. 040 1. O84 ? 
2 2685. 581 37224. 84 b 2Fas—w 'Fix 4 a 
» 2685. 893 37220. 51 b ?Po—w *Dixg : vote 
3 2686. 418 37213. 24 ann 
15 2686. 889 37206. 72 a 4*Gy—y *Ghxy 1? (1. O21) 1. 03 u 1. 02+ 
25 2687. 071 37204. 20 b Fay Fey 5 . 147 1. 186 1. 039 0. 072 1. 700 
18 100 2687. 494 37198. 34 a*Hs—z *Hix, ; p : 
18 80 2688. 147 37189. 31 a*Hy—z *Hiy . ; ee 3 
. 16us l 2689. 431 37171. 55 : , ees 
10 2690. 208 37160. 81 ¢ *Gys—w ! Diy } (1. O78) 1. 11 u 0. 97 
, O7us 100 2692. 120 37134. 43 a *Dy—z *Dixg } 1. 576) 1. 71 w 1. 444 
006 2692. 716 $7126.21 |b ?Gys—w A Fy 
797 20 2693. 628 37113. 64 b *Pa—z 'Gixy 
I 2694. 199 37105. 77 ob *Fag—z 2Dix . 
l 2694. 598 37100. 28 b ?Pow—r 'Pixg 
} 2 2696. 495 37074. 18 b 4Fa.—z 2Ghuy 
18 20 2697. 07 37066. 3 a *Sou—y?Pix 1 0. 516 1. 284 1. 800 0. 258 1. 026 
10 2697. 12 37065. 6 a *Gys—y *Gix 
10 2698. 167 37051. 21 a Sac—y ®P5: 
20 2698. 268 37049. 82 : w 1. 16? 
15 
3h 2698. 767 37042. 97 w----- 
5 2699. 228 37036. 64 b?Dag—r *Pix } . 480 (1. 102) 1. 582 0. 240 ? 
50 2700. 163 37023. 82 b?Hs—y *Hix, 6 (0. 031) 1. O77 1. 108 . 172 1. O92 
143 25 2700. 999 37012. 36 b *Fig—y *PRY? ; acu m 
10 2702. 639 36989.90 = b?Hs.—y *Hin, ; f = 
02 
15 2704. 192 56968. 66 b?Hy—y 4H 3h 5 0. 146 0. 958 0. 812 . 073 1. 469 
80 2704. 585 36963. 29 a *lec—z 2K% 1 (1. 08) 1. 07 w 0. 99+ 
80 10 2704. 830 36959. 94 b?Hya—y tHe, 5 . 160 0. 98 1. 14 0. O80 1. 70? 
? I 2705. 568 36949. 86 b ?Fa.—zx 2F% J - 
lh 2706. 936 36931. 19 
2 50 2707. 310 56926. 09 6b? Ay —y Hi, 6 (0. 04) . 95 0. 99 . 19 0. 97 
66 2 2707. 555 36922. 74 - -* 
l 2709. 000 369038. 05 a Say—y ®Pix " 
28 100 2710. 228 36886. 33 b *Paw—y * Psy 6 0. 281 1. 576 1. 295 . 709 1. 436 
300 2712. 409 36856. 67 a *Do—z *Dix { . 201 1. 637 1. 838 . 101 1. 335 
20) 
10) 2713. 071 36847. 68 b *Gy,—ar *Hix ; vee ef 
1) 2713. 585 36840. 70 b 4*Fas—z Alix 6 . 247 1. 203 0. 956 1. 120 1. O8O 
l 2714. 717 36825. 34 “ . . . 
l 2715. 101 36820. 13 j : 
50 2716. 132 36806. 16 a *Gas—z 2Fiy, 6 19 0. 65 . 84 0. 46 0. 745 
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TABLE 6. The second spectrum of ruthenium (Ru i1)—Continued | 
2 3 } 5 65 7 S 9 10 
Intensity \ (air o Term combina- Type Ag i] gq Strong p Strong » In 
tion 
A A 
10 2716. 580 386800. 09 bh 4F y FS. 6 0. 03 1. O9 1. O6 0. 10 1. OR) 
l 2716. 78 36797. 38 a‘H z Hix, 
SO 2717. 447 36788. 35 a‘*H, 2 (HG 6 0. 05 1. OS 1. O38 22 1 O06 
10 2717. 855 36782. 82 b Fo y *F ing 5 0.514 0. 94) 0. 427 257 1. 712 
20 2719. 717 36757. 64 a *Gyug—z PFGiy 5 238 1. 218 982 118 2. 220 
I 36749. 28 bh 4F z 2D; 
l 36725. 72 a‘éG y AGG : 
20 36711. 63 bh 25 y Fa 5 165 0. 947 L. £49 O82 1. 524 
2 36702. 38 h 2G v4 F 5 
2 3669 L. 36 bh 2P, 2i 
100 36688. 21 b?H y 4H 5 127 1. 053 1. 180 O62 1. 878 
300 36680. 13 a*D z §D5 { 1. 448 1. S40 3. 288 72 l. 116 
20 36659. OS b4p r {Fe 
2 36636. 76 b ‘Ds: , He 
2 36627. 38 b*D r Fay 6 0. 36 1. 344 0. 98 90 
Du 2730. 413 36613. 65 b *P r 4—D 
5 2730. GSS 36609. 97 h 4] , +] 
3 2731. 389 36600, 58 bh4*p r 4F 
3 2732. 685 36583. 22 
00 2734. 345 36561. O1 a*D z §p 7.6 1. 576 1. SS4 0 1. 580 
10 2736. 456 36532. 81 bh 2] 4?D 
60 2736. 826 36527. 87 atG z 2] j 0. 11 1. 26 1. 15 0. 53 0. 62 
l 2737. 384 36520. 42 
10 2737. 606 36517. 46 ¢ 2G r?H 7 1. O78 1. O74 0 1. O56 
20 2737. 783 36515. LO b 4F 1 IF Sie 5 0. OF 1. 09 1. 02 0. 17 1. 26 
20 2739. 372 36493. 92 b?H, maar 6 0. 12 0. 95 1. O7 53 1. O1 
Sh 2741. 178 36469. SS 
25 2742. 401 36453. 61 a‘ z 4] % 6 195 1. 245 1. O50 1. O73 1. 148 
l 2742. 691 36449. 76 e ?D: v 2Ghy 
25 2743. 513 56438, 84 aéH 2 4H; 5 0. 265 0. 761 1. O26 0. 134 1. 953 
5 2743. 700 36436. 36 bh 4F z 415 ) 
60 2743. 93 36433. 25 a*D z §—; 7 1. 637 1. 643 0 1. 640 
I 2744. 29 36428. 5 b 2p v Fig 
100 2745. 158 36417. 00 b *P y *P 
50 2745. 827 36408, 15 a‘*H, z 4H 
20) 2746. 068 36404. 94 a*D 2 8): 7.6 1. 843 1 £0 0 1. 82? 
25 2746. 695 36396. 63 bh ‘F, 2 4G 5 219 1. 192 0. 973 0 110 1. 958 
10 2747. 563 36385. 13 h4*Dp r 4F j 297 1. 361 1. OO4 148 0). 322 
50 2747. 963 36379. 83 a*D > De. 7 . = 32973 0 3: 
5 2749. 100 36364. 79 d?p 2 
1d 36363. 07 7 0 0. 93 
Dh 36357. 65 a *as r ‘]): 
10 36325. 02 b 4F, y 4B 6 0. OS 1. O83 1. 11 0. 21 1. O7 
30 36320. 57 a*D z §—; | 0. O4 1. 59 1. 55 Louw 1. 39-4 
50 36316. 40 a‘éDp > 4p j 0. 175 1. 399 1.574 OSS 0. 961 
20 3. 508 36306. 58 b 4F, > 4G 6 0. OS 1. 26 1. 18 35 1. 22 
2 53. 83 36302. 32 
20 2755. 258 36283. 52 bh2?H r 4G3 | 0. 09 1. O83 b. ie 39u 0. 64 
1 2755. 519 36280. OS h4+p , 4H 
| 2755. 748 36277. 07 
1D 2756. 411 36268. 34 
lh 2756. 899 36261. 92 
5 2758. 347 36242. 89 b?D s | 1. OS7 0. 975 2. O62 544 131 
2 2758. 692 36238. 35 a‘éG , 4D 
10 » 2759. 159 36232. 22 bh2G r 2k | 0. 06 1. OS 1. 14 21 87 
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Strong 


lL. OS) 
l. O64 


2. 280 


1. 580 


. 62-4 


lL. O56 


|. 26 


07 
39-4 
961 


64 


n 














[ntensit\ 


g 2760. 115 
i) 2760. 745 
3 2761. 473 
10 2762. O71 
20) 2765. 134 
100 2765. 429 
15 2765. S89 
Rg) 2766. 563 
950 2768. 926 
| 2770. OSS 
0) 277 O60 
5 2771 O4 
3 2772. O11 
100 2772. 459 
10 2774. 197 
30 2775. 631 
50 2777. 401 
50 27744. 5A 
| 2778. 024 
150 2778. 388 
30) 2778. 975 
15 2770. 406 
Sh 2780. 817 
8 2782. 23 
10 2782. 361 
S 9783. 69 
} 2783. 799 
10) 2784. 516 
10 2785. 186 
1) 2785. 741 
25 2785. 87 
15 TRG. B83 
20) 2787. 25 
200 2787. 823 
25 2788. 729 
ITSO, 453 
5 2790. OS2 
15 2790. 240 
a) 2791. 193 
100 2792. 32 
5 2704. 002 
20h 2794. 278 
10 2795. 349 
a) 2745. 555 
5 2795. 993 
Ss 2796. 461 
2/1 2797. O77 
10 2797. 747 
25 29708. 779 
bat 27909. 587 
15 9799. UIT 
| 2800. 408 
2 2800. 586 
30 2802. 152 
3 2802. 717 


TABLE 6. 


o 


A 


36219. 67 
36211. 41 
36201. 86 
36194. 03 
36153. 93 


36150. O08 
36144. 06 
26135. 26 
26104. 42 
26080, 28 


36076. 62 
36071. 36 
36064. 24 
36058. 42 
36035. 83 


36017. 21 
35994. 26 
35992. 46 
US6. 19 
USL. 48 
5073. 87 
5968. 30 
5950. 05 
35931. 79 
35930. 10 


35912. 95 
35911. 54 
35902. 29 
3589038. 66 
B5S886. 51 


35884. S4 
35878. 24 
S5867. OS 
35859. 71 


S5838 7D 
35830. 67 
5828. 65 
s5S816. 41 


ay iD. 20 
35776 87 
s-. 763. 17 


35709. 03 


35704. 70 
35698, 56 
35696, 29 
35676. 34 
35669, 15 


The second spectrum of ruthenium (Ru it) 


Term combina- 


b? Hg. 
a *Gaug 
d 2D. 
b *Pix 
a ‘Dox 
b *Po. 
? Dax 
a*D 
4h ie 
h ‘[) 
¢ 2Ggx 
1 Hy. 
b Dor. 
a ‘DD, 
b 2G gg 
h if 
4 D 
a ) 
b 4B, 
I if, 
42D 
c2Qy 
a Mag 
b 2P, 
h if 
d2D, 
4h 1F 
bh 2Py 
d 2 Dig 
a*D 
64D 
h Py, 
ee? 
b 2F, 
h ‘FP, 
42D 
hb Fy, 
h4p 
a Ss 
h 2G 
a Myron 
c 2G 
b 2P, 
a *T5 
a ig 


1 


tion 


Gx 


2 Ding 
‘Dix. 


2G3 


ye 


25 


‘Psu 
G3. 


‘H., 


‘Dix, 


Po, 
“TT Gu, 


Continued 


5 6 7 8 
Type Ag g g 
| 0. 260 lL. O21 1. 281 
} 0. 160 1. 206 1. 046 
t . 322 1. 643 1. 321 
} » a 1. 339 1. 711 
7 1. 18 1. 20 
) 146 1. O85 0. 639 
6 200 0. 764 - . 964 
| 1. 452 1. 189 2. 641 
| (1. 067) 1. 10 
6 0. 11 1.10 0. 99 
} 210 1. 843 1. 622 
5 0. 03 1. 23 i. 26 
5 0. 614 0. 415 1. 029 
y 354 1. 107 . 463 
5 228 1. 053 . 825 
6 12 1. OS . 96 
7 1. 167) 1. 1& 
5 (1. 093) 1. 14 
t . 194 1. 227 1. 033 
6 . 292 0. 944 0. 652 
6 0. O7) 0. 80 0. 87 
| 0. OS) 1. 68 1. 60 
} 0. 333 (1. 397) 1. 063 
5 0. 943) 0. 99 
D . 182 1. 120 . 938 
| 134 (1. 200) 1. O64 
6 254 (0. 883) 1. 127 
) 545 0. 649 1. 194 
5 268 1. 197 0. 829 
0. 90) 91 
6 173 1. 168 . 995 
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Strong p 


130 


0. OSO 
. 160 


. 186 


. 708 


. 728 
w 
0. 38 


. 105 


. l4w 
. 308 
.177 
. 114 


. 03 


0. O97 
. 743 


0. 10 
. 2lw 
167 


0. O92 
. OOS 


. 8389 


. 274 


10 


Strong n 


. 200 


. 864 


. 463 


. 95-4 


. 05-4 


= 
wa 


. 393 
” 


. 557 


. 381 


. 86h 
. 082 














Taste 6. The second spectrum of ruthenium (Ru i1)—Continued q 
" 
l 2 3 | 5 6 7 s 9) 10 
Intensity X (air) a Term combina- Tyre Ag g q Strong p Strong » Int 
tion 
A AK 
j 
1 2803. 535 35658. 75 b 4Fog—y Pig 6 353 0. 943) 1. 296 883 ° y 
50 2804. 903 35641. 36 b?Hys—ar *Giug 5 198 0. 954 1. 152 . O99 2. 043 
100 2806. 77 35617. 65 a‘D ‘Pixs 6 2. 654 0 2. 654 1. 327 1. 327 
25 2807. 200 35612. 19 ¢ 2 Dow—w * Ding 1 lL. 16 1. 14 v lL. OG) 
. ~ of orn fa *Ci, z §—; 
9 ? wn } ‘ z 
5 2807. 350 35610. 29 le 2D r 2D; 
20 2807. 554 35607. 70 b 2Fyy—ar AF Sy 5 0. O8 ) 1. O09 1. O1 0. 20w 1. 29—9 
S 2810. 249 35573. 56 b4*Dyay—ar Gig 
250 2810. 649 35568. 50 a *se—z 2 KGa 7 0. 90) 0. 91 0 0. 94 
10 P81. 545 35557. 16 b 4D z *Ph., 1. 200) 1. 18 1. 19 
10 2811. 931 35552. 28 
l 2812. 6S9 $5542. 70 b?Hy—a *Diy 
100 2813. 311 35534. 84 b *P. z ‘Sin ! 0. 234 1. 564 1. 798 0. 118 1. 213 
60 2813. 694 35530. 00 b 4F z 1Giu 5 230 0. 418 0. 648 115 0. 994 
2 2814. 952 35514. 13 b 2] r 4G 
2 2815. O85 35504. SS d?Dp yp 2} 
l 2816. 281 35497. 37 b?H z 2Hi: 
1 2817. 255 35485. 10 b *Doy—z *Pix 
50 2817. 591 35480. 87 » 2Fou—y * Ding 5 . 212 . 945 733 . 106 1. 263? 
100 2821. 34 35433. 72 b*Hya—z *Hiy 6 95 1. O1 29 0. 98 
l 2821. 883 35426. 90 b 4Deay—z 2Pix 
150 2822. 542 S54 18. 65 a‘D z 4Pj 6 528 lL. 185 1. 713 795 1. 448 
10 2822. 778 354.15. 67 a*Dgu—z *F ix 
75 2823. 176 35410. 68 b 2Gac—y Hing 5 O05 0. 93 0. 98 19v 1.17 
l 2°23. 596 35405. 41 ¢ *Dys—t PDing - 
30 2825. 071 35386. 93 = 5 168 1. 202 1. O34 O85 1. 454 
10 2825. 162 35385. 79 
10 2825. 484 35381. 75 ¢ *Gay—at FHi 5 0. 942 0. 96 1. 02 
50 2826. 230 35372. 41 a ty z TS 
50 2826. 674 35366. 86 atG z Tin 5 278 1. 234 956 0. 140 2. 485 
2 2827. 053 35362. 12 b 4B, 2 (Hi 
2 2827. 433 35357. 36 b?D l 
2 2827. 604 35355. 23 : 
20 2829. 092 35336. 63 a *Gyu—z 2Giy 7 1. 021) 1. O1 0 0. 975 
50 2231. 84 35302. 3 b *Fug—y *P iy 5 0. S84 0. 405 1. 289 0. 442 2.615 
50 2833. 806 35277. 86 b?Hs—a 'F4 / 12 1. O4 1. 16 52lw 0. 52 
2h 2834. 232 35272. 55 b *Poay—y ? Diy 
Zh 2836. 319 35246. 60 b 4*Deuw—y *7Ghy 
3 2837. 112 35236. 75 € *Gauy—w I Fi, 
3 2837. 898 35226. 99 a*Gyc—y *F iy 1. 167) 1. 39 w 1. 28h 
10h 2838. 875 35214.87 a *Gow—y *Pix 
10 2839. 385 35208.54 bb *Dyy—ar Dix, 
50 2841. 147 35186. 71 b 2G, y *Hin 7 1. O67) 1. 315 0 1. 116 
LOO 2841. 680 35180. 11 a‘De z Pin 6 235 1. 339 1. 574 0. 592 1. 456 
75 2844. 716 35142. 56 » 2F y *Ghy 7 1. O83 1. O7 0 1. 029 
15 2845. 229 35136. 23 b *Dow—az * Dix 6 240 Q. 145 0. 385 0. 120 0. 265 
10 2847. O87 35113. 30 a ‘*Gy, 2 *Gihy 
20 2847. 574 35107. 29 a *Fa.—z *Pix 1 . 850 . 86 et $25 12 
50d 2849. 289 35086. 16 15 15 
15 2849. 591 35082. 45 b 2Gau—y 2F iy j 314 0. 883 1. 20 157 ? 
2 2850. 024 35077. 12 
20 2850. 697 35068. 83 d *Dy—w ? Dixy 5 265 0. 803 lL. OGS 135 1. 466 
5 2852. 993 35040. 62 
10 2853. 126 35038. 98 b 2Piy 1 *D3x 5 12 1. 19 1 31 Tw 1. 48 
1) 2854. 722 35019. 39 a‘’p z *Pix, 5 1. 730 0. 030 1. 700 . 865 2. 565 
15 2854. 980 35016. 23 b 2F xr G3 
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043 
327 
OBA 


29-9 


94 
19 


263? 


QS 


M09 
15 


yz 


Sh 


16 
56 


~~ 











Intensity 


30 
30 
lo 

6 
20 


200 
th 


20Hu 


20 


10) 
10 
10 


10 


9 
A (alr 
A 
PSH. 357 
9855. 710 
2856. 570 
2857. 240 
29857. 590 
Q9RAT 780 
2858. 042 
PR5S. 5OO 
2859. 622 
2859. 749 
2861 10S 
261. GSO 
2862. S48 
2863. 261 
9RO5. 544 
2866. 096 
2266. 276 
2867. 107 
2868. 832 
2870 573 
2871. 493 
2873. 329 
2873. 752 
ISTS. O56 
2879. O78 
ISSO. OD! 
2SS80. 440 
2882. 112 
S85. 473 
2887. OS4 
2RSO, 396 
2890. 476 
2891. 219 
2ROL. OST 
284. 148 
2897. 713 
2808. 242 
28090, 502 
29000. 441 
2001. 431 
2902. 026 
2902. 572 
2903. 304 
2006, 937 
2907. 160 
2909, 329 
29009, 750 
2O090. S45 
2912. 185 
2912. S04 
2913. 462 
2913. 803 
2913. 999 
2916. 370 
2917. 416 


TABLE 6. 


A 


35011. 60 
35007. 28 
34996. 74 
34988. 53 
349084. 25 


34981. 92 
34978. 72 
34973. 04 
34959. 39 
34957. 84 


34941. 23 
34934. 17 
34920. 00 
34914. 96 
34887. 15 


24880. 43 
34878. 24 
34868. 13 
34847. 16 


34826. 03 


34814. 87 
34792. 63 
34787. § 
34735. 48 
34723. 15 


34710. 94 
34706. 74 
34686. 61 
34646. 20 
34626. 87 


34599. 1 
34586. 2 
34577. 3 


34568. | 
34542. 36 


1499. 86 
$493. 57 
1478. 58 
$467. 42 
$455. 66 


34448. 59 
34442. 11 
34432. 36 
34390. 40 
34387. 76 


34313. 38 
34309. 36 
34307. 06 
34279. 17 
34266. 88 


The 
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second spectrum of ruthenium (Ru it) 
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2Dow, 


‘Dx. 


Cr gig 


‘[) 


GG 
G; 


Has, 


Dic 


! 


tion 


‘erm combina- 


‘Dix 
2D; 


Hy 
*'Cibis 
ah 
1G 


‘His 


Guy 
2Gjug 
Ging 


0 6 
Ty pe Ag 
1, 5 0. 162 
} 0. O04) 
a) 0). 240 
} 
6 1. 165 
D 0. 376 
a) . 212 
a) . 14 
6 0. O06) 
5 (0. OS) 
6 0. 12 
t 301 
5 _ 188 
i 
} . 380 
4 . 228 
t 0. OL) 
} 0. 353 
5 
i) 14 
} 0. 364 
4 0. O04) 
6 0. 196 
t 
5 . 407 
5 . 384 
t 143 
5 
1* 
6 . 261 
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. 942) 


. 205 


. 958) 


. 866 


. 658 
. 954 
Ol 


. 789 


. O88 


. 083) 
. 4038 
. O9S 


. 942 


. 397) 


. OSO 
. 102) 


. 655 
. 403 


. 965 


. 884 


. 337) 
. 06 

. 962 
. 647 


ae 


. 926 


q 


0. 98S 
. 965 


1. 03 


2. 031 


1. 034 
1. 166 
1.15 
0. 87 
1. 14 
1. 16 
1. 28 
0. 797 
1. 130 
1. O17 
1. 308 
0. 99 


1. 302 


1. 35 
0. 82 
1. 248 
1. 231 
1. 02 
1. 158 
1. 054 


1. 566 


0. 783 


1. O71 


9 


Strong p 


. O81 
0. llw 


. 120 


0. 583 
192. 
. 106 
070 


0 

0. 36 
. 27w 
. 42 
. 150 


. 093 
w 

0. 190 

. 114 


. 25w 


. ie 
w 


0. 180 


0 

0. 17 w 
. 882 
. 19w 


. 203 
. 192 


. O72 


0. 13w 
. 651 


10 


Strong n 


Nee 


. 026 


. 234 


04-4 


. 06? 
. 35 
. 072 


142 


. 426 
. 08 


. 104 
. 201 } 

















TAaBLe 6. The second spectrum of ruthenium (Ru 11)—Continued { 
” 
I 2 3 1 5 6 7 8 9 10 ] 
Intensity X (air) a Term combina- | Type Ag q q Strong p | Strong » Inte 
tion 
A K ; 
50 2918. 521 34253. 90 h 2F r 4*G3y 5 0. 04) 0. 94 0. 98 Llu 1. 09 
lh 2920. 514 34230. 53 1 
10 2920. 949 34225. 43 b?Dag—r *Ging 6 0. 402 110 0. 70 1. O0O6 0. 90 i 
3 2921. 880 34214. 53 b4*Diug—r *Dix 
lh 2922. 192 34210. 87 b 4D, > 2FS . 
25 2922. 347 34209. 06 d?Da,—w ? Dixy 6 131 1. 20 1. O7 0. 328 ? | 
10 2923. 112 34200. 11 7 0 0. 87] 
10 2923. 906 34190. 82 b?Fsu—ar *D3y 6 122 1. O86 0. 964 0. 42) 1. 025 
3H 2926. 519 34160. 29 ¢ 2? Dag—w IF yy 
200 2927. 535 34148. 44 a*Doa—z *Div 5 193 1. 337 1. 144 096 1. 626 
2 2928. 063 34142. 28 b*D r ‘Di: 
2 2928. 487 34137. 34 a °F, z ‘P 
20 2931. 209 34105. 64 b 2Fas—r 4 Fhxy 5 0. 09) 0. 97 1. O6 23w 1. 29 
3 2933. 059 34084. 13 b 4D xr Dix 5 0. 151 1. 333 1. 182 O75 1. 56? 
50 2933. 232 34082. 12 b 2Fs: > 2H j 0. 066) 1. O83 1. 02 . 23w 0. 784 
l 2933. 710 34076. 57 
10 2935. 517 34055. 59 b ?Pog—liy 5 0. 333 0. 876 1. 209 . 166 1. 375 
25 2937. 043 34037. 89 c?D w * Din. 5 242 . 794 0. 552 12] 0.915 
lh 2938. 967 34015. 61 
5 2941. 979 33980. 79 b?Dug—y *Dsy, 
50 2942, 244 33977. 73 a *Dys—z 4Fin 6 724 1. 192 168 1. O89 830 
10 2945. 098 33944. 81 e2G y * THs 7 0. 942 . 95 0 . 970 
500 2945. 661 33938. 32 a*Dy,—z *Dixy 7 1. 403) 1. 397 0 1. 400 
lh 2948. 279 33908. 18 b ?Fac—y *Gin 
l 2948. 847 33901. 65 
20 2950. 028 33888. O8 a‘*G > AH Guy 5 211 1. 231 1. 020 0. 106 2. 181 
10 2952. 246 33862. 62 b 2?F x ‘F; 5 13 1. O8 1. 21 470 1. 54 
10 2952. 706 33857. 35 | 
25 2954. O84 33841. 55 a *Gys—z 'G3xy 5 . 184 1. 167 0. 983 . 098 1. 811 
I 2954. 855 33832. 72 e*Dag—y *Giu 
25 2960. 216 33771. 45 { 0. 180 1. 367 1. 547 0. O91 0. 917 ' 
20 2061. 538 33756. 38 bh 2F zr Gx 6 . 248 0. 953 0. 705 623 830 
10H 2962. 08 33750. 20 
60 2963. 398 33735. 19 a‘Dp = Dixy 5 1. 200 1. 200 000 . 601 1. 800 
200 2965. 554 33710. 67 a‘D z IFS. 1 1. 337 1. 28 w 1. 14-4 
5 2965. 879 33706. 98 a *Tes—y Hay 
10 2966. 399 33701. 07 5 0. 18w 1. 35 
10 2968. 022 33682. 64 d?D y ?Pixg 1 0. 08 1. 21 1. 29 12u 1. 09+ 
th 2070. 672 33652. 59 b?D > *P ix, 6 0. 21 0. 965 Be 31 t 
2h 2971. 057 33648. 23 a 4G, 2 4G 5 333 0. 645 0.978 166 1. 811 | 
15 2072. 466 33632. 28 b*D r D3. 6 0. O7 1. 35 1. 28 17 1. 32 | 
10 2972. 578 33631. 02 b2G, / TEx 
8 2972. 98 33626. 47 bi Ad £37 6 0. 08 1. O8 1. 05 16 1. 06 
= h > 
Sh 2973. 826 33616.90 / A 
f ’ j 2 ¥ 
1 2974. 643 33607. 67 bh 2F 1 4G 6 0. 05) 1. OS 1. O82 19 1. 05 
1 2974. 882 33604. 97 15 0. 158 OSO ? 
200 2976. 578 33585. 82 a‘*D) Fou } 0. 05) 1.38 1. 33 1Ou 1.144 
30 2977. 219 $3578.59 a‘D ‘Fin, 5 0. 46 0.47 0. 01 220) 0. 700? 
25 2077. 471 33575. 75 b2?H ‘] }? 0. 04) 1. O8 1.12 Zlu 1. 33 
50 2978. 638 33562. 5 h2u, | 5 0. O9) 0. 93 1. 02 12 1. 45 
60 2979. 713 33550. 49 a‘D ‘T) 7 1. 337) 1. 331 0 1. 334 
70 2979. 946 $3547.87 a‘D iP { 0. 14 1. 26 1.12 0. 22u 0. 90 | 
25 2980. 966 33536. 39 b 2G, r 4G 5 0. OR) 1. 09 1. Ol Tu 1. 36—? 
10H , 2982. 026 33524. 47 h 2P r 4F; 
l 2983 751 32504. 75 a ‘Cy , IP 
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TasLe 6. The second spectrum of ruthenium (Ru u)—Continued 
f 
2 3 ! 5 6 7 s 9 10 
10 I 
‘trong , Intensity d (air o Term combina- | Type Ag g g Strong p Strong n 
Mt tion 
A K 
s 
09 5h 2085. 676 334383. 49 c Diss x *P ss 1 0. 940 0. 809 1. 749 . 471 0. 339 
+ 10 9991, 453 33418. 82 a *Gyg—z *Gbrg 5 (0. 06) 1. 18 1. 24 . 28w 1. 52 
90 70 2991. 621 33416. 95 a *Dig—z *Dixg 7 (1. 188) 1. 156 0 1. 172 
ii 15 2992. O83 33411. 79 a 4Gya—z 'Gix : 
20h 2992. 601 33406. O1 a *F3.—z *Piy } 0. 440 (1. 14) 1. 58 0. 220 ? 
° 10 2993. 142 33399. 97 b?*Deay—a *D3xy 
87] 5h 2994. 650 33383. 15 a *"su—y *Hi - 
025 20h 2996. 002 33368. 09 b *Paus—ar ‘Dixy 7, 6 (1. 188) 1. 21 0 1. 20 
Bic 1 2996. 331 33364. 42 b 2Py z *Fiu 
626 3h 2998. 126 33344. 45 b 2*Pos—z 2Pixy 
75 2998. 886 33336. 00 a*Dos—z *Dix Unaffected 0 0 
50 2999. 789 33325. 96 b*Hya—z 2F hu 1 0. 248 0. 960 1. 208 0. 124 0. 092 
29 10 3000. 465 33318. 46 b*Dys—a * Dix, 
5G? 25 3005. 141 33266. 61 ¢ *Day—a ?F ing 1 1. 17 (1. 137) w 1. 054 
7 l 3005. 891 33258. 31 ¢ *Dyy—vw fF Sy 
7 
10 3006. 716 33249. 19 7 0 1.2 
375 Qh 3007. 88 33236 32 a F 2, z *F Su 6 . 21 0. 86 (1. 07) 0. 105 ? 
915 30 3010. 450 33207. 95 a 4*Gy z Hiv 1 (0. 09) 1.13 1. 04 . 40w 0. 644 
lh 3014. 360 33164. 88 b 2*Gys—y 2Gig 
30 3015. 410 33153. 3: fb 2Hsy—z *Giy 5 (0. 06) 1. O5 0. 99 0. 29w 1. 32 
1b 2Dn y Gi =f 
aa 5 3015.956  33147.33 | b2Fag—z?Piy | 4 0.200 0951 1.151 100 0.651 
100 30 3017. 816 33126. 90 b 2Fau—y *Ging 5 (0. 06) 1. O9 1.15 . 22w 1. 37 
lh 3018. 690 33117. 31 b*Dow—y *Gix : 
l 3019. 642 33 L06. 87 b?Dnug—y *Dix 
5 3019. 76 33105. 57 a *Foa—z *Dix i 0. 270 0. 871 1. 141 . 134 0. 466 
- | 3020. 549 33096. 92 b 2Gau—y Hin, 
lh 3021. 23 33089. 47 {6 *Has—z wes | 
S11 la “Soy u Ding } 
th 3021. 977 33081. 29 a 2Asc—y *Hexy 
10h 3022. 61 33074. 36 € *Dys—2iy 6 (0. 789) 0. 85 wv 0. 82 
7 | 10h 3022. 946 33070. 68 — -b 2 Deow—y *Dix i 0. 133 1. 106 1. 239 0. 066 91 
a 8 3023. 841 33060. 90 b 2*Fy.—z 2F ix 5 . 232 1. O87 0. 855 . 117 1. 667 
g00 Sh 3026. 817 33028. 39 c 2 Dow Y 2F iu, | 
14 30 3027. 780 33017. 88 a *Dos—z *Dix 5 1. 144 0. 000 1. 444 . 573 1. 716 
| 3028. 27 33012. 53 
5 3028. 928 33005. 37 b4*D3y—y *D3xu 6 (0. O04) 1. 38 1. 34 13 1. 364 
04 5 3032. 250  32969.21 = a *Su—ar * Pix 
: 7 3032. 442 32967. 13 b 2*Foau—y ?Ghg 5 (0. 06) 0. 99 1. 05 . 15w 1. 20 
11 5 3032. 624 32965. 15 
8 3033. 911 32951. 16 a *As—y ?Giy 
9 15h 3035. 792 32930. 75 b ?Fo.—ar *Diy 7 (0. 945) 0. 95 0 0. 967 
6 75 3036. 463 32923. 47 b*Hsy—z Alix 5 (0. 07) 1. OS 1. 15 0. 36w 1. 51 
2h 3037. 116 32916. 40 a Sea—y *Dix 
15 3038. 758 32898. 61 b ?Pras—y *Dixu 5 0. 410 1. 187 0. 777 . 205 1. 392 
F lh 3039. 50 32890. 58 b *Gyy—ar 'Gig 
2 3040. 096 32884. 13 c¢ *Doay—a ‘Pi 
14 l 3042. 367 32859. 58 z *Din—e *Dox 
K)? lh 3043. 342 32849. 06 b*Hsu—z *1 ix 
_ 10 3046, 132 32818. 97 a*Dys—z *Diry 5 . 155 1. 193 1. 348 . 077 1. 580 
: lh 3046, 827 32811. 49 ¢ 2? Dyg——w t Fix 
l $047. 223 32807. 22 a 4Ga--z Hix 
. | 25 3047. 702 32802. 07 a*Dey—z *Diry, 5 (1. 337) 1. 38 w 1. 50 
9 50 3049. 246 32785. 46 b*Ha—z *1hy, 5 0. O08) 0. 98 1. 06 0. 38w 1. 44 
l 3053. 80 32736. 57 
5H 3054. 563 32728. 39 2 *Fixc—e *Da ; 
100 3056. 855 32703. 85 b?*Hs—y ‘Fix 1 (0. O8) 1. O7 1 15 0. 35w 0. 72-4 
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The second spectrum of ruthenium (Ru 11)—Continued 
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Term combina- 
tion 


a*Dy.—z *Fix 
a ‘Days Z ® Fix 
b *Darg —Yy ‘Fix 


b *Dig—y *P 3x 
b TH sis Z §Ging 
a Dor Zz ®Fixg 
é 2Gas Ir Dis 


b *Ggu—y *Ging 
b H 416 Hew 


o 
te 
~~! 
=a 
w 
2 
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5 6 
Type Ag 

5 (0. 05) 
5 0. 123 
6 (0. 04) 
5 (0. 08) 
4 0. 174 
6 (0. 05) 
4 (0. 06) 
5 0. 140 
7 0 

5 (0. 08) 
5 

5 0. 457 
“5 . 304 
5 (0.07) 
4 (0. 03) 
4 (0. 03)? 
4 0. 242 
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‘ @ e< _ 
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TasLe 6. The second spectrum of ruthenium (Ru t1)—Continued 
l 2 3 4 5 6 7 8 9 10 
Intensity » (air) o Term combina- | Type Ag g g Strong p | Stron 
tion a) ef 
1 A : 
3 3270. 155 30570. 78 | ¢?2Dey—zr 'Ging : : 
1 3272.010 | 30553.45 | a *Su—y *Dix 
10 3273. 612 30538. 50 | b 2*Pas—y *Dix : 
15 3275. 729 30518. 76 b *Dguy—z *Giy 
2 3276. 783 30508. 94 
2 3281. 193 30467. 94 ¢ *Gas—ar ‘Fix 
8 3281. 518 30464.93 | b*Dsy—y *P3y : epee. 
60 3282. 973 30451. 42 b *Dys—z #7 Ding 5 0. 367 (0. 965) 1. 332 0. 184 ? 
30 3284. 343 30438.72 = ¢ *Dis—l fy Se acnd wales 
2 3285. 762 30425.58  b 2Fy—y *F iy 
l 3286. 114 30422.32 ¢ *Day—r *Fig le 
2 3286. 840 30415. 60 | b2Dey—z *Diy a 
300 3294. 232 30347. 35 | a*Dox—z *Fixg 5us | - . 49w 1. 82— 
5 3298. 321 30309.73 a *ley—z 2Tix a 
2 3300. 115 30293. 25 
ib 3301. 14 30283. 85 J b 7H 515 2 ‘Hiss 5 (0. 04) 1. 10 1. 06 . 16w 1. 26—? 
le *#G3u—y 2Gihuy 
15 3301. 219 30283. 12 . 
2 > 4 ° 
3 3308. 102 | 3022011 {0 *Hes—z (His 
\d 2Dirg —w 4 2hy 
10 3308. 885 30212. 96 c *Gyus—y *Gig 
10 3309. 233 30209. 79 b *Foays—y *Pix 
15 3310, 283 30200. 21 | d?Doy—r Gh | 
10 3312. 831 30176. 98 6 0. 265 . 397 ? 
1 3313. O10 30175. 35 ~ ‘ , 
1 3313. 738 30168. 72 ¢ *Gyas—y *Gény : 
3 3317. 015 30138. 91 | ¢ *Gye—z *Hixg | 
3 3317. 402 30135. 40 b 2*Gay—ar * Ding 
5 3317. 591 30133. 68 § a*Dyc—z ®Fing 
3 3322. 784 30086. 59 
2 3324. 553 30070. 58 a *Dy—z °F ix 
l 3324. 800 30068. 35 9 ¢ 2?Gyus—z 2His 
25 3325. 407 30062. 86 b *Prg—y *Pig 5 0. 108 1. 189 1. 297 0. 052 1. 459-— 
l 3335. 852 | 29968. 73 a 2Isu—z ‘Tix sane 
l 3338. 329 29946. 49 bb 2Fag—z 2 Diy vale 
75 3339. 788 29933. 41 a*Dys—z ®F iy 5 . 126 1. 177 1. 303 . 064 1. 492 
l 3340. 404 29927. 89 ; 
2 3340. 789 29924.45 | a *Sup—liy, 
1 3341. 34 gogi9. 51 {> {Dawu (r 
C *Xagig Ir k 4! 
2 3342. 665 29907. 65 ¢ *Dys—r AF iy 
50 3343. 200 29902. 87 b 2Fyas—z *Ghy 6 . 10 1. OS 0. 98 . 34 1. 03+ 
100 3350. 246 29839. 98 b?Hya,—z *Hiy 7 (0. 958) . 964 0 0. 937 
8 3353. 29 29812. 89 b 2*Fyac—z *Gjy 
15 3359. 230 29760. 17 b *Gay—z 27Gh ee 
15 3359. 324 29759. 34 b 2*Fag—y Pig i . 20 0. 86 (1. 058) 0. 099? ? 
a 2B ox, Z ®F ing 
5 3359. 985 29753. 49 h *Pay-—z *Pig 3 
8 3362. 140 29734.42 a4Dog—z Fix 
2 3365. 384 29705.76 = d?Das—w *Dixg : 
40 3369. 288 29671. 34 a *Dow—z *Fing 5 1. 06 . 00 1. 06 53 1. 59 
30 3370. 061 29664.53  ¢ *Gay—y *Gh 6 0. 08) . 95 1. 03 . 27 0. 99+ 
2 3373. 472 29634. 54 b 2? Dag—y 4 Fig : 
2 3376. 624 29606. 88 a Te.—z 21h 
' 
3 3380. 356 29574.19 | b*Dyy—z *HGy? 
2° 3382. 229 29557. 82 b *Dyg—z Si 
3 3387. 963 29507. 79 b*Dag—z *H Sr ! 0. 392 (1. 344) 0. 952 0. 196 ? 
2 3389. 208 29497. 0 b *Gauy—z 2F ing 
1 3393. 904 29456. 14 b ?*Gag—z ‘Th 
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The second spectrum of ruthenium (Ru 11)—Continued 
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Type Ag g g 
5us ; <4 
6 (0. 08) 1. 07 1.15 

6 (0.06) | (0.945) | 1.00 
5 | 0.339 0.954 | 1.293 
4 (0. 04) 1. 04 1. 00 

7 (1.067) | 1.03 
6 0. 121 1.082 0.961 
6 12 (1.067) | .95__ 
6 (0. 07) 1.11 1.18 
6 (0. 05) 1. 05 1. 10 
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3757. 
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Self-Ignition Temperatures of Materials From Kinetic- 
Reaction Data 


Daniel Gross and A. F. Robertson 


Results of experimental determinations of the kinetic constants of the self-heating 
reaction are presented for wood fiberboard, cotton linters, sugar pine, cork, crepe rubber, 
GRS rubber, natural, synthetic, and blended foam rubber (with and without additive), 
various oils (raw linseed, cottonseed, rapeseed, sperm, olive, castor, and neatsfoot) applied 
to cotton gauze in a ratio of 1 part of oil to 6 parts of cotton by weight, ammonium per- 
chlorate, and nitrocellulose plastic. Under the assumption that self-heating follows a 
first-order reaction, these constants were used to calculate the critical radii of spherical 
piles for each of four surface temperatures likely to be experienced in long-period storage. 
Calculated self-ignition temperatures of piles of 's-inch-diameter to 22-inch-diameter spheres 
of wood fiberboard and \-inch-diameter to 2-inch-diameter spheres of cotton linters were 
in reasonable agreement with previous measurements by N. D. Mitchell (National Fire 
Protection Association Quarterly 45: 162, 1951). 


1. Introduction 


Previous workers studying the self-ignition be- 
havior of materials have experienced difficulties in 
the determination and_ specification of a unique 
ignition temperature. Because of the critical heat 
balance that exists during the ignition process, the 
measured “ignition temperature” is dependent upon 
the characteristics of the testing apparatus, the 
degree of control of the ambient temperature condi- 
tions, the time of exposure, and other experimental 
factors. It is highly dependent also upon the size 
of the specimen under test. This is due to the fact 
that self-ignition of a material is only possible when 
the rate of heat generation within the material ex- 
ceeds the rate at which heat is lost by radiation, 
convection, and conduction to the surroundings. 
Since the heat generated is a function of the volume 
of the material, while the heat lost is a function of 
its surface area, the measured ignition temperature 
should decrease as the linear dimensions increase. 

The incidence of fires resulting from self-ignition 
makes desirable an understanding of the relationship 
between the self-heating reaction and the physical 
parameters involved. The present paper proposes 
and illustrates the manner in which kinetic-reaction 
data obtained under conditions of negligible heat 
loss may be used to predict critical surface and 
ambient temperatures as a function of size and shape 
for a variety of materials. 


2. Analysis and Nomenclature 


In the analysis, the following symbols and units 
will be used: 


A heat-generation coefficient cal see~!em-3 


B radius or half-thickness em 

B.  eritical radius of sphere em 

c specific heat cal g-!° K-! 

e  Napierian base 

E activation energy keal mole~! 

h surface heat-transfer co- cal see~!'em-2° K-! 
efficient 

k thermal conductivity cal sec™lem=!° K-! 


m coefficient in eq (1) - 
, vA 4 
Rh gas constant 0.001987 2 keal mole o K 1 





t time sec 
T temperature within mate- ° K 
rial 


T, self-ignition (critical am- ° K 
bient) temperature 

Ts critical surface tempera- ° K 
ture 

To critical center temperature ° K 

‘m mean temperature=(To+ °K 


Tp)/2 
x space coordinate em 
p density g em- 


Analyses of self-heating have been made under the 
assumption that heat is generated within the body 
according to a first-order reaction law. Although 
this may be an oversimplification for complex mate- 
rials and over extended temperature ranges, the 
assumption appears to be applicable for many mate- 
rials in the relatively low-temperature ranges where 
self-heating becomes important and the ignition 
process is initiated. Provided the rate of the self- 
heating process within the material is not limited by 
diffusion of the reactants and/or products, or by 
depletion of the fuel supply, a steady-state condi- 
tion can exist when the heat generated within the 
material is equal to the heat lost to the surroundings. 
The particular condition of interest is the one that 
produces the maximum self-heating which is still 
steady state. Either a slightly higher ambient tem- 
perature or a slightly larger bulk size results in a 
greater quantity of heat being generated than can 
be dissipated, and ignition results. Since this steady- 
state condition is critical with respect to ignition, the 
corresponding surface temperature and size of mate- 
rial are referred to as the critical surface temperature 
and the critical size. The ambient temperature cor- 
responding to given critical steady-state conditions 
with given heat losses is considered the critical ambi- 
ent or self-ignition temperature. There exists for 
each material and for each set of heat-loss conditions 
a series of corresponding self-ignition temperatures 
and critical sizes. 

Consider a homogeneous material of half-thickness 
B generating heat in accordance with the first-order 
reaction law Aexp—E/RT and losing heat from its 
surface. The steady-state heat-conduction equation 
which applies for the condition when the heat lost 
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through the surface is equal to the heat generated 
within the material may be written: 


Yl dT e E 
I + ae RT—(), (1) 


The critical condition of ignition in the theory of 
thermal explosions originally proposed by Frank- 
Kamenetzky [1]! requires the solution of this non- 
linear equation. Chambré [2] presented analytic 
solutions in terms of known functions for the sphere 
(m=2) and cylinder (m=1) to supplement that for 
the semi-infinite slab (m=0). Recent analyses by 
Genensky [3] for a sphere and by Thomas [4] for 
all three geometries have been made employing the 
assumption common to the previous solutions, 
namely, that temperature differences within the ma- 
terial are small in comparison with the temperature. 
Enig, Shanks, and Southworth [5] have provided 
solutions in the form of convenient tables of the 
center and surface temperatures for sphere, cylinder, 
and semi-infinite slab geometries without this restric- 
tive assumption. From the tables it is useful to 
note that for a given surface temperature 7',, the 
critical radius for the cylinder and the critical half- 
thickness for the semi-infinite slab are given very 
closely by 0.775 B, and 0.514 B,, respectively, where 
B. is the critical radius for the sphere. Assuming 
that the heat generated within the material at a 
mean temperature 7’, is lost from its surface accord- 
ing to Newton’s law of cooling, the following relation 
may be used for calculating the critical ambient 
temperature for a sphere: 


where 7), is a mean temperature, here assumed as 
T+ T3/2. 

Since the analysis by Thomas employs Newtonian 
surface cooling as a boundary condition, the ambient 
temperature is obtained directly. From the graphi- 
cal relations given, it is again useful to note that for 
a given ambient temperature 7, the critical radius 
for the cylinder and the critical half-thickness for the 
semi-infinite slab are given very closely by 0.782 B, 
and 0.520 B, respectively. 

The experimental method used in this study for 
determination of the constants of the first-order 
reaction involved self-heating conditions which 
differed from those of the self-heating process 
described. Instead of a fixed ambient temperature, 
the temperature within the adiabatic test furnace 
surrounding the specimen was controlled to maintain 
negligible loss of heat. Under these conditions, the 
specimen temperature remained uniform throughout 
its mass and the heat generated within the body 
increased its temperature according to the relation: 


' Figures in brackets indicate the literature reference at the end of this paper 
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d ar 1 
. Ae FT. 
“ dt ’ 


(3) 


By transposing terms, taking logarithms and 
plotting In d7T/dt versus 1/7, the resultant line hag 
slope of —(//R)and intercepts the In d7'/dt axis ath 
Ajpe. For each material the determination of th 
values of A and £, as well as of its more commoy 
thermal properties, permits the solution of th 
general self-heating problem including the Case 
where size and temperature become critical and 
ignition results. 


3. Experimental Apparatus 


In the determination of the kinetic constants, yg 
was made of an improved version, figures 1 and 2.9 
an adiabatic furnace previously described (6) 
Although the basic operation remains the same 
refinements have been made and are here sym. 
marized: 

Addition of guard heaters and associated control 
system. In order to reduce heat losses and to maip. 
tain minimum temperature gradients within the 
furnace enclosure, the Dewar flask has been enclosed 
within a close-fitting cylindrical shell. The shel 
consists of two concentric stainless steel cylinders 
containing two electric heating elements and an 
insulating fill. The shell contains the bottom and 
center guard heating elements while a top-guard 
heating element has been placed within the Dewar 
flask plug composed of several layers of asbestos 
board. A thyratron control system as shown ip 
figure 3 permits adjustment of the guard-heater 
evele. 

Use of multiple ~junction thermocouple. The sens- 
ing element for the servo controller consists of eight 
No. 30 gage (B&S) copper-constantan  thermo- 
junctions arranged in series to measure the mean 
temperature-difference between the specimen and 
its ambient medinm. 





recording 


equipment, > 


Figure 1. Adiabatic furnace with control and 





Fi 





The changes and additions described permit im- 
proved performance of the system. The controller 


(3) can sense temperature differences of the order of 
0.01° C while the rate of temperature drift as ob- 

thms an served with an inert specimen has been reduced to 
line has, approximately 0.25° C/hr over the temperature 
t AXis at |p range 30° to 300° C. 
10. Of the 
€ commo, 4. Results 
Of the , : 
the ca 4.1. Determination of A and E 


itical and Table 1 lists some thermal and kinetic constants 
for a number of materials. Some of the values of & 
and ¢ were determined during this investigation, 
while others were obtained from handbook sources. 
The values of A and E were determined over the 
indicated temperature ranges from experimental 
temperature-versus-time data and the use of eq (3). 
For some of these materials, there is reasonable agree- 
ment with published values obtained by other 
methods. These values were determined over lim- 
ited temperature ranges before loss of weight became 
appreciable. An analysis that includes the effect of 


S 


stants, use 
1 and 2 of 
ribed (6), 
the same, / 
here sum- 


‘ed control 


to main. : . . 
rithia 7 weight loss may be applied for the higher tempera- 
n endlial ture ranges. Analysis of several cellulosic materials 

} has actually shown that when weight changes are 


The shell 
cylinders 
S and an| Fieurr 2. Adiabatic furnace, showing Dewar flask and speci- 


taken into consideration, the values of the kinetic 
constants were not much different from those 
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Figure 3. Circuit diagram of furnace-control system. 
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TABLE 1. 


Material Temperature 


range 
Cc 

Wood fiberboard 80-225 
Cotton linters 180-260 
Sugar pine 140-310 
Cork 170-225 
Crepe rubber 210-310 
GRS 220-330 
Foam rubber, 100% natural 100-180 
Foam rubber, 100% synthetic 115-210 
Foam rubber, blended, with additive 110-190 
Foam rubber, blended, without additive 110-200 
Raw linseed oil 90-120 
Cottonseed oil * 75-120 
Rapeseed oil 90-140 
Sperm oil * 120-160 
Olive oil * 115-160 
Castor oil « 140-160 
Neatsfoot oil 140-170 
Ammonium perchlorate 205-270 
Nitrocellulose plastic 135-170 


p c 
gem calg"'°K cal sec 
0. 25 0. 33 0 
32 32 
. 36 . 32 
13 4s 
v2 5 

1.10 5 
108 5 
129 5 
118 5 
104 5 
309 34 
338 34 
316 34 
313 34 
308 34 

-314 34 
316 . 34 
91 35 
S82 34 


Physical and kinetic properties 


Property 


“em Kk 
00012 
00010 
00024 
00010 
O00 282 
000299 
OO0096 


. 000100 


OO009S 
OO00095 
oOoorl 
OOOL1 


. 00011 


Ooo11 
Ooo11 
ooo11 
00011 


00051 


keal mol: 
25. 7 


t Ww -_ n~ 
SEISSsoZz Se 
SHO wu- 


nage 


wrest 


_ 
te 


cal sec 


1 


9 


NK APN Kawa & eel tS 


' 
cm-3 
97X< 10° 
30X10"! 
11x 106 
5OX 105 
06 X 108 
11x10" 
48x 10 
96 X< 10° 
25x 10 
10x10! 
SOx 10" 
42X10 
13 105 
56x 10° 
41x10" 
33X10!) 
82x 10! 
22x 10" 


3X 10!" 


sApplied to cotton gauze in ratio of 1 part of oil to 6 parts of cotton by weight. All property values refer to oil-cotton combination. 


4.2. Determination of Critical Size and Temperature 


The critical surface temperatures of spherical piles 
of wood fiberboard and cotton linters were computed, 
using the measured values of the kinetic and thermal 
constants, and the tables given in reference [5]. 
These computations were performed over the radius 
range corresponding to the specimen sizes reported 
by Mitchell [7]. Equation (2) was then used to cal- 
culate the corresponding critical ambient tempera- 
tures. The value of the surface heat-transfer 
coefficient was assumed as 4.4 10-* cal/see cm? °K. 
Using this value for the surface heat-transfer coeffi- 
cient, the critical ambient temperature was also 
determined from the relations given in reference [4]. 
The results of these calculations are shown in table 2 
for two materials for which published experimental 
data [7] are available. With the exception of the 
smallest size specimens, for which the assumption of 
an arithmetic mean temperature for eq (2) introduces 


Comparison between measured and calculated 


T, 


TABLE 2. 
self-ignition temperatures, 


Size Wood fiberboard Cotton linters 
Self-ignition temperature  Self-ignition temperature 
Thick- Half- Meas- Calculated' Mes- Caleulated> 
ness thick- ured ured 
2B. ness 
B. 
Ret7 | RefS5and) porg | Retz  Ref5and) por 
eq 2 eq 2 
in. cm c Cc c ( ( c 
1g . 159 317 183 260 304 257 307 
hy . 318 252 218 243 270 274 291 
le 635 212 215 223 248 267 272 
l 1.27 202 200 202 222 251 253 
2 2.54 180) 180 180 201 233 233 
4 5. 08 1530) 159 1) 
s 10.2 130 140 140 
12 15. 2 122 129 129 
22 27.9 109 114 1i4 


« Specimens 44, 44, and ‘4 in. thick were cubes; specimens | in, and larger 
were octagonal prisms 


b Based on spheres of radius equal to half-thickness. 





considerable error in the calculated  self-ignition 
temperature, the results are seen to be in reasonable 


agreement. 


Table 3 lists the critical radius B, of a sphere corre- 
sponding to several critical surface temperatures for 
each of the materials listed in table 1. 
been plotted in figure 4 showing this relationship for 
In this temperature range (20° to 


seven materials. 


100° C), of practical interest for ordinary storage, the 
critical size calculations were based upon the assump- 
tion of the applicability of the kinetic constants whieh 
were determined over higher temperature ranges. It 
was further assumed that the thermal properties re- 
mained constant over the whole temperature range. 
Critical dimensions for cylinders and slabs may be 
estimated from the previously mentioned constants of 


0.775 B, and 0.514B,, respectively. 


weight 
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TABLE 3. 


Material 


W ood fiberboard 
Cotton linters 
Sugar pine 

Cork 

Crepe rubber 
GRS rubber 


Foam rubber, 100° 
natural 
Foam rubber, 100°; 


synthetic 
Foam rubber, blended, 
with additive 
Foam rubber, blended, 
without additive 
Raw linseed oil * 
Cottonseed oil * 
Rapeseed oil * 
Sperm oil * 
Olive oil * 
Castor oil * 
Neatsfoot oil * 
Nitrocellulose plastic 
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All property values refer to oil-cotton combination 
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Figure 4. Computed critical radii and critical-surface tem- 
peratures for spherical piles of various materials. 
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critical surface temperature. For a finite h, the self- 
ignition temperature may be estimated from eq (2) or 
may be obtained directly from the relations presented 
in reference [4]. 


5. Conclusions 


Recent analyses of self-ignition have been made 
under the assumption that self-heating follows a 
first-order reaction law. Experimental determina- 
tions of the kinetic constants of this reaction have 
been made for a number of materials and these con- 
stants have been used to calculate self-ignition 
(critical ambient) temperatures as a function of size 
andshape. The calculated self-ignition temperatures 
are in fair agreement with published experimental 
data for two materials, wood fiberboard, and cotton 
linters. 


The experimental determinations and most of 
the computations were performed by Joseph J. 
Loftus to whom the authors are deeply appreciative. 
Thanks are due also to James S. een Jr., who 
designed the thyratron control circuits and to 
William H. Bailey, who assembled the furnace and 
made the drawings. 

An important contribution to this field has been 
made by Kazuo Akita of the Fire Research Institute, 
Tokyo, in a paper to be published in J. Chem. Soc. 
Japan, Section of Indust. Chem. 
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Some Studies of Atmospheric Transmittance on Mauna Loa 


Ralph Stair and Russell G. Johnston 


Studies are described of the spectral intensity of solar radiation and its transmission 
through the atmosphere above a new 11,140-foot-altitude station on Mauna Loa, Hawaii, 


during May and June 1957. 


microns to 2.5 microns for different air masses from about 3.0 to 1.0. 


The spectral measurements covered the region from 300 milli- 


The amounts of ozone 


and of water vapor which were present in the upper atmosphere were calculated from spec- 


tral data obtained on several different days. 


1. Introduction 


During the past few years measurements on the 
spectral distribution of the radiant energy from the 
sun have been carried out at locations in Colorado 
and New Mexico [1-3]! because of high altitude and 
creat distance from commercial activities. How- 
ever, in each case there was a considerable amount 
of interference resulting from the presence of rela- 
tively great amounts of dust or water vapor much 
of the time. To get away from conditions of this 
type it was proposed that a station be set up at a 
high altitude in the Hawaiian Islands for this work. 

‘Although few weather data were available for high 
altitude locations in the Hawaiian Islands, such data 
as existed indicated the superiority of a location on 
the north slope of Mauna Loa, on the island of 
Hawaii, for use in solar radiation and atmospheric 
transmission studies such as were being made. 
Furthermore, the available data and information 
indicated that an observing station situated at an 
altitude above 11,000 ft on this mountain would be 
above the normal atmospheric inversion layer much 
of the time and hence above most of the atmospheric 
dust and water vapor. Furthermore there are but 
few urban or rural land areas nearby to produce city 
and country air-pollution products. 

For many years the U. S. Weather Bureau had 
been interested in a high altitude weather observing 
station on Mauna Loa. The first positive action in 
this direction occurred in 1951 when, through the 
efforts of various persons and agencies, a small 
shelter was erected at an altitude of approximately 
13,450 ft near the Mokuaweoweo crater at the sum- 
mit of Mauna Loa. As this location was difficult 
of access because of a rough lava terrain and the 
absence of an improved roadway it became advisable 
to establish the Mauna Loa Observatory at the ter- 
minus of the existing road just above 11,000 ft (at 
approximately 11,140 ft). This altitude has a fur- 
ther advantage over that of the summit location in 
affording conditions of atmospheric pressure (and 
hence oxygen concentration) more in keeping with 
the normal requirements of the average individual. 
Many observers would be unable to live and work 
at a much higher level. 

Late in 1955 arrangements were set up between 
the U.S. Weather Bureau and the National Bureau 





Figure in brackets indicate the literature references at the end of this paper. 





of Standards for the latter institution to establish 
an observing station (with funds supplied by the 
U.S. Weather Bureau) on Mauna Loa, for the joint 
use of the two laboratories. Accordingly, through 
contacts with the Territorial Government of Hawaii 
a small parcel of land (4.05 acres) on the north slope 
of Mauna Loa, the land marker being set at an alti- 
tude of 11,134 ft, was transferred back to the Federal 
Government (Department of Commerce) and the 
Mauna Loa Observatory was constructed by the 
National Bureau of Standards during the early part 
of 1956. 

Since this is a new station and location for spectral 
solar radiant energy measurements a few comments 
relating to the building and its facilities are in order. 
The main section of the Observatory (see figs. 1 and 
2) consists of a one-story concrete block building (20 
by 40 ft) partitioned into a 16- by 20-ft laboratory, 
three bedrooms, and a kitchenette, all furnished for 
pleasant and comfortable living and working quar- 
ters for six observers. A concrete platform (15 by 
45 ft) on the south side, and a tower (8 by 8 ft) 
at the east end of the building, provide space for 
solar radiation measurements and for other experi- 
ments requiring outside exposures. Electrical power 
for the laboratory and housing facilities is provided 
by diesel-electric generators at 110 v, 60 eps single 
phase. Water is provided through a roof catchment, 


storage, and distribution system, while the fuel for 
the heating, cooking, and refrigeration is provided 
A radio 


through the use of butane from tanks. 





Figure 1. Mauna Loa Observatory; altitude 11,140 ft, looking 


southwest. 
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Fiaure 2. Instrumental setup on outside platform in front of 


Mauna Loa Observatory. 


transmitter and receiver linked with a similar sta- 
tion at the U. S. Weather Bureau in Hilo furnishes 
continuous two-way contact with the outside world. 
A small library, supplemented by TV and radio 
broadcasts help in furnishing recreation at evenings 
and during bad weather. Other data about the 
Mauna Loa Observatory have been published else- 
where [4]. 

The observatory is reached by a cinder road, built 
and kept in repair by the Department of Institutions, 
Territory of Hawaii, and extending a distance of 
about 45 miles from near Hilo to just above the loca- 
tion of the observatory. Through most of the dis- 
tance up the mountain the grade is gradual so that a 
stock automobile (having high clearance) can ma- 
neuver the 11,000-ft climb in less than 2 hr. It is, 
however, recommended that a 4-wheel-drive, high- 
slung vehicle be employed for general use in this area. 

In the following paragraphs the spectroradiometric 
setup is briefly described and some results are given 
on atmospheric scattering and the amounts of ozone 
and water vapor above Mauna Loa on several days 


in May and June 1957, as deduced from measure- | 


ments of solar spectral radiant energy. 


2. Instruments and Method 


The instrumental setup (see fig. 2) was essentially 
the same as employed at Sunspot, N. Mex., in 1955 
[3]. New PbS cells (and associated housings and 
electronics) having higher sensitivities and less noise 
were employed in the present case. A further im- 
provement was obtained by operating the PbS cells 
at lower temperatures (the ambient air temperature) 
rather than at a controlled elevated temperature 
as was the case in 1955. Asa result the infrared data 
could be repeated more closely from time to time dur- 
ing each day with only small sensitivity drifts as the 
air temperature slowly changed. A type—935 photo- 
tube was employed in the ultraviolet’ spectral 
measurements. 

An additional instrument for obtaining a measure 
of the amount of sky scattering was incorporated 
into the equipment. This consisted of a photoelec- 
tric photometer for comparative brightness measure- 





| 
| 
| 
| 








ment of the sun and an area of the sky nearby and 
surrounding the sun. Readings were alternate} 
taken for the intensity of the solar disk and for that | 
of an arbitrary area of the sky consisting of an anny 
lar ring of a few degrees immediately surrounding the 
sun. In each case the radiant energy was confined 
essentially to a narrow band in the Vellow-grogy 
region of the visible spectrum through the use of g 
yellow-green filter (a dark-shade welding glass). 

The electronics of the instrument for the measure. 
ment of sky scattering were arranged so that the samp 
tuned amplifier and recorder could be readily shifted 
from the output of the spectroradiometer phototube 
to that of the auxiliary instrument without disrupting 
the recorder time scale. Thus the two records a 
peared on the same chart in proper sequence. Cali. 
bration was accomplished through an on-the-spot 
comparison with an Evans [5] visual photometer 
Thus the portrayed data are given in terms of sky 
brightness within an annular conical angle near the 
sun relative to that of the solar disk 

The use of electrical power from a small diesel. 
electric generating plant caused a number of prob. 
lems in the practical operation of the equipment, 
First, it is very difficult to keep a small diesel-electric 
generating plant (in this case a 5-kw unit) in continu. 
ous operation at precisely 60 cps with the voltage 
constant, especially at this altitude and with consid- 
erable variation in the load. The power require- 





ments of all other station activities were dependent 
upon the same power source. Furthermore, the need 
for regular servicing of the equipment required the 
alternate use of two similar units, which usually 
varied in output voltage and frequency. Because 
these difficulties had been anticipated previous to the 
setting up of the equipment, crystal-controlled 
electronic power supplies had been procured to handle 
the power requirements of all motors and amplifiers 
of the setup in which either voltage or frequency 
changes might result in erratic operation. Greatly 
improved operation resulted through the use of the 


crystal-controlled electronic power units. However, | 


there remains no question but that commercial power 
would be preferred. Hence, a real need at the Mauna 
Loa Observatory is an extension of power lines for 
the remaining 12 miles up the mountain to the 
observatory. 

In addition to the equipment employed in the 


spectral radiant-energy and sky-scattering measure- | 


ments the standard U. S. Weather Bureau station 
instruments were available and were regularly read 
and the data transmitted to the Hilo office, and in- 
cidentally employed in the evaluation of the surface- 
humidity values reported herein. 


3. Spectral Solar-Energy Distribution 


On those days when the sky was clear, during 4 
period of about 8 weeks in May and June 195/, 
spectral data were obtained either through the 
shorter-wavelength region, 300 to 550 mu, by using 
a 935-type photoemission tube, or else through the 
enete Lol wae region, about 330 my to 2.5 4, 
by using a PbS cell. These data were evaluated in 
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a manner similar to that previously followed and 
were found to be in general agreement with the 
earlier data [3] both as to spectral quality and abso- 
jute magnitude, except that they show a_ higher 
degree of scattering as a result of difficulties in 
keeping the power supply constant, or as the result 
of changes in atmospheric humidity. 


4. Atmospheric Transmittance and Ozone 


Following the same procedures as employed in 
evious work [1-3] the atmospheric transmission 
for the atmosphere above Mauna Loa (see 
figs. 3 and 4) was determined on the basis of two 
sets of ozone transmission coefficients—those of 
Fabry and Buisson |6], and those of Vigroux [7] 
for the average data for 4 days in May (May 21, 
98, 29, and 30). 

In figure 3 the Fabry and Buisson ozone trans- 
mission coefficients were employed in the calculation 
of the two curves representing the amounts of ozone 
required for a specific amount of optical absorption. 
On the basis of these data, about 0.24 em of ozone 
(ntp) was present above Mauna Loa at the time of 
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Figure 3. Atmospheric transmittance above Mauna Loa Ob- 


servatory, altitude 11,140 ft; also determination of total ozone 
using Fabry and Buisson transmission coefficients. 


Upper curve, scattering; lower curves, ozone; circles, measured atmospheric 
transmittance; air mass, 1.00; wavelength seale expanded as a function of 
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FicurE 4. Atmospheric transmittance above Mauna LoaYOb- 
servatory, altitude 11,140 ft; also determination of total ozone 
using Vigroux transmission coefficients. 


Circles, measured atmospheric transmittance; air mass, 1.00. 


the measurements. This value is in close agreement 
with the results of previous measurements of high 
altitude ozone at Climax, Colo. [1], and at Sunspot, 
N. Mex. [2, 3]. , 

In figure 4 the same experimental data (for wave- 
lengths shorter than about 337 my) are replotted, 
this time along with ozone curves calculated on the 
basis of Vigroux’s transmission coefficients. New 
data by Inn and Tanaka [23] agree closely with those 
of Vigroux within the spectral region of 300 to 350 
mu. It is to be noted on the basis of these data 
that the calculated amount of ozone (ntp) becomes 
about 0.32 em. Furthermore, if the data obtained 
at Climax and Sunspot are evaluated in terms of 
Vigroux’s coefficients of ozone absorption, then the 
previously published ozone values should be in- 
creased to about 0.30 em. Since the calculated 
value of atmospheric ozone depends upon the lab- 
oratory measurements of the absorption coefficients 
of ozone, and since there is considerable spread 
between the values obtained through using the two 
sets of data, it seems advisable to report both values 
until further experiment definitely determines the 
absorption coefficients of ozone. 
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5. Water Vapor Concentration Above 
Mauna Loa 


In previous work at Sunspot, N. Mex. [2], changes 
in upper-air water-vapor concentrations resulted in 
marked absorption variations in the infrared solar 
spectrum (see fig. 5). As a matter of fact the large 
fluctuations noted at Sunspot were a primary reason 
for seeking the new location on Mauna Loa where a 
lower water-vapor concentration would interfere less 
with the measurement of the spectral intensities in 
the infrared. But water vapor is such a_ highly 
absorbing constituent of the atmosphere that even 
at the very low concentrations ordinarily present | 
above Mauna Loa there is still a marked absorption 
within the principal water-vapor bands, even near 
noon on days having a relatively low total water- 
vapor concentration. Then small changes in water 
vapor produce large changes in total optical absorp- 
tion within the regions of the principal absorption 
bands. 

The spectroscopic method of evaluating the 
amount of the water vapor was originated by Fowle 
nearly a half century ago [8-11]. However, its use 
in routine studies has been greatly neglected until 
recent years when the U.S. Weather Bureau [12-17] 
and others [18, 19] have set up various types of 
hygrometers, filter radiometers, and spectroradiom- 
eters for water-vapor studies. 


| 
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It should be noted that the water-vapor molecyl 
has three principal modes of vibration resultin i 
three primary absorption bands located at 2.66 2.74 
and 6.26 uw [20]. The resulting combinations and 
overtones are responsible for the production of many 
bands (of lesser absorption) extending to shorts 
wavelengths. 
centered at 0.942, 1.135, 1.370, and 1.872 yg. Which 
increase in intensity of absorption with wavelength 
Each of these bands produces a marked amount of 
optical absorption in the solar spectrum. 

The relationship between the apparent absorption 
of water vapor (as measured with a prism instry. 
ment) and the amount of the vapor is not a simple 
one since each absorption band really consists of 
many fine lines of uneven intensity irregularly spaced, 
Furthermore, the attenuation at any wavelength dye 
to water vapor is a complicated function of the pres. 
sure, temperature, and concentration of the vapor 
per unit volume. 
plications does not prevent the use of the absorption. 
spectra method through the adoption of an empirical 
relationship between the optical absorption and the 
absolute amount of water in the atmosphere. 

The spectroradiometric method for determining 
the water-vapor concentration in the upper atmos- 
phere is rapid and practical for routine work—so 
long as the sun is not clouded over. It remains the 
most applicable of some 15 methods of water-vapor 
or moisture evaluation recently described [21]. This 
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Ficure 5. 


July 15, 1953, approximately 10:00 a. m. 
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Representative data taken from a section of the recorder chart while employing a PbS cell as a detector at Sunspot, 


Ne w Me Tico, 


The numbers on the chart have no significance except to identify the curves and 
data illustrated show the relative intensities of some of the principal water-vapor absorption bands 


in 1958. 


umplifier sensitivity settings. The 
The time interval represented by the data approximates 18 min. 


Principally among these are the bands | 
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as developed by Fowle, permits calibration 
f total band absorption or of absorp- 
at a single wavelength. In practice a combina- 
he two (the depth of a prismatic band) is 
has been found to yield more consistent 
18] when a prism spectroradiometer is 


method, 
either in terms 0 


tion 
tion of t 
simpler and 
results [10, 


employed ee 
Any one of the water-vapor bands may be chosen 


for use in determining the total amount of precipi- 
table water above the observer. In the present case, 
as noted above, the available bands for this purpose 
are located at 0.94, 1.135, 1.370, and 1.872 uw. (See 
figs. 5 and 6 which are representative of field data 
obtained at Sunspot, N. Mex. and on Mauna Loa, 
respectively.) It is important for highest accuracy 
to choose that band which is located in the part of 
the spectrum where the recorded intensities of solar 
spectral radiant energy are such that, in combination 
with the response properties of the spectroradiometer 
employed, a nearly horizontal curve results when the 
envelope of the recorder trace is drawn. Further- 
more, greater sensitivity in detecting changes in the 
water-vapor content of the upper atmosphere can be 
obtained if the maximum of absorption falls between 
25 and 75 percent. In the present case the band at 
1.135 u best meets these requirements. Hence this 
band has been chosen for use in evaluating the 
precipitable water above Mauna Loa. 

Although the spectroradiometric method is rapid 
and highly practicable, its use entails certain assump- 
tions relating to the characteristics of the shapes of 
the water-vapor absorption bands as affected by 








such things in the environment and instrument de- 
sign as temperature, pressure, slit width, ete., which 
may and have been studied in the laboratory [8, 10, 
11, 19], and their effects evaluated. Principal among 
these in importance is instrumental slit width or 
spectral purity. Fowle [10] determined a calibration 
employing specific slit widths through the use of a 
spectrobolometer by observing the depths of the 
water bands resulting from the absorption of a known 
amount of water vapor in a long tube. He further- 
more determined that under the usual conditions of 
temperature and pressure of the upper atmosphere 
that in the practical case their effects could be 
neglected. This calibration of the total precipitable 
water in centimeters is given in figure 7 in terms of 
the intensity ratio for the water-vapor band at 1.135 
» When the spectral purity (slit width in wavelength) 
was 0.04 uw, as employed in the Mauna Loa work. 
We have used this calibration in the evaluation of 
the Mauna Loa data. 

While this calibration and its application to the 
field data may not be of the highest accuracy owing 
to some question as to the total width of the band 
and its resultant effect upon depth measurements 
of the center of the band, it appears that the result- 
ing errors are small since they would appear in a 
similar manner in both the field and laboratory data 
and would thus for the greater part cancel out. 

The data for 4 days (June 7, 8, 12, and 13) have 
been analyzed and the amounts of precipitable water 
above the station plotted in figures 8 to 11, inclusive, 
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Ficure 6. Representative data taken from a section of the recorder chart while employing a PbS cell as a detector at Mauna Loa 
in 1957. 

_ June 12,1957, approximately 10:00 a. m. The numbers on the chart have no significance except to identify the curves and amplifier sensitivity settings. 

They should not be compared directly with those of figure 5 because of possible differences in spectrometer adjustment and the use of a different PbS cell.) The 


data illustrated show the relative intensities of some of the principal water-vapor absorption bands. The time interval represented by the data approximates 18 
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of the time of day. On these same 


- a function - : 
yoo are illustrated the relative amounts of pre- 


ater (absolute surface humidity) based 
ments of relative humidity and tempera- 


cipitable w 
= oe dhe station. These data are necessarily on 
= relative scale (but the same relative scale for the 
4 days) as no means were available for determining 
the water at higher altitudes by this method. It is 
sionificant that on any one day there is a general 
correlation between the water-measurements at the 
station and the total amount of water as integrated 
for all altitudes by the spectroradiometric method. 
This probably indicates that the major part of the 
water above the station Is near the station level. 

A study of figures 8 to 11, coupled with the gen- 
eral observations of the sky at the station over a 
riod of several months and with the experience of 
other observers in Martian studies [22] at this sta- 
tion, signifies that as a general rule there is a very 
low water concentration above Mauna Loa from 
late afternoon until late in the morning of the next 
dav. Under certain weather conditions the low 
water concentration continues undisturbed through- 
out the 24-hr period. On many days, however, the 
inversion layer is disrupted near noon or in early 
afternoon thus allowing a mass of humid atmosphere 
to sweep up the mountain and enshroud the station 
either temporarily for an hour or two or until late 
afternoon or early evening. Except during unfa- 
vorable weather, however, the mass of humid atmos- 
phere again descends the mountain during the late 
afternoon or early evening so that a picturesque 
sunset may be viewed—and on especially favorable 
oceasions in such a way that the brilliant “green 
flash” may be observed immediately following the 
disappearance of the edge of the sun from view. 

The amount of atmospheric scattering in the 
visible spectrum as measured in the annular conical 
angle surrounding the sun (figs. 8 to 11, inclusive) 
seems to have no definite relationship to the total 
precipitable water above the station. This is in- 
teresting and indicates the importance of additional 
study of atmospheric transmittance and scattering 
at high altitudes. It was noted that the sky ap- 
peared less blue when there was a minimum amount 
of water present; also that the sky usually appeared 
much bluer when viewed at the 5000- to 8000-ft 
levels than when viewed from the altitude of the 
Mauna Loa Observatory (11,140 ft). 
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Central Notations for the Revised ISCC-NBS 
Color-Name Blocks 


Kenneth L. Kelly 


Nickerson and Newhall published, in 1941, the central notations of the original ISCC- 
NBs (Inter-Society Color Council-National Bureau of Standards) color-name blocks which 
were used in the preparation of the soil and rock color-name charts. In 1955, the ISCC-— 
NBS color-name blocks were revised to accord more closely with usage in the textile and 
other fields (NBS Circular 553). The central notations of these revised color-name blocks 
have been computed and are given in the present paper in tabular form. A color chart 
showing the central colors of the 267 ISCC-NBS color-name blocks would serve for rapid 
determination of the ISCC—NBS color designation, especially in field work where speed and 
ease of operation are more important than high accuracy. 


Journal oO 


hd i 4 ° . . 
1. Introduction | Che chroma, C, of the centroid [1] is defined as: 
8 
In 1941 Nickerson and Newhall published a | 2 C. 
paper [7|’ entitled Central Notations for ISCC-NBS | | C? (cos 0)dCdé 
Color Names based on the first description of the a = oO 
[ISCC-NBS method of describing colors published a ee — 


in 1939 [2]. Since then, the ISCC-—NBS color-name 
charts have been revised to accord more closely with 
usage in the textile and other fields and have been 
published in 1955 under the title The ISCC—-NBS 
Method of Designating Colors and a Dictionary of 
Color Names [3]. The present paper contains the 
central notations of the revised color-name blocks. 


2 "C: 

| | CdCc6 

6 Ca 

7_2(C2+C,C,+€% sin (rAH/100) 

3(C,+C,) (#4 H/100) 

where C, and C, are the upper and lower chroma 

limits of the block, respectively. er: 
By this method, the central notation (H, V/C) of 

an elementary block like moderate red, hue range 

1R to 7h, value range from 3.5 to 5.5, and chroma 

range from 7 to 11, is determined as: 


, 


2. Computation of Centroids 


The color-name blocks into which the color solid 
has been divided are bounded at top and bottom by 
planes of constant Munsell value, at two of the sides 
by planes of constant Munsell hue, and at the other | H=(1R+7R)/2=4R; V=(3.5+5.5)/2=4.5, 
two sides by concentric cylindrical surfaces of constant 


Munsell chroma. The shape of such a block might | and 
be called a sector of a right cylindrical annulus (like 7 2(11°7+-7 *11+-7?) sin (67/100) 


a piece of pie with the point bitten off). The 
ns er — is ie — ni: nga 2247 X 0.1874 
two-value, two-hue, and two-chroma limits; the | =- : —— = 
most complicated of these blocks may be broken 3X18 X 0.1885 
down into three blocks of this simplest or elementary | |The central notation of a compound block, that is, 
shape. The general method is to find the centroid | one in which two or more elementary blocks are 
of each such block and then find the volume-weighted | joined together along common constant-hue planes, 
centroid of the set of one to three such elementary | such as yellowish white or moderate blue, is found 
centroids. by computing the volume-weighted mean of the 
The hue and value of such an elementary centroid | central notations of the simple blocks. The centroid 
are simply the averages of the hue and value limits, | of a two-component block, for instance, will lie on 
respectively. The chroma of the centroid is not | the straight line joining the centroids of the com- 
entirely a function of the chroma limits, but depends | ponent blocks. Therefore, it is necessary to graph 
on the number of Munsell-hue steps AH in the block. | the respective centroids on circular graph paper, 
Let the cross section of the block at constant Munsell | connect them by a straight line divided inversely 
value be laid out in cylindrical coordinates with the | according to the two volumes, and read off the hue 
radius being Munsell chroma, C, and since there are | and chroma of the resultant centroid. For a three- 
100 Munsell-hue steps in the complete hue circuit, the | component block, the resultant centroid for two of 
angle A@ in radians will be given in terms of Munsell | the adjacent blocks is found and then combined by 
hue, 7, as 27A/7/100. Further let 7=0 correspond | the same method with the centroid of the third 
to the Munsell hue, H. of the centroid of the block. block. The value notation of the resultant centroid 
J is computed as the volume-weighted mean of the 
Figures in brackets indicate the literature references at the end of this paper. | Value notations of the separate centroids. 


3(11+7) (62/100) 


9.1 
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3. Peripheral Color-Name Blocks 


The shape of the psychological color solid may be 
described as a grossly misshapen grapefruit, set so 
that the pithy core is vertical with an irregular 
outer surface or skin, determined by the MacAdam 
theoretical pigment limits [4]. The bounding sur- 
face has been determined in Munsell renotation 
terms by Nickerson and Newhall [8]. This solid, 
cut into 267 blocks each with a color name attached, 
forms the ISCC-NBS system of colornames. Now 
consider the outermost or peripheral blocks, those 
that extend to the outer surface or skin. These 
represent colors of maximum saturation, maximum 
lightness or minimum lightness as determined from 
the MacAdam limits. These peripheral blocks, of 
which there are 120, have several simple boundary 
surfaces (plane and cylindrical) and are closed by 
the outer complicated curved surface of the color 
solid. The centroids of these peripheral blocks have 
been estimated graphically by plotting the Mac- 
Adam limits on the color-name charts. 


4. Table of Central Notations of Color- 
Name Blocks 


In table 1 the ISCC-—NBS color-name blocks are 
identified by number and name, and the Munsell 
renotation of each centroid is given. The central 
notation of each of the 120 peripheral blocks, less 
certain because it was estimated rather than com- 
puted by formula, is marked with an asterisk. 
The numbers used to identify and order the color 
designations in table 1 are taken from table 1 of 
NBS Circular 553 [3]. Thus, for any hue name, the 
order of modifiers is vivid; brilliant, strong, deep, 
very deep; very light, light, moderate, dark, very 
dark; very pale (very light grayish), pale (light 
grayish), gravish, dark grayish, blackish; -ish white, 
ight -ish gray; -ish gray, dark -ish gray, and -ish 
black. Not all of these modifiers are required with 
every hue name, as for instance there is no very 
dark pink or very pale olive. The notations of the 
centroids are given to one decimal place except that 
the value or chroma notations of the centroids of 
some of the peripheral blocks are given as whole 
re rs followed by a plus or minus sign, such as 
4.5 YR 6.8/16+, indicating uncertainty in the next 
decimal place. 


TABLE 1. Central notations of ISCC-NBS color-name blocks. 
C553 
num- ISCC—NBS color name Central notation 
ber 


1 | vivid pink 1.5R 7 /13+* 
2 | strong pink 1.5R 7.5/9.1* 
3 | deep pink LOR 6.0/11.1 
4 light pink 2.5R 8.6/5.2* 
5 | moderate pink 2.5R 7.2/5.2 
6 | dark pink 2.5R 6.0/6.0 
7 | pale pink 2.5R 8.8/2.3* 
8 | grayish pink 2.5R 7.2/2.3 
9 | pinkish white 7.6R 9.2/1.0* 
10 | pinkish gray 7.6R 7.5/1.0 
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TARLE 





| C553 
} num- 
| ber 


70 


notations 


blocks. 


1. Central 


ISCC-—NBS color name 


vivid red 
strong red 
deep red 

very deep red 
moderate red 


dark red 

very dark red 
light gravish red 
grayish red 
dark grayish red 


blackish red 
reddish gray 

dark reddish gray 
reddish black 

vivid yellowish pink 


strong yellowish pink 
deep yellowish pink 
light yellowish pink 
moderate yellowish pink 
dark yellowish pink 


pale yellowish pink 
gravish yellowish pink 
brownish pink 

vivid reddish orange 
strong reddish orange 


deep reddish orange 
moderate reddish orange 
dark reddish orange 
grayish reddish orange 
strong reddish brown 


deep reddish brown 

light reddish brown 
moderate reddish brown 
dark reddish brown 

light grayish reddish brown 


gravish reddish brown 

dark grayish reddish brown 
vivid orange 

brilliant orange 

strong orange 


deep orange 
light orange 
moderate orange 
brownish orange 
strong brown 


deep brown 

light brown 
moderate brown 
dark brown 

light grayish brown 


grayish brown 
dark grayish brown 
light brownish gray 
brownish gray 
brownish black 


vivid orange yellow 
brilliant orange yellow 
strong orange yellow 
deep orange yellow 
light orange yellow 


of ISCC 
Continued 


NBS 


Centrs 


5.0R 
4.0R 


1.5R 5 
3.5R 
4.9R 
4.6R 
3.5R : 


OP wat 


~- 


3.5R 


6.5R 5 


6.0R 3 


6.0R 04 


8.0R 7 


8.0R 7.5 
5.6R 6. 


O.7YR 
O.7YR 


7.0R 6. 


O.5R 4 


O.5YR 
O.5YR 


O.S8R 1. 
O.5YR: 
V.6R 3. 


O.5R 1 


3.0YR : 


9.9R 3. 


9.5R 2 
4.5YR 
4.5YR 
£$5YVR 


5YR 5. 
4.5YR § 
OYR ! 
5.OYR 3.! 


4 
4 
L6VR 
} 


rrooc 


_ 
>) 
nates 


6Y 
5) 
SYR 
6.4Y 
5) 


VOYR 
9OYR 
VOYR 
9.0OYR 
90YR 


4 
5.0R 2. 

I 

} 


a to 


OY 
5Y 
5YR: 
5Y 
y 


color-name 


> 


| 
il notation 


9.0* 
0/12.4 
8.6/4.8* 
7.2/4.8 
0/6.0 


.0/12.0 
O.5DR 5.5 
9.5R 4. 


9.1 
0.9.1 
5.5/6.0 
3.0/11.2* 


+ 


CIO ee SIU 


4/2. 
0/2. 
6.6 
S.( 

6.5 


0.8/0.8* 
7.2/16 
8.4/1: 
7.2/1% 
6.0/1: 
8.6/8. 
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lor-name 
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Otation | 

















TABLE 


C553 


num- 


| ber 


96 


OS 
99 
100 


101 
102 
103 
104 
105 


106 
107 
108 
109 
110 


Ce Oho 


116 


118 
119 


t 
= 
— 


bo bo bo bo bo 
ee 


126 
127 
128 
129 
130 


1. Central notations of ISCC-NBS color-name 


blocks. 


ISCC—NBS color name 


moderate orange yellow 
dark orange yellow 

pale orange yellow 
strong ve lowish brown 
deep yellowish brown 


light yellowish brown 

moderate yellowish brown 

dark yellowish brown 

light grayish yellowish 
brown 

grayish yellowish brown 


dark grayish yellow brown 
vivid yellow 

brilliant yellow 

strong yellow 

deep yellow 


light yellow 
moderate yellow 
dark yellow 
pale yellow 
gravish vellow 


dark grayish yellow 
yellowish white 
yellowish gray 

light olive brown 
moderate olive brown 


dark olive brown 

vivid greenish yellow 
brilliant greenish yellow 
strong greenish yellow 
deep greenish yellow 


light greenish yellow 
moderate greenish yellow 
dark greenish vellow 

pale greenish yellow 
gravish greenish vellow 


light olive 
moderate olive 
dark olive 

light grayish olive 
grayish olive 


dark gravish olive 
light olive gray 
olive gray 

olive black 

vivid yellow green 


brilliant vellow green 
strong vellow green 
deep vellow green 
light vellow green 
moderate yellow green 


pale vellow green 
grayish vellow green 
strong olive green 
deep olive green 
moderate olive green 


dark olive green 

grayish olive green 

dark grayish olive green 
vivid yellowish green 
brilliant yellowish green 


Continued 


Central notation 


9. 5YR 2 2. 


_——_ 
_ 
—_ 


wosI6 6 
vit 


ne 
wt ee 


es ated ee 
[Re ole ole « 


S 
to S 
= oa 
Sc = 
— 
S So 
- 


“100 G2 =100 
ead 
oS 
= 
S 
san 
C 


te 
_ 


yp 
— 
or 
1 
ge 
* 


yD 
S 
—_ 
ne WO 
o 
eo" 
o 
+ 


° 2 

~ 

—_ 
wet. 


non ~1d 
NmwNAT 


wo 
- 


tho 
~ 
—) 
—, 
=> 


- 
— 
— 
~ 
2 
— am 2.9 or be 
we Oo 
@i am Pee On 
“IS Oren 
ono- > 
' SS me ee tS 
~ mm =o 2 
—- —Ocw 
_ * 
* or 
4 
* 


~ 
— 
~ 
“2 
_ 

‘ 
~ 
S 
i, 
og 
= 


5. 0G Y 4, 0/4 
5.0GY 8.6/ 
5.0GY 6.0 


3.5GY 8.7 
4.4GY 6.1 
5.0GY 3.: 


NS NS me ee 
* 


WINS ¢ 
* 


or site 
No aj More 
cin 
i @ 


» 
5.0GY 2.4+-/7.1+* 
5.0GY 3.5 
».0GY 1.8/3.7* 
».0GY 3.5/2.2 
r.0GY 2.0 2.0 | 
0.5G 7 /15.5+ 


0.5G 8.0/9.1* 
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TABLE l. 


C553 
num- 
ber 


ps fh ft fh pet 
WAnAWwW 
of Aoi 


136 
137 
138 
139 
140 


ST ee 


146 
147 
148 
149 
150 


or on ot on on 
OFA 


ourc 
GO -1 


159 
160 


161 
162 
163 
164 
165 


166 
167 
168 
169 


— 
~I 
~ 

— 


sho 


“IsJ) se] s+) *] 
CBS ore Ww 


~IsJ e736) 


—_ 
_ 


80 


181 
182 
183 
184 
185 


186 
187 
188 
189 
190 


ISCC-—NBS color name 


Central notations of ISCC-—-N BC color-name 
blocks.— Continued 


Central notation 





strong yellowish green 


| deep yellowish green 


very deep yellowish green 
very light yellowish green 
light yellowish green 


moderate yellowish green 
dark yellowish green 

very dark yellowish green 
vivid green 

brilliant green 


strong green 
deep green 
very light green 
light green 
moderate green 


dark green 

very dark green 
very pale green 
pale green 
grayish green 


dark grayish green 
blackish green 
greenish white 
light greenish gray 
greenish gray 


dark greenish gray 
greenish black 

vivid bluish green 
brilliant bluish green 
strong bluish green 


deep bluish green 
very light bluish green 
light bluish green 
moderate bluish green 
dark bluish green 


very dark bluish green 
vivid greenish blue 
brilliant greenish blue 
strong greenish blue 
deep greenish blue 


very light greenish blue 
light greenish blue 
moderate greenish blue 
dark greenish blue 
very dark greenish blue 


vivid blue 
brilliant blue 
strong blue 
deep blue 

very light blue 


light blue 
moderate blue 
dark blue 
very pale blue 
pale blue 


grayish blue 
dark grayish blue 
blackish blue 
bluish white 

light bluish grav 


0.54 5.5/9.1 | 
0.54 3.5/11.5+* | 


| 0.5G 2.0/9.2+-* 


0.5G 9.1/5.1* 
0.5@ 7.5/5.1 


0.54 3.5/5.1 


| 0.5 1.5/4.9* 
6.5G@ 5.2/18+* 


| 6.06 8.6/5.1* 
6.0G 6.5/5.1 
| 6.04 4.5/5.1 


4.5B 6.5/5.2 
4.5B 4.5/5.2 
4.5B 2.7/5.0 
4.5B 1.3/4.9* 
| 3.0PB 4.2/15+* 
| 8.0PB 6.4/11.0* 
3.0PB 4.2/11.0 


0.9PB 4.2/3.0 


| 
| 
0.54 5.5/5.1 | 


6.0G 7.4/9.1* 


6.0 4.5/9.1 | 
6.04 2.3/9.1* | 


6.0G 2.7/5.0 
6.0G 1.3/4.9* 
9.0G 8.7/1.8* 

9.04 6.5/1.8 

9.0G 4.5/1.8 


9.04 2.8/1.6 
9.0G 1.1/1.4* 
7.0G 9.2/0.8* 
7.04 7.5/0.8 

7.04 5.5/0.8 


7.04 3.5/0.8 

7.54 0.9/0.7* 
4.5BG 5.3/15.5+-* 
4.5BG 7.3/9.0* 
4.5BG 4.5/9.0 


4.5BG 2.3/9.0* 
4.5BG 8.5/5.0* 
4.5BG 6.5/5.0 
4.5BG 4.5/5.0 
4.5BG 2.7/5.0 


4.5BG@ 1.3/4.9* 
4.5B 5.2/13+* 
4.5B 6.8/9.0* 
4.5B 4.5/9.0 
4.5B 2.5/9.0* 


4.5B 8.2/5.2* 


3.0PB 2.0/9.0* 
2.4PB 8.2/7.2* 


2.7PB 6.5/7.3 
2.9PB 4.2/7.2 
2.8PB 1.7/5.0* 
0.9PB 8.4/3.0* 
0.9PB 6.5/3.0 


9.5B 2.6/1.9 
9.5B 1.1/1.5* 
9.5B 9.1/1.0* 
9.5B 7.5/1.0 











TABLE 1. Central notations of ISCC-NBS color-name 
blocks——Continued 

C553 

num- ISCC-—N BS color name Central notation 
ber 
191 | bluish gray 9.5B 5.5/1.0 
192 | dark bluish gray 9.5B 3.5/1.0 
193 | bluish black 9.5B 1.0/0.7* 
194 | vivid purplish blue 8.0PB 2.5/22+ * 
195 brilliant purplish blue 8.0PB 6.0/10.8* 
196 strong purplish blue 8.0PB 4.1/11.4 
197 | deep purplish blue 8.0PB 1.6/9.1* 
198 | very light purplish blue 7. PB 7.8+-/6.1+* 
199 | light purplish blue 7.5PB 6.0/6.8 
200 | moderate purplish blue 8.0PB 3.5/6.9 
201 dark purplish blue 7.9PB 1.1/4.7* 
202 ~-very pale purplish blue 7.0PB 8.2/4.1* 
203 pale purplish blue 7.0PB 6.0/4.1 
204 grayish purplish blue 7.2PB 3.3/4.0 
205 vivid violet 1.0P 2.7/21+ * 
206 brilliant violet 1OP 5.9/11.1* 
207 strong violet 1L.0P 3.5/11.1 
208 deep violet LOP 1.4/10.3* 
209 very light violet 1.0P 7.9/7* 

210 | light violet 1.0P 6.0/7.2 

211 | moderate violet 1OP 3.5/7.2 

212. dark violet 1L.0P 1.3/5.0* 
213 | very pale violet 1.0P 8.2/4.1* 
214 pale violet L.OP 6.0/4.1 

215 | grayish violet 1.0P 3.3/4.1 

216 | vivid purple 6.0P 3.5/20+ * 
217 brilliant purple 6.0P 6.7/11.1* 
218 strong purple 6.0P 4.5/11.1 
219 | deep purple 6.0P 2.8/10.2 
220 «very deep purple 6.0P 1.2/10.2* 
221 ~=very light purple 6.0P 8.0/7.2* 
222 | light purple 6.0P 6.5/7/2 

223 moderate purple 6.0P 4.5/7.2 

224 | dark purple 6.0P 2.7/5.1 

225 very dark purple 6.0P 1.1/4.9* 
226 | very pale purple 5.4P 8.4/3.3* 
227  ~=pale purple 6.8P 6.4/3.0 

228 = grayish purple 7.8P 4.5/2.9 

229 dark grayish purple 10.0P 2.8/2.0 
230 | blackish purple 10.0P 1.0/1.4* 

5. Applications 
For nonborderline colors it is much easier to 


determine the ISCC-NBS color name from a color 
chart containing 267 samples, one for each of the 
central notations of the 267 ISCC—NBS color-name 
blocks, than it is to determine the Munsell renotation 
of the color and then to find the ISCC-NBS designa- 
tion from the color-name charts. In using such a 
color chart, the observer simply picks out the 
particular color most closely duplicating that of the 
specimen to be designated; he does not have to 
check for conformity to hue, lightness, and satura- 
tion limits separately. The production of such a 
color chart is now under consideration. It should 
be pointed out, however, that if the specimen color 
falls about equally near to two or more of the 267 
central colors, it would be impossible to determine 
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| of the centroid locations. 








TABLE 1. Central notations of ISCC-NBS color-name 
blocks.—Continued 
C553 
num- ISCC-—NBS color name Central hOtations 
ber "a 
231 purplish white 9.0P 9.1/1.0% 
232 = light purplish gray 9.0P 7.5/1.0 
233 =o purplish gray 10.0P 5.5/1.0 
234 dark purplish gray 10.0P 3.5/1.0 
235 purplish black 10.0P 1.0/0.7* 
236 ~—s~vivid reddish purple LORP 4+4+-/19+% 
237 ~~ strong reddish purple LORP 4.5/11.1 
238 deep reddish purple LORP 2.8/10.3 
239 very deep reddish purple LORP 1.3/10.3* 
240 ~~ light reddish purple LORP 6.0/7.2 
241 moderate reddish purple LORP 4.5/7.2 
242. dark reddish purple LORP 2.7/5.1 
243 ~=~very dark reddish purple 1LORP 1.1/5.0* 
244 pale reddish purple LORP 6.0/4.1 
245 = grayish reddish purple LORP 4.5/4.1 
246 = brilliant purplish pink 4.0RP 7.9+/11+* 
247 = strong purplish pink L.0RP 7 /14.5* 
248 deep purplish pink t.0RP 6.0/12.0 
249 light purplish pink LORP 8.3/7.1* 
250 moderate purplish pink LORP 7.0/7.1 
251 = dark purplish pink 6.0RP 6.0/7.2 
252s pale purplish pink LORP 8.5/3.5* 
253) gravish purplish pink LORP 7.0/3.5 
254 ~=vivid purplish red 7.0RP 4+ /17++* 
255 = strong purplish red 7.0RP 4.5/11.9 
256 deep purplish red 7.ORP 2.8/11.0 
257 ~=very deep purplish red 7.0RP 1.3/9.0* 
258 moderate purplish red 7.0RP 4.5/9.0 
259 dark purplish red 7.0ORP 2.7/6.3 
260 very dark purplish red 7.ORP 1.2/4.9* 
261 light grayish purplish red 7.ORP 6.0/4.0 
262 grayish purplish red 7.O0RP 4.5/5.2 
263 white 3.0Y 9.25/0.06 
264 light gray 3.0Y 7.5/0.06 
265 medium gray 3.0Y 5.5/0.06 
266 dark gray N 3.5 
267 — black N 1.25 





accurately which ISCC—NBS designation to apply to 
it. If necessary, such borderline cases would have 
to be resolved by recourse to the method described 
in NBS Circular 553 [3]. 

These central notations have a number of appli- 
cations. Those published by Nickerson and Newhall 
|7] were used in designing the color-name charts used 
in the description of the colors of soils [5] and rocks 
6]. In both of these color charts. the Munsell 
Color Company painted samples representing the 
central notations of each of the color-name blocks 
used. Matching to one central color rather than | 
to several limit colors is a less exact but quicker 
method and was developed for field use where spee 
and ease of operation are more important than 
high accuracy. 


=“ —— — - 
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Viscosity of n-Hexadecane 
Robert C. Hardy 


Attempts to prepare high-purity n-hexadecane from commercial cetane by simple labo- 


ratory procedures were unsuccessful. 


Fractional distillation at reduced pressure, of ma- 


terial previously treated with silica gel, produced a few small fractions of about 99.5-mole- 


percent purity. 


The kinematic viscosities of the original material and of the fraction of 


highest purity, 99.63-mole percent, were found to be the same, 4.4635 centistokes at 20°. 


All other fractions had lower viscosities. 


The viscosity of NBS standard sample of n-hexa- 


decane, 99.94-mole percent, was found to be 4.4642 centistokes or 3.4540 centipoises at 20°C. 
n-Hexadecane is not likely to be suitable for use as a second calibration standard for viscom- 
etry until a simple, easily defined routine of purification is developed. 


1. Introduction 


Most viscometers must be calibrated by the use of 
one or more liquids of known viscosity. The vis- 
cosity of water at 20°C has been accurately deter- 
mined [1], and water at this temperature is widely 
accepted as the calibration standard for viscometry. 
However, each viscometer has practical usefulness 
for measuring only a limited range of viscosities, 
which in general, does not include the viscosity of 
water. Consequently, most instruments must be 
calibrated with liquids other than water. 

For about 30 years the Bureau has supplied a series 
of oils of known viscosity for use as viscometer cali- 
brating liquids. The American Petroleum Institute 
with the cooperation of the American Society for 
Testing Materials provides a similar but more limited 
service, used principally by the petroleum industry. 
These oils are multicomponent solutions whose prop- 
erties may change slightly over extended periods of 
time, depending on many factors including the con- 
ditions of storage. If deterioration due to age or 
contamination is suspected, fresh samples must be 
obtained from the supplier. There is no process by 
which their original viscosity can be reestablished or 
reaffirmed which does not require redetermination of 
the viscosity. 

This suggests the desirability of using pure chemi- 
cal compounds as calibration standards since if de- 
terioration occurred or was suspected, the liquid 
could be reestablished as a standard by any purifi- 
cation process which would return it to its original 
purity. Chemical compounds to be suitable for use 
as calibration standards should be liquid at normal 
room temperatures, readily available, easy to purify, 
reasonably stable. and  nonhygroscopic. They 
should have a low vapor pressure, and a surface ten- 
sion value in the range of about 25 to 30 dynes/em. 
The viscosity of the first of these compounds which 
would serve as a second calibration standard, prefer- 
ably should be several times that of water. Cetane 
(n-hexadecane) meets the requirements relative to 
viscosity, surface tension, and vapor pressure. Pre- 
sumably, it is also suitable with respect to stability 
and hygroscopicity. Also, a fairly pure grade (about 
95%) is commercially available as ASTM Reference 
Fuel Cetane. 

a 


' Figures in brackets indicate the literature references at the end of this papers 
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Although information developed prior to actual 
work on this project indicated that suitable purifi- 
cation of the commercially available cetane would be 
difficult, it was decided to proceed with the work 
since no more suitable commercial material was 
known to be available. 


2. Materials and Test Procedures 


Attempts to prepare highly purified n-hexadecaen 
from commercial cetane by recrystallization, by silica- 
gel filtration, and by the hydrocarbon-urea-clathrate 
complex, were unsuccessful. A first attempt to purify 
ASTM Reference Fuel Cetane by fractional distilla- 
tion at reduced pressure, using a 10-ft-long rectifying 
column | in. in diameter was also unsuccessful. Prac- 
tically all of the fractions had greenish-golden tints, 
apparently due to cracking of some of the impurities. 

A successful fractional distillation was made using 
a Podbielniak column at 75-mm-Hg pressure. The 
starting material was ASTM Reference Fuel Cetane 
which had been filtered through silica gel to remove 
some impurities which it was thought may have been 
responsible for pyrolysis in the first attempt. The 
fractions from this distillation were not colored and 
there was no evidence of cracking. Because of some 
flooding of the column and a drop in the readings of 
the differential refractometer on the last few frac- 
tions, the distillation was stopped when 22 fractions 
had been separated. These last fractions, 17 to 22, 
were returned to the still, the distillation resumed, 
and 8 more fractions separated. 

The freezing points of the Reference Fuel Cetane 
before and after silica-gel treatment, and of certain 
of the fractions, were determined and the purity of 
the samples calculated [2, 3]. The kinematic vis- 
cosities of these samples and of NBS standard sample 
No. 568-5S (99.94 + 0.04-mole-percent n-hexadecane) 
were determined at 20°C. Freezing point, viscosity, 
and density determinations were made on a composite 
of fractions 6 to 16, inclusive, and on a second, more 
recent (1957) lot of ASTM Reference Fuel Cetane 
(after silica-gel treatment). This lot was from the 
same supplier as the first lot used for the distillation. 
The appreciably higher purity of this material is 
understood to be the result of improvements in the 
production process. The densities of the composite 
sample and of the second lot of reference fuel (after 
treatment with silica gel) were determined at 20° 





and 25° C with a dilatometer (capacity about 50 cm*). 
The densities of the NBS standard sample, and of the 
composite sample also, were determined at room 
temperature with two small dilatometers of approxi- 
mately 2-cm* and 1l-cm* capacities. The values 
thus obtained were corrected to 20° C using the aver- 
age coefficient of thermal expansion derived from 
the measurements at 20° and 25° with the large 
dilatometer. 

Kinematic viscosity measurements were made 
with four Cannon master type U-tube viscometers, 
the same instruments used in determining the 
viscosities of the calibrating liquids supplied by NBS. 
Calibration of the instruments has been described 
elsewhere [4]. The instruments were used in pairs, 
and measurements were made with each instrument 
on each sample, except that only two instruments 
were used for the NBS standard sample. Also, after 
measuring three of the samples, one of the instru- 
ments was found to contain a small obstruction in the 
capillary, e. g., a lint fiber apparently introduced 
when charging the instrument. Observations in 
these three instances were discarded. The values 
reported for each sample and each instrument repre- 
sent the average of at least two observations of the 
flow time. The estimates of the standard deviations 
are based on the individual determined values, not 
the averages for each instrument. 

The temperature of the viscometer thermostat was 
measured with a platinum-resistance thermometer 
and a Mueller G—2 bridge which were calibrated 
prior to start of the viscosity measurements. The 
viscosity work was carried on somewhat inter- 
mittently over a period of 16 months. To detect 
any significant changes in the instruments or condi- 
tions, the instruments were thoroughly cleaned with 
chromic acid solution and the flow time for freshly 
distilled water at 20° C was determined for each in- 
strument before the first measurements were made 


and before resuming measurements after an inte, 
ruption. These tests did not reveal any significay 
changes. 

Times of flow were measured with electric gf 
clocks operating on a crystal-controlled 60-cps 
current. 


3. Results 


The results of the viscosity measurements togethe 
with data on freezing points and estimated purity gp 
presented in table 1. This table also includes resyly 
of the density determinations and estimates of the | 
standard deviations of the viscosity measurement 
for each sample. 

Considering only the samples of purity greater tha 
99.3-mole percent, the average value for kinematig 
viscosity at 20° C was found to be 4.4620 centistoke 





equivalent to +0.1 percent at the 95-percent cop. 
fidence level. 

The kinematic viscosity at 20° C of the sample gf 
| highest purity, i. e., the NBS sample with purity ¢ 
| 99.94+0.04-mole percent, was found to be 4.4642 «@ 
| with a standard deviation of 0.00099 es which js 

equivalent to £0.04 of 1 percent at the 95-percent 

confidence level. It is believed that the true valy 
| for n-hexadecane of absolute purity would fall withiy 
| these limits. Using the value 0.77370 g/cem* for the 
| density at 20° C, the absolute viscosity of n-hexa- 
decane was found to be 3.4540 centipoises (ep) 





4. Discussion 


shown in figure 1. Obviously for these samples there 
is no simple relation between kinematic viscosity and 
mole-percent purity as estimated from the freezing 
The failure to effect suitable purification by 


| A plot of purity versus kinemaiic viscosity i 
| 
| 
| 


point. 


(es) with a standard deviation of 0.0022 es which is} 











TABI E 1 Re sults of measurements 
Kinematic viscosity 
Density 4 bsolute 
Freezing | viscosity 
Material point Purity Instr. 1 Instr. 2 Instr.3 | Instr. 4 Mean 
Determinations at 20° C 
( Vole ‘ gcm cs cx cs cs ce s cp 
Reference fuel, lot 1 17. 262 13. OF 4. 4636 4. 4635 4. 4044 4. 4626 4. 4635 8. 210 
Reference fuel, lot 1 « 17. 308 04. 59 4. 4619 4. 4614 4.4615 4. 4606 4. 4614 5.3 
Distillate fraction 3 17. 377 4.44 4. 4578 4.4574 4. 4570 4. 4564 4.4572 8 
Distillate fraction 5 17. 856 17. 83 
Distillate fraction 6 4. 4505 4. 4591 4. 4585 4. 4590 2 
Distillate fraction 13 8. 090 09. 52 4. 4508 41. 4506 4. 4504 4. 4586 4. 4503 . 3 
Distillate fraction 16 18. O76 WO. 42 4.4599 4. 4505 4.4554 4. 4502 a 
Distillate fraction 21 18. 104 09. 63 4. 4636 4. 4636 3. 4633 4. 4635 1.9 
Distillate fraction 23 18. 100 09. 60 4. 4632 4. 4634 4. 4634 4. 4628 4. 4632 2.9 } 
Distillate fraction 24 18. O74 Wy. 41 4. 4625 4. 4027 4. 4625 4. 4622 4. 4625 2. 1 
Composite of fraction 6 to 16, 
inclusive P 18. O48 09, 22 0. 77370 4. 4610 +. 4606 4. 4601 4. 4505 4. 4602 &1 3. 4509 
Reference fuel, lot 2 17. 837 97. 69 77355 4. 4598 4. 4592 4. 4599 4. 4602 4. 4598 7.1 3. 4499 
N BS standard sample 568-58 18. 147 9. O4 7737 4. 4651 4. 4634 +. 4642 9.9 3. 4540 
Determinations at 25° C 
Composite of fraction 6 to 16, 
inclusive 18. O48 Ww. 22 0. 77021 3. 9724 3. 9722 3. 9723 3. 9717 3. 9721 3.8 3. 0504 
Reference fuel, lot 2 « 17. 837 97. 69 77013 3. 9719 3. 9715 3. 9718 3. 9710 3. 9716 3.4 3. 0586 


* Percolated through silica gel. 
> Calculated from freezing points and cryoscopic data given in reference [3] 
¢ NBS standard sample is the same as the A PI standard 
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Freezing point from reference [3] 
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KINEMATIC VISCOSITY AT 20°C,cSs 


Fieure 1. Kinematic viscosily versus purity. 


NBS Standard Sample; @, lot 1 of Reference Fuel; , lot 1 of Reference 
Fuel after silica gel treatment; , lot 2 of Reference Fuel after silica gel treatment, 
, distillate fractions 
recrvsiallization suggests the formation of solid solu- 
tions during the freezing. If this be irue, the caleu- 
lations of puriiy from the freezing poinis may be 
expecied to be in error, pariicularly if subsiantial 
amounts of impurities are present. Explanation of 
the scaiier of the poinis in figure 1 probably is to 
be found in this and the possibility that two or more 
impurities are present in different proportions in the 
different samples. 


TABLE 2. Absolute viscosity * of n-heradecane at 20° C 
Observer Viscosity Density Refer- 
ence 
cp gcm 
Ubbelohde and Agtha 3. 322 5 
f 3. 499 determined ln sono 
Evans, E. B o £60 emoethed 70. 7752 6 
Nederbragt and Boelhouwer b(3. 44 77385 7 
API Project 44 3.474 77342 Ss 
Schiessler et al 3. 453 7737 i) 
rhis work $454 re yi 
* Viscosity values reduced to basis consistent with the value 1.002 ep for 
water at 20° C 
>b Extrapolated from 25° (¢ Unreduced value. The basis of their cali- 
bration is not given 


Values for ihe viscosity of n-hexadecane at 20° C 
reported in the literature are presented in table 2. 
Various insiruments were used, bui, with one ex- 
ception, all were of the kinemaiie type and all were 
calibrated with waier. Evans [6] used mercury also 
as a calibraiing liquid. The values in table 2 have 
been reduced to the common basis of 1.002 cp for 
the viscosiiy of waier at 20° C insofar as information 
is available. Nederbragt and Boelhouwer [7] did 
not indicaie the calibration temperature or the value 
they used for the viscosity of waier. The API 
Project 44 [8] selected value is based on the work of 
Ubbelohde and Agiha (as reported by Engler and 
Hofer) [5], of Evans [6], and of Nederbragt and 
Boelhouwer [7], together with information from API 
Research Projects 42 and 44. This selected value 
appears to be high, but within its estimated uncer- 








tainty, i. e., 0.003 to 0.030 cp, it agrees with the 
value reported here and the value reported by 
Schiessler and co-workers [9]. 

In every instance, the value used for the density 
of the sample was a factor affeciing the value found 
for the absolute viscosity. Thus uncertainties in 
the density of the sample, though small, do contrib- 
ute to the uncertainty of the value found for absolute 
viscosity. In the present work, with the two small 
dilatometers, the densities of the NBS standard 
sample and of the composite sample were found to 
be 0.77344 g/cm’, which agrees with the API Project 
44 [8] selected value of 0.77342 g/em’. With the 
large dilatometer, the density of the composite sam- 
ple was found to be 0.77370 g/cm’ which agrees with 
the value 0.7737 g/em® reported by Schiessler and 
co-workers [9]. The values obtained with the small 
instruments were interpreted as indicating that the 
densities of the two samples were essentially the 
same. However the value obiained with the large 
instrument was believed to be more accurate and was 
considered to apply to both samples. 

lt may be noted that the values found for the 
kinematic viscosity of the original maierial and of 
several other samples are included in the range of 
values ascribed to samples of puriiy greater than 
99.3-mole percent. In faci, the average value for all 
the samples examined is 4.4611 es with a standard 
deviation of 0.00227 es and the total spread of values 
determined with the four instruments is only 0.18 of 
a percent. These figures cover material ranging in 
purity from 93.6- to 99.94-mole percent. However, 
since only two lots of ASTM Releccnes Fuel Cetane 
from the same source are involved, the data do not 
warrant the assumption that they could be applied 
to any sample of n-hexadecane of purity in this 
range. 


5. Conclusions 


Under present conditions, cetane or n-hexadecane 
is not satisfactory as a calibration standard for vis- 
cometry because of the difficulty of purification. 
Although many laboratories are equipped for purifi- 
cation work of this type, it is not to be expected 
that most laboratories interested in viscosity deter- 
minations would care to use a standard requiring 
such extensive preliminary treatment and testing. 
To attempt to overcome these difficulties by arrange- 
ments to supply certified samples of n-hexadecane 
from a central source would appear to have little if 
any practical advantage over the present system of 
supplying calibration oils by NBS and API. 

It is desirable that a calibration standard should 
be a material that can be purified by common labo- 
ratory techniques with a high degree of assurance 
that the purified material will have a definite vis- 
cosity. Also, it would be desirable, if not absolutely 
necessary, that the purity be verifiable by some test 
that is at least as sensitive as viscosity. At present, 
water seems unique in the ease with which these 
requirements are met. 

The objections to the use of water as a calibration 
standard are that its viscosity is low and its surface 
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tension is high, both in comparison with many organic 6. References Journc 
liquids, such as lubricating oils. Corrections can be 
made for the difference in surface tension between [1] J. F. Swindells, J. R. Coe, Jr., and T. B. Godfrey, J Re 
water used as a calibration standard and a liquid Research NBS 48, 1 (1952) RP2279. “7 
whose viscosity is to be measured, the uncertainty in | [2] A. R. Glasgow, Jr., A. J. Streiff, and F. D. Rossini, S | 
such corrections being well under 0.1 of a percent iCrnearch NBS 38, 355 (1945) RP1676; AgTm y 
such correc +4 a | ; 1D1015-55. 
when properly calculated [10] or determined experi- | [3] A. J. Streiff, A. R. Hulme. P. A. Cowie, N. C. Krousk 
mentally. Thus, if n-hexadecane were used as a and F. D. Rossini, Anal. Chem. 27, 411 (1955), 
calibration standard, the uncertainty that would | [I J. Pe ee Rar yp Cottington, 
arise due to the uncertain purity, and hence viscosity, 15} Knater and Hiséer, A ni Erdat 1. 53 hg 
of a given sample would be greater than that pres- | [6] E. B. Evans, Inst. Pet. Tech. 24, 38 (1938). 
ently resulting from the use of water as a calibrating | [7] G. W. Nederbragt and J. W. M. Boelhouwer, Physic, } 
liquid and evaluation of a surface-tension correction. . See ny righ nD ee rye | 
Similar difficulties probably would be encountered in 18} jin ony hems Pace Ay edeeceiees” oon pe Pl 
purifying other hydrocarbon liquids with viscosities compounds, Table 20c (1955), Table 20a (19586) TiO, 
equivalent to or higher than that of n-hexadecane. [9] R. W. Schiessler and F. C. Whitmore, Ind. and E; Cous 
Chem. 47, 1660 (1955) and private communication, Vri : 
{10} G. Barr, Proce. Phys. Soec., London, 58, 575 (1946). rie 
simil 
The author acknowledges the advice and coopera- 1 
tion of other members of the staff, and in particular - 
the work of A. B. Bestul and Vera Belcher in attempt- —_ 
ing purification by recrystallization; of Thomas W. — 
Mears and George S. Saines (present address Penn- of 4 
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freezing points, silica-gel treatments, and experiments = 
with the hydrocarbon-urea clathrate complex; the ay 
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System Calcia-Titania, Showing the Compound Ca.Ti;O, 
Robert S. Roth 


Vol. 61, No. 5, November 1958 Research Paper 2913 


The compound Ca,TisO% has been found to be a stable phase in the system CaO-TiOo. 


It melts incongruently at about 1755-4 


10° C, presumably to CaTiOs; plus liquid. 


This 


compound is shown in a revision of the phase diagram for the system CaQ-TiOnz. 


1. Introduction 


Phase equilibrium diagrams for the system CaQ- 
TiO, were published almost simultaneously by 
Coughanour, Roth, and DeProsse [1]' and by De- 
Vries, Roy, and Osborn [2]. The diagrams are very 
similar, both showing the existence of two com- 
pounds, CaTiO; and Ca;Ti,O;. DeVries et al. [2], 
discussed the reactions between the two compounds 
and, on the basis of some changes in the X-ray pat- 
terns of the mixtures, concluded that a small amount 
of solid solution existed in both compounds. They 
published an X-ray pattern of Ca;Ti,0; which was 
quite similar to that reported by Fisk [3] for the same 
compound. Coughanour, Roth, and DeProsse [1] 
did not find solid solution between the two com- 
pounds but reported an X-ray pattern for Ca;Ti,O, 
which contained many more lines than those of the 
other publications. 

Ruddlesden and Popper [4] studied the compound 
Sr,Ti,O; and have indicated it to be body centered 
tetragonal with a=3.90 A and c=20.38 A. They 
have also found a compound Sr,TiO, [5] which con- 
tains alternate layers of perovskite and SrO, and 
have proposed that Sr,Ti,O; contains double perov- 
skite layers interleaved with one SrO layer. 


2. X-Ray Diffraction Data 


The X-ray pattern previously published by Coug- 
hanour, Roth, and DeProsse [1] for Ca,Ti,O0; can be 
partially accounted for on the basis of a body cen- 
tered tetragonal cell with a=3.832 A and e=19.505 
A. The pattern published by DeVries, Roy, and 
Osborn [2], which contains fewer lines, can be almost 
completely accounted for by a body centered tetrag- 
onal unit cell with a=3.8 A and c=19.5 A. These 
values correspond to double perovskite layers inter- 
leaved with CaO layers. 

In the SrO-TiO, system Ruddlesden and Popper 
(4] apap the possible existence of phases with 
triple or quadruple perovskite layers interleaved 
with SrO layers. They found some evidence for a 
compound Sr,TizOy, which is also body centered 
tetragonal, with a=3.90 A and c=28.1 <A, and 
contains triple perovskite layers. 

Many of the extra lines in the pattern of Coug- 
hanour, Roth, and DeProsse [1]? can be accounted 


_— 


' Figures in brackets indicate literature references at the end of this paper. 

? This pattern is represented in the ASTM index of X-ray diffraction powder 
patterns by card No. 6-0698. As this pattern is a mixture of two compounds it 
should be eliminated and replaced by the two patterns listed in this paper. 





for by assuming a mixture of two body centered 
tetragonal compounds, Ca;Ti,O; and Ca,Ti;Oj, the 
latter having a c-axis of 27.147 A. In the present 
study the two compounds have been isolated and the 
partially indexed X-ray patterns of both compounds 
are shown in tables 1 and 2. Although these two 


TABLE 1. X-ray diffraction powder pattern for the compound 
Ca; TivO; (CuK; radiation) 
. } 
DeVries, 
Fisk Roy, Present work 
Osborn hkl 
d d dobs I/TIo 1/d2ovs 1/d2 eal e | 
| 
9.77 9.73 > 100 0. 0106 0. 0105 002 
(4. 88) 4. 887 4. 873 b 100 . 0421 0421 | 004 
3. 762 3. 762 . . 0707 . 0707 101 
3. 302 3 . 0907 0918 | 103 
: 0946 | 006 
2. 730 2. 733 100 . 1339 . 1338 105 
2. 70 2.710 2. 708 66 . 1364 . 1362 110 
. 1467 | 112 
2. 51 2. 481 2. 478 7 . 1629 eee. 
2. 439 b5 . 1681 . 1682 008 
- " 1783 | 3914 
> 19 «|f (2-254 2. 255 11 . 1967 1969 | 107 | 
oor 2. 081 23 -2308 | .2308 | 116 | 
1.93 1. 9509 1. 9505 b 100 . 2629 . 2629 | 0-0-10 | 
1.91 1.9141 1. 9159 27 . 2724 . 2724 200 | 
1, 8860 1. 8871 14 . 2808 2810 | 109 
. 2829 | 202 =| 
. 3044 118 
1. 7836 6 3143 3145 
1. 7363 3317 ; 
1, 7080 3 3428 . 3432 211 
1. 6581 4 3637 . 3642 213 
i . 3671 206 
1. 6258 b 37 3783 . 3785 0-0-12 
Be : . 3862 1-0-11 
1.57 |f 1.5826 1, 5831 38 . 3990 3991 | 1-1-10 
—" tf ia 1. 5690 27 . 4062 4062 | 215 
0. 4406 208 
1, 4598 4 0.4693 | . 4693 217 
. .5123 | 10-13 
: : 5147 | 1-1-12 
1. 3912 1. 3933 b 32 .5151 | .5152 | 0-0-14 | 
1. 358 1. 3662 1. 3667 14 5354 | .5353 | 20-10 
1. 3539 10 . 5455 5448 | 220 
1. 3423 1, 3444 5 . 5533 5534 | 219 | 
, ‘ ; . 5554 222 =| 
. 5869 224 
6156 | 301 | 
666 | 303 
. 6395 226 
1, 2304 5 . 6510 6509 2-0-12 
: 6514 | 1114 | 
| 
anes oi : lf .6586 | 21-11 
1.2303 | 1.2318 5 ae eet 
1. 2191 b15 . 6729 . 67: 0-0-16 
. . 6787 305 
1, 210 1.2119 1, 2109 4 6819 6811 | 310 
/ ‘ . . 6016 312 
Kin donee 7131 | = 28 | 
SS ese: 7231 | 314 | 
eS 307 =| 
sa Ce Se 316 | 
Ba. a. Cos J 


See footnotes at end of table. 
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TABLE 1. Y-ray diffraction powder pattern for the compound TaBLeE 2. X-ray diffraction powder pattern ‘for the compound 
Ca;Ti,O; (CuK; radiation) —Continued CayTizO;) (CuKa; radiation) 
dobs I/To 1/d2obs 1/d2ent ® hkl o 
DeVries, ’ | 
Fisk Roy, Present work | 
Osborn hkl 13. 5 6 0. 0055 0. 0054 002 
6. 78 7 0217 0217 ood | 
4. 525 11 O48S O4S8Y 006 | 
d d dove T/T 1/2 obs 1/d2-a1* 3. 705 i) O64 O66 101 bh 
a 7 . OR805 103 "3 
7876 2.0-14 | OS6S OOS 
1. 1234 2 7¥24 | 3. 337 5 OSUS 
1. 1124 4 8081 8077 2-2-10 1022 105 c 
0. 8091 1-1-16 2. 725 100 1347 1348 107 
. S250 309 2. 718 100 1353 1357 0-0- 10 
1. 0995 4 8272 8277 1-0-17 | } 
S43 318 2. 706 16 1365 1365 110 
R516 0-0-18 | 1420) 112 
1582 114 4 
SAND) 321 | 2. 422 s 1705 . 
1. 0559 4 SOY } 2. 367 4 1785 . 1782 109 
gO90 323 
W233 2-2-12 Is54 116 
Y310 3-0-11 2 283 ; 1919 e 
1u4 0-0-12 " 
9319 11S 2.212 { 2044 
1. 0292 6 O44 W439 3-1-0 2.116 10 2234 2234 11S 
M53 24-16 
1. 021 1, 0251 { 9516 O511 325 2.075 3 2323 2325 1-0-11 Fiat 
ONTS 1-1-18 | 1. 939 33 Jt) 260 0-0-14 fr¢ 
ie 1.915 65 2726 ac 1-00 : 
1.0142 $27 ‘ 200) 
0. G9 168 10 1. O168 1. O1l70 1-0-19 785 a2 
97534 71 1.0512 1.0514 0-0-20 | 1.874 3 2848 
1. 0572 3-0-1383 
1. 0506 3-1-12 2048 214 
1. 832 7 2979 2076 1-0-13 
1. 0600 22-14 1. 764 6 3214 3219 206) 
GAT 2 1. O4O4 1. OR07 1) 
{ 1. 0083 329 j 
95359 3 1. 0997 1. 1002 2-1-17 0. 3319 1-1-12 
| 1. 1002 102 1.713 { (). 3409 3427 211 
1. 1241 2-0-18 
1. 1317 hoe 3474 0-0-16 
} 3535 213 
1. 1604 411 1. 4 5 sH12 3504 208 
1. 1814 $13 3736 1-0-1 
1. 1843 1) $753 21 
0. 91768 7 1. 1875 1. 1876 11-20 “2 | 
2 1. 1962 $-1-14 | 1. 576 15 12s HI25 1-1-14 
1078 217 
1. 2034 3-2-11 | 1. 65 i) 1085 4088 24)-10 
. 91125 1. 2043 1. 244 0-15 | 1306, 0-0-18 
1. 217 2.16 4512 219 
l 415 | 
. 90346 6 1. 2251 1. 2259 330 04 1-0-17 
| Hs5 2-0-12 
. 90236 1 1. 2281 1-0-21 1X34 1-1-16 
332 1. 40 1 (57 WSS 2-111 
6 10S 1. 362 21 S388 A390 24)-14 
334 
SS673 643 1. 2718 0-0-22 5428 0-0-20 
| 1. 351 +) ATT 461 22) 
. 1. 2866 4] 1. 326 1 
SSOSS s 1, 2806 1. 2804 2.1.19 1. 323 
1. 3205 $36 1. 256 ' a 
SHOLS { 1. 3238 1. 3238 2.4)-20 
a 1. 3296 3-2-13 1. 224 1 
1.213 12 Ps 
S5US2 5 1, 3526 1, 3425 4.0.10 1. 109 7 
1. 3540 3-1-16 1. ONT 3 
1. 3621 $20 1 O52 
{ 1.3707 419 
SSS $ 1. 3715 1. 3726 30-17 | 1.044 3 
| 1.3726 122 1. 028 5 
1. 3919 2-1-20 | 1. 026 +] 
| 1. 023 1 
1. 3941 338 
1 2065 9.2.18 
1. 4042 124 ‘ , 
. J The 1/d? values have been calculated on the basis of a tetragonal unit cell | 
ean wes | 1.8 | 611-2 | with a=3.827A and ¢=27.147A 
: 1. 4564 126 >» Due to the very large unit cell, no attempt was made to index the rest of this 
1. 4586 1.0-23 | Pattern : 
1. 4682 £-0-12 | patterns can be almost completely accounied for on 
1. 4758 4-1-1] the basis of two tetragonal unit cells, there are still 
o aioaa - anna aoe | ore some very small peaks which cannot be indexed on 
size 22 Lc | Lae | 002 this basis. It seems possible thai these caleia- s 
303 428 : : . 
tiiania compounds are not completely isosiructural 
500a3 ; aig | cae | saun with the corresponding strontia-titania compounds, 
oe | son li is possible that these compounds have lower 
O62 2.2.20 - .. vim , 
symmetry, just as CaTiO; has a lower symmetry Fi 
than SrTiQsg. F 
* The 1/d? values have been calculated on the basis of a tetragonal unit cell De\ ries, Rov, and Osborn [2] have published parts 
with a=3.832A, c=19.505A | _ lle . : > “a d 
» These (001) peaks are probably abnormally strong due to preferred orientation. ! of X-ray diffraction patterns of a series of mixtures 8 
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COMPOUng 


from Ca;Ti,O; to CaTiO,; after melting and quench- 


\\ N A ing. Part a (of fig. 3, ref. [2] reproduced here as fig. 1) 
\ A.J \ 2m / \y~_. | represents a portion of the pattern of Ca,Ti,O,; and 
: part e that of CaTiO;. Figure 2 shows essentially 


= 4 





the same portions of the X-ray patterns of the two 
compounds CagTi,O; and Ca,Ti;O;o, as found in the 
present work. Although all of the patterns in 
figure 1 are very poorly resolved it can be seen in 
part b that there are three peaks in the region of 46 
to 48°20. The composition of this mixture is very 
close to 4CaO:3TiO, and this pattern represents the 
compound Ca Ti,O;9 with a small amount of CasTi,O,. 


} / . . . 
“aN e \ AN Parts c and d are mixtures of CaTiO,; and Ca,Ti,O,9 
= Sf 
/ “ a and both show the center of the three peaks at 47° 286, 
\ - while this peak is eliminated in part a, Ca,Ti,O;. 
\ \ " Py These published X-ray diagrams support the theory 


of an intermediate compound (Ca,Ti;Oj.) much 

aE " better than the theory of solid solution proposed by 

the authors of the published patterns. 

Fieure |. X-ray diffractometer patte rns of a series of mixtures 
from CasTizO; to CaTiO; (after melting and quenching). 


3. Experimental Procedures 


Taken from DeVries, Roy, and Osborn [2] 


1=51.29/48.71 wt. % Ca0/TiO2g (CagTinO:) 
b=48.65/51.35 wt. ©, CaO/TiO2 

‘ 46.00/54.00 wt. ©; CaO/TiOg 
d=48.60/56.40 wt. ©, CaO/TiOeg 
e=41.24/58.76 wt. © CaO/TiO, (CaTiO 





It might be thought possible that Ca,TigOjo is only 
a metastable compound having no true equilibrium 
position in the CaO-TiO, phase diagram. If such 
were the case, an equi-molar mixture of preformed 
CaTiO, and Ca,Ti,O; would not form the Ca,Ti,O;, 
phase upon heating. Such an equi-molar mixture 
was prepared, using the same raw materials as in the 
previous study (Coughanour, Roth, and DeProsse 
[1]). The specimen was heated in a conventional 
Pt-wound quench furnace to 1,650° C, held for 2 hr, 
and quenched in air. The compound CajTi,O,) was 
the only product identified in this specimen. It 
must therefore be concluded that the compound 
Ca,TisO,) is the equilibrium product of the compo- 
' sition 4CaO:3TiO,, at least at 1,650° C. 

From the data of Coughanour, Roth, and De- 
Prosse [1] as well as that of DeVries, Roy, and 
Osborn [2], it can be seen that the compound Ca,Ti,- 
| Os can net melt congruently. It must either dis- 
~V\ J “/u\,. | sociate before melting or else melt incongruently at 

a temperature somewhat higher than the incon- 
; aa ae ; a a | gruent melting point of Ca,Ti,O;; given in both 
papers as 1,750° C, 

Specimens with a calcia-titania ratio of 4:3 were 
pressed into pellets, calcined in a conventional Pt- 
wound furnace at 1,250° C, ground, repressed, and 

| _ reheated to 1,600° C and slow-cooled. hey showed 

only the compound CayTi,Oj. Specimens were then 

| unit eel ground and placed in a Pt pln 2, heated in the 
~~ Pt-wound quench furnace, and quenched into water 
from temperatures of 1,500° C, 1,550° C, 1,600° C 








| 
for on : | and 1,650° C. All of these quenched specimens still 
re still alll . _ showed only the compound Ca,Ti,O;o. Specimens 
ced on , e \ \ | were then heated and quenched in a graphite 
calcia- J av susceptor induction furnace, from temperatures of 
ctural < ~~ eS eet lr «| «1,700° C and 1,725° C. These two specimens also 
yunds, | showed only the compound CayTi,Oj. One more 
lower —— ee ee 7 : | specimen was placed in a smooth Pt tube, so that 
metry | the temperature could be read more precisely with 


, » ) = — ne ° 8 — A =» = 
Ficurr 2. X-ray diffraction patterns, from the present work, | gy optical pyrometer, heated to 1,745° C and 
of the two compounds Ca;Ti,O; (above) and CayTizOw ‘ 


; ie ee ae 
| parts below) covering essentially the same part of the diffraction | quenched. This specimen still show ed only ¢ ay Tis- 
xtures spectra as figure 1. | Oy. As the reported solidus for CasTi,O, is 1,750 
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Figure 3. Revised phase equilibrium diagram of the system 
CaO-TiQ»2, showing the compound Ca,Ti2O; melting incongru- 
ently at 1,740° C and the compound CayTizQOy melting 
incongruently at 1,755° C. 

(After Coughanour, Roth, and DeProsse [1]). 


@ Solidus temperatures determined in the present study. 


C, it must be assumed that Ca,Ti,O,, either does not 
dissociate to two solids, or that such dissociation 
cannot be quenched. 

In order to check the solidus temperatures of the 
two compounds Ca;Ti,O; and Ca,TizO, fragments 
of the two materials, previously heated to 1,650° C 
for 2 hr, were set side by side on an iridium button in 
the graphite induction furnace. The Ca,Ti,O, spec- 
imen definitely began melting before the CayTizOjo 
specimen. However, the fluid liquid quickly spread 
over the iridium button, obscuring the results of 
melting of the second specimen. As it appeared 
that the iridium button was hotter than the inside 
of the graphite crucible, a new experiment was per- 
formed utilizing only an iridium crucible with an 
iridium button inside. Using a flat button and a 
cover on the crucible and sighting through a hole in 
the cover, perfect black-body conditions were ob- 
tained and the specimen could not be seen. By 
utilizing a concave button and no cover, the speci- 
men could be watched during the whole experiment. 
The temperature of the solid specimen at beginning 
of melting appeared 25° to 30°C cooler than the 
true temperature, but the liquid formed appeared 
to show the true temperature. The solidus temper- 
ature for Ca;Ti,O; was found to be 1,740°C and 
that of CayTi,O,) was 1,755° C. These results were 
rechecked by again heating fragments in the iridium 
crucible using a flat button and the crucible cover. 
The Ca;Ti,O,; specimen when heated to 1,740° C 








showed a rounding of the corners and a small umMoun 
of crystal growth on the iridium button. The 
Ca,yTi,0,9 specimen showed no sign of melting 
this temperature but did show some melting Whep 
heated to 1,755° C. 

These measurements have been used to revise the 
phase diagram of the binary system CaO-TiO, fipy 
proposed by Coughanour, Roth, and DeProsse nt 
The revised diagram is shown in figure 3. 





4. Summary 


The compound Ca,Ti;Ojo has been shown to exig | 
as a stable phase in the system CaO-TiO,. Bot) 
compounds are pseudo-tetragonal with only a fey 
weak lines indicating a departure from tetrago. 
nality. The c-axis=19.505 A for Ca;Ti,O, ang 
27.147 A for CaTizOy. The dimensions indicate 
that Ca,Ti,O; has double perovskite layers inter. 
leaved with CaO, while CasTizOjo has triple peroy- 
skite layers interleaved with CaO. The new com. 
pound, Ca,TisOy, has been found to melt ineop. 
gruently at about 1,755° C, while the incongruent 
melting point of CasTi,O; has been shown to be 
about 1,740° C. The revised phase diagram for the 
system CaQ-—TiO, is presented. 


The assistance of Charles R. Drew in preparing 
and heating the specimens and of Sylvanus F. 
Holley for operating the induction furnace is grate. | 
fully acknowledged. 
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} The two previously published diagrams of Coughanour et al. [1] and DeVries 
et al. [2] disagree on the temperatures of the melting point of CaTiOs (1,915° C 
and 1,970° C) and on the temperature of the eutectic between CaO and CasTig0; 
(1,725° C and 1,695° C). No attempt has been made in the present study to 
reconcile these two discrepancies 
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